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GEOLOGY 


TECTONICS OF A METAGABBRO LACCOLITH IN THE 
ADIRONDACK MTS., AND ITS SIGNIFICANCE IN DETERMINING 
TOP AND BOTTOM OF A METAMORPHIC SERIES 


BY 


D. DE WAARD 


(Communicated at the meeting of February 25, 1961) 


ABSTRACT 


A series of high-grade metamorphic, metasedimentary, Precambrian rocks in 
the southwestern Adirondack Mountains is folded to a closed, overturned ‘downfold’. 
The determination of whether the downfold represents a syncline or the plunging 
front of an anticline depends on the determination of the direction of top and bottom 
in the folded rock series. One small metagabbro mass has the plano-convex shape 
of a laccolith; its shape enables determination of top and bottom of the intruded 
series prior to the folding and indicates that the nature of the downfold is a syncline. 

A tectonic analysis is given of the laccolith and its contact relations. Using the 
axial symmetry of a laccolithic mass a composite section is constructed normal to 
the foliation to illustrate the possible shape and the position of the gabbro laccolith 
in the metamorphic series. 


Introduction 


The object of this note is to demonstrate the use of a laccolithic structure 
in the determination of top and bottom in a folded series of metamorphic 
rocks. The series in question consists of a succession of sheets of various 
gneisses, marbles and amphibolites, with sills and small bodies of 
metadiabase and metagabbro. The series is folded, almost isoclinally, 
into an overturned fold with an axial plane dipping about 35° N. Detailed 
mapping revealed a structure which can be described as a ‘“‘down fold” 
or a “synform”’ 1); depending on top and bottom direction in the sheets 
the structure is either a syncline, the simplest interpretation, or the 
plunging front of an anticiline, a not uncommon type of structure in 
Alpine tectonics. Too little is known of the geologic relations in this 
region to indicate which of the two types of structure might be expected ; 
the synclinal or anticlinal nature of the down fold can be determined 
only by ascertaining the top and bottom directions in the rock series. 

In metamorphic rocks, features which indicate top and bottom of a 
series are mainly limited to the inherited structures of sedimentary and 


1) J. G. Ramsay, Superimposed folding at Loch Monar, Inverness-shire and 
Ross-shire, Q. J. Geol. Soc. London, 113, 1958, 271-307. 
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igneous origin. SHROCK !) lists a number of such palimpsest structures, 
none of which were observed in this area of high-grade metamorphism. 
Not mentioned by SHrock is the shape of an intrusive body as a possible 
indication for top and bottom of the intruded series. 

Most of the intrusive bodies in the area in question are diabase sills; 
metagabbro occurs in thicker sills, in lenticular structures, and in one 
body which has a short and thick, plano-convex shape with concordant 
contacts, located in the north limb of the down fold. Absence of orientation 
in the blastophytie texture of the metadiabase and metagabbro, except 
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CHARNOCKITE META GABBSRO 
GARNET-B/IOT/TE GNEISS METAD/ABASE 


MARBLE AND AMPHI/IBOLITE DOMINANT, INCLUDES OTHER METASED/IMENTS 


Fig. 1p Geologic map of the metagabbro laccolith. The topography of this area 
is shown in map, fig. 2. Strike and dip symbols indicate the orientation of foliation 
planes; arrows refer to linear mineral parallelism. 


1 ‘ wae x : 
) R. R. Surocx, Sequence in layered rocks, McGraw-Hill, New York, 1948 
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locally such as along contacts, indicates that little internal deformation 
took place in the rock during metamorphism, and that the original shape 
of the intrusions largely has been preserved. Assuming (1) that the 
plano-convex gabbro intrusion is a true laccolith, and (2) that the convex 
side of a laccolith always represents the top of the intrusion, then the 
position of this structure indicates top and bottom directions in the 
series which, in turn, enables the determination of the synclinal or 
anticlinal nature of the down fold. 


The laccolith 

The gabbro laccolith is located at 43° 38’ 29” N, 74° 41’ 25” W, in the 
Adirondack Mountains, West Canada Lakes quadrangle, New York. 
Topographically, the outcrop of gabbro stands out conspicuously as a 
round hill amidst a strongly linear topography controlled by the north- 
ward dipping series of gneisses. 


METERS 100 


Fig. 2. Topographic expression of the metagabbro laccolith (shaded eee), Radial 
lines indicate the location of sections used to construct composite section, fig. 6; 
section A—B is shown in fig. 3. 
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A detailed geologic map (fig. 1) illustrates the structural relations of 
the rocks in this area. The foliation planes of the rocks NW of the gabbro 
body appear little or not affected by the presence of the intrusion; 
anywhere else the foliation tends to conform to the orientation of the 
gabbro contact, and to diverge from the general trend of the foliation in 
the area. Except for the NW contact the foliation planes of the gneisses 
near the gabbro border dip between 25° and 35° towards the center of 
the intrusion, which suggests the shape of the gabbro mass approximates 
a segment of a sphere. The NW contact is straight and parallel to the 
foliation of the gneiss, which indicates a flat contact plane dipping about 
35° N. Section, fig. 3, illustrates the shape of the remaining segment of 
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Wig. 3. Section A-B across the laceolith. Orientation and loeation of the syminetry 
axis was determined to construct a composite section (fig. 6) perpendicular to this 


section and containing the axis. 


this plano-convex gabbro mass as based upon the contact relations in 
the outcrop; a possible reconstruction of the removed part of the laccolith 
is shown in broken lines. 

Intrusion occurred prior to the period of metamorphism in the area, 
and, most likely, prior to the folding of the intruded beds. The shape of 
a laccolith is most simply explained by a doming-up of overlying beds 
during the intrusion at a time that the intruded strata were in normal, 
probably flat-lying position !). Since the convex side of a laccolith is in 
contact with the overlying beds, the present position of the metagabbro 


1) G. KX. GitBeRT, Report on the geology of the Henry Mountains, U.S. Geo- 
graphical and Geological Survey of the Rocky Mountains region, Washington, 1877 
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mass 18 overturned in the northern limb of the down fold, and the 
downfold, therefore, is a syncline. 


Tectonic analysis 


Diagrams, fig. 4 and 5, illustrate the effect of the convex contact of the 
laccolith on the orientation of adjacent foliation planes. Poles of unaffected 
foliation planes (dots) are concentrated around an average of 253-36; 
poles of foliation planes from near the convex contact (circles) are scattered, 
towards and around the center of the diagram. Similarly, unaffected 
lineations (crosses) are grouped together while those found near the 
gabbro contact (asterisks) deviate to various degrees from the average. 

Diagram, fig. 5, attempts to elucidate the deviation of linear and 
planar structures as shown in fig. 4. From the average orientation of 
unaffected foliation planes, S, rotation occurs radially; small circles show 
the degree of rotation, and arrows indicate the sense of rotation of the 
poles of foliation planes. The axes of rotation are situated in the great 
circle normal to 8S, which is the average foliation plane, and to the flat 
base of the laccolith. Lineations, being lines in foliation planes, rotate 
with the foliation from the average, B, in paths along small circles, some 
of which are shown and marked a, b, c, d, e, and f, corresponding to 
the similarly marked planes and axes of rotation of the foliation. A 
dashed line indicates the effect of 60° rotation of B along various paths. 

The dispersion of poles as shown in fig. 4 corresponds to the shaded 
area in diagram, fig. 5. The restriction of the rotation to this area is, 
of course, due to the fact that data are obtained from the present outcrop 
which is only one section of the laccolith. Other sections would show 
different degrees and other directions of rotation. In the present section 
a maximum of rotation is observed of 55° (c, d, and e) about rotation 
axes oriented between NE-SW and E-W. There is a corresponding 


rotation of linear structures. 


Composite section 


Construction of composite sections is based upon the symmetry of the 
structure, and on the assumption that this symmetry persists throughout 
the structure. Composite sections serve to give a more complete picture 
of the structure by projecting the observed surface data on the plane of 
the section in accordance with the symmetry properties of the structure 
in question. The symmetry of the gabbro body as observed at the surface 
in the shape of the contact and contact relations point to the axial 
symmetry of a laccolith, the principle of which is shown in diagram, 
fig. 5. Transfer of surface data to the plane of the composite section is 
based upon the symmetry property that all sections which contain the 
axis of symmetry are identical. Since contact relations are concordant 
this also holds largely for the adjacent wall rock. 
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Fig. 4. Pole diagram illustrating the observed effect of the proximity of the convex 

contact on the orientation of foliation planes and mineral lineations in the area. 

Poles of unaffected foliation planes (dots) and lineations (crosses) are found in 

groups; poles of foliation planes from localities near the contact (circles), and the 

lineations from those localities (asterisks) are dispersed by rotation. (Equal area 
projection on the lower hemisphere). 


Fig. 5. Analysis of diagram, fig. 4. Foliation planes near the convex contact of 
the laccolith are rotated from the unaffected average, S, about axes of rotation 
situated in the plane normal to 8, along paths radiating from §S, in a sense indicated 
by arrows, and to the amount shown by small circles. Similarly, lineations in the 
foliation plane rotate along paths radiating from the unaffected average, B; an 
amount of 60° is indicated by the intersecting dashed line. Due to the restriction 
of data to the present surface (which is one section through the laccolith), the sense 
and the degree of rotation which can be observed is restricted to the shaded area 
in the diagram, and is indicated by full-drawn lines and arrows. The observed maxi- 
mum deviation of 55° means that a maximum angle of 55° is found between a contact 
on the convex side and the floor of the laccolith, 
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Fig. 6. Composite section across the laccolith. Surface data along sections (locations 

shown in map, fig. 2) which contain the axis of symmetry are projected by rotation 

on the plane of the composite section which is perpendicular to the average foliation 

plane with a strike parallel to that of the foliation. The section is inverted to demon- 

strate the constructed shape and the contact relations of the laccolith in the position 
of intrusion. 


The location of the symmetry axis was found by using the orientations 
of foliation planes near the contacts to determine the radii and centers 
of circular sections through the laccolith at right angles to the axis. The 
surface geology of each section shown by lines radiating from the symmetry 
axis in map, fig. 2, are combined in fig. 6, and sketched in full to a 
composite section which, in restored position, illustrates the possible 
shape of the laccolith, and its position in the series of metamorphic rocks. 

The section demonstrates the different plastic behavior of the various 
rock kinds during metamorphism and deformation. The  laccolith, 
apparently, stood out rigidly during deformation, pushing away the more 
mobile material in the layers above it. Particularly the biotite-garnet 
gneiss appears to have the tendency to be squeezed out near the top of 
the laccolith and to accumulate alongside the gabbro mass. 


Geologic history 
The sequence of events in the area can be summarized as follows: 
1. Sedimentation of alternating beds of various chemical composition 
corresponding largely to that of charnockite, biotite-garnet gneiss, marble, 


and amphibolite. 
2. Intrusion during the period of sedimentation of the gabbro laccolith, 


var 


lenticular gabbro masses, and diabase sills. 
3. Period of metamorphism and orogeny. Metamorphism of sediments 
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to charnockite, biotite-garnet gneiss, marble, and amphibolite, and 
metamorphism of the basic intrusives to metagabbro and metadiabase. 
Synmetamorphic folding of the rock series; development of foliation 
parallel to the stratification, formation of a mineral lineation parallel 
to the fold axis of the syncline. A rigid structural behavior of larger 
metagabbro and metadiabase masses; development of mineral orientation 
in smaller masses, and locally along shear zones and contacts of larger 
intrusions. During folding the laccolith rotated 144° with the northern 
limb of the syncline to an overturned position. 
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FROSTING OF QUARTZ GRAINS 
BY 


PH. H. KUENEN 1) ann W. G. PERDOK 2) 


(Communicated at the meeting of March 25, 1961) 


Abstract 


It is shown that the frosted appearance of many quartz grains is not 
produced by circular impact markings of a few microns, but by diffraction 
of light on minute irregularities. Impact cones on rounded quartz grains 
are much larger, nearly one tenth of the grain diameter. Wind action 
gives coarse frosting on protruding parts of a grain, but this is subdued 
by the action of dew or under water. Aqueous transportation can also 
frost protruding parts. Reaction solution (etching) tends to produce 
frosting, as does alternate wetting and drying by dew. Non-reaction 
solution (physical solution) tends to round and de-frost, and immersion 
in a saturated solution can also cause de-frosting by regrowth. 


Introduction 

Most grains of quartz sand are frosted, if only slightly, and in the 
desert the surface of the particles is usually as opaque as ground glass. 
Most authors have attributed this surface finish to mechanical attack, 
by scour under water and more especially to eolian impacting. However, 
SoRBy (1880), Russet (1939), TRicarr (1958) and others have quoted 
cases in which chemical corrosion must have been responsible. Polishing 
is generally believed to result from mechanical rubbing under water. 
In earlier papers one of us has pointed out that experimental frosting 
in a wind tunnel produces a coarser finish than natural frost and that 
it is possible to polish sand by rolling in dry condition, but that under 
variable conditions of rolling under water grains remain opaque (IKUENEN, 
1959, 1960). The conclusion was drawn that both frosting and de-frosting 
must be attributed mainly to chemical causes. 

This surmise was confirmed by noting the distribution of the surface 
finish over irregular or very small particles. Experimental mechanical 
polish or frost are both restricted to convexities of the surface, while 
hollow and even flat parts are left unchanged. In contrast, natural frost 
covers the entire grain uniformly, even in deep depressions. Besides, 
grains in the process of becoming de-frosted attain the same degree of 
transparency in concavities as on bulges. A weighty argument was found 


1) Geological Institute, University of Groningen. 
2) Institute for Crystal Physics, University of Groningen. 
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in the frosted surface of very small, entirely angular particles of 
1/59 millimetre occurring in many desert sands, for these have not under- 
gone any detectable impact-abrasion. Alternate wetting by dew and 
drying was suggested as the cause of desert frosting, solution and 
deposition by under- and oversaturated ground water or river water can 
account for aqueous frosting as well as polishing. 

In the present paper the nature of frosting is examined and its causes 
are dealt with briefly. A more extensive treatment will be given elsewhere. 


Optical nature of frosting 


CAILLEUX (1942), led by his conviction that frost is due to wind action, 
claimed that under powerful magnification the surface is seen to consist 
of minute impact scars of circular or linear shape of 5 to 10 microns 
diameter. For several reasons this view cannot be upheld. 

In the first place slight etching of quartz by HF produces a frost that 
presents the same picture under the microscope. In the second place the 
minute angular desert chips likewise show this circular pattern. In the 
third place in deep pits the frost is also visible as a circular marking. 
In the fourth place electron micrographs of frosted grains show complicated 
irregularities of 1 micron or less, but fail to reveal any circular shapes of 
a few microns diameter. In the fifth place obvious impact cones are 
found on many rounded river grains (probably with an earlier eolian 
history), but these are much larger and their size is related to the 
dimensions of the particle, just as with pebbles, namely roughly one 
tenth of the diameter. 

The true cause of the circular pattern shown by frosted grains is the 
diffraction of the light by small surface irregularaties. It becomes even 
more marked when the surface is coated with a strongly reflecting bright 
film of aluminium. 


Two kinds of solution 


To understand the non-mechanical origin of frosting and de-frosting 
it is necessary to realize that two entirely different processes are designated 
by the term ‘solution’. One can be called reaction solution (etching or 
corrosion). It consists of a chemical reaction between a solid and a liquid, 
like the solution of CaCO 3 in HCl or quartz in HF. Evaporation will not 
cause the original solid to be reformed. The other process can be indicated 
by non-reaction solution (physical solution). It does not involve a chemical 
reaction and evaporation will cause the original solid to crystallize again 
(provided no complications are introduced), 

The present authors claim that reaction solution tends to attack the 
solid over its entire surface especially at points where imperfections in 
the crystal lattice occur. Because of the exothermal nature of most such 
reactions, the local production of heat will favour attack in concavities. 
This should counteract to some extent lack of renewal of the solvent in 
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surface depressions !). On the other hand non-reaction solution will 
favour assault on angularities and crystallographic edges because here 
the faces of greatest solution velocity are potentially present, while in 
concavities the endothermal nature of the solution process as well as 
lack of renewal tend to decrease the rate of solution. 

As a result of this difference, reaction solution will usually frost the 
surface of a grain, for instance in a weathering process, whereas non- 
reaction solution will tend to round and smooth a particle, for instance 
a grain of quartz in water undersaturated in silica, e.g. sea water. 

In a saturated solution showing slight temperature variations, such as 
many rivers, sharp angularities will tend to dissolve while narrow cracks 
and indentures will tend to fill up because there is less renewal of the 
liquid in the latter. 

In the case of dew the solution will round off the angularities and 
re-deposition will tend to fill in cracks and hollows thus softening the 
eolian mechanical roughness. But evidently it leaves a finely granulated 
(hence frosted) surface owing to the evaporation. 

These various contentions were verified experimentally in a preliminary 
fashion, but a full confirmation will require a much more extensive 
treatment. 

That physical solution can round a crystal is easily demonstrated and 
that industrial etching is usually attained by reaction solution is well 
known. The writers have shown that a mechanically frosted crystal of 
alum is smoothed and becomes transparent again by immersion in a 
saturated solution for a few hours. 

The present authors have investigations under way on the quantitative 
significance of solution for the rounding of quartz grains, on the chemical 
and physical nature of quartz solution, and on the surface finish caused 
by various kinds of solution and deposition. 
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weak one a typical frosting. 
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STRESS RELEASE IN 
PLASTICALLY DEFORMED SILVERCHLORIDE BY ANNEALING 
BY 


R. G. DE LANGE!) anp W. G. BURGERS 


(Communicated at the meeting of March 25, 1961) 


Summary 

It is shown that the reduction of the internal stresses in the individual 
crystallites of a strained silverchloride testpiece on annealing takes place 
at a different rate. Single crystals recover at a faster rate than poly- 
crystalline testpieces. 

l. Introduction 

The glide mechanism of silverchloride crystals and the character and 
pattern of the internal stresses produced by plastic deformation was 
studied by Nye (1949) by means of the photoelastic effect. Reerystallized 
polycrystalline plates of silverchloride, deformed by simple extension, 
show birefringent bands in individual grains. By combining their direction 
with the orientation of the crystal lattice (deduced from a Laue photo- 
graph), Nye found that glide took place in ¢110> directions, the glide 
plane being not always a crystallographic plane, lying somewhere between 
{221} and {111}. This was interpreted as pointing to so-called “pencil 
glide” 2). 

Analysis of the pattern of birefringence led to the conclusion that the 
glide packets were bent about axes Lying in the glide planes and contained 
a set of dislocations having either all the same sign or at least an excess 
of them. 

Annealing of the deformed plates at a temperature of about 350—400° C 
brought about a general reduction of the birefringence, without altering 
its general distribution. As Nyx points out, this might occur by a (partial) 
straightening-out of the glide-packets, involving a removal of the dis- 
locations to their edges. While this is considered possible in deformed 
single crystals, it is thought that in a polycrystalline aggregate, where 
each grain is embedded by surrounding grains, such a process could not 
occur without setting up stresses in neighbouring grains, which would 
require definite atomic re-arrangements. 


') At present at the Metal Research Institute T.N.O.-Delft. 

*) According to a recent paper by CARNAHAN, JOHNSTON, STOKES and Li (1961) 
pencil glide oceurs only at relatively high temperatures (26° C and —72° C); at 
—196° C, however, deformation apparently takes place on fewer slip systems. 
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silverchloride by annealing. 


Rig. 1. 


PLATE I 


(100) 
{110} @ 
al} «A 


G experimental pole of glide-band 


Birefringent bands in plate-shaped silverchloride crystals after 6 % 
extension; crossed nicols. Direction of extension — (x 170). 


a: normal to plate ~ //<111>: poles of glide-bands near {111} 
b: normal to plate ~ // <100>: poles of glide-bands near {110} 


PLATE II 


a as deformed; analyzer only b as deformed; crossed nicols 


d annealed — 230° in 
30 minutes; crossed nicols 


c annealed — 220° in 
10 minutes; crossed nicols 


e annealed — 350° in 
75 minutes; crossed nicols 


Fig. 2. Reduction of birefringence in plastically deformed silverchloride by 
annealing at increasing temperature. The reduction is different for individual 
crystallites: most pronounced for the crystallyte marked by the arrow. (~ 80 xX). 


PLATE II 


a b as deformed; crossed nicols 


d = c+ 25 minutes 
at 210°; crossed nicols 


c annealed — 210° in 
10 minutes; crossed nicols 


f annealed + 330° in 
55 minutes; crossed nicols 


e annealed — 270° in 
45 minutes; crossed nicols 


Fig. 3. Same as Fig. 2, other test-piece. (~ 80 x). (The white spots at the left 
on the photographs of Fig. 3 and at the bottom of those of Fig. 4 are due to a damaging 
effect of the tweezers used in handling the test-piece). 


PLATE V 


PLATE IV 


PLATE 1 


— 180° in 8 minutes — 200° in 11 minv 


ce annealed -—> 190° in b 
10 minutes; crossed nicols 
eS .- - = 
f rgd 
as deformed —> 180° in 13 minutes > 20( | minut 
7 ) i j n 5 
Fig. 6. Difference in rate of stress-reduction for two deformed silverch] single 
crystals on annealing: crystal b is “slower” than crystal a (= crystal o. 5). 
t=" 80 x) 
os 
annea! 
minutes 
Fig. 4. Same test-piece as Fig. 3, for 


(Position in mi roscop 


as deformed —> 180° in 8 minutes — 200° in 11 minutes 


Fig. 5. Influence of crystalsize on the reduction of stress birefringence in strained 
silverchloride due to annealing: the reduction increases with grainsize (~ 80 x). 
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Nye’s remark suggests the possibility that individual crystallites in a 
strained polycrystalline testpiece “straighten out’ at a different rate, 
depending on the lattice conditions existing at their boundaries, which, 
in their turn, depend on the orientation relationship existing between 
adjacent grains. If this were so, one might expect a difference in rate of 
reduction of stress birefringence for individual crystallites on annealing 
a strained polycrystalline silverchloride testpiece. Moreover it might be 
anticipated that for comparable conditions of annealing, the stress reduc- 
tion would be faster in a single crystal than in a polycrystalline testpiece. 

The present investigation was primarily undertaken to find out to 
what extent these suppositions came true. As for some of the larger 
crystals the direction and position of the glide-bands could be determined 
with reasonable accuracy, the results obtained regarding this point are 
also given. They confirm on the main the conclusions already reached by 
NYE concerning the geometry of the glide process in silverchloride. 


Experimental 


Plates of silverchloride, measuring about 307 0,25 em, varying in 
crystalsize from polycrystalline (with about 20 crystallites) to nearly 
monocrystalline (some plates actually consisted of only one crystal) were 
prepared by recrystallization of rolled silverchloride plates. The method 
used was similar to that applied by Burcers and Tan Korn Htok 
(1946) and by Nye (1949). In the present research the original plate of 
1 mm thickness (obtained from Harshaw Company) was rolled to the 
final thickness of 0,25 mm, annealed at 270° C for 14 hour (to produce a 
fine-grained material), then strained by 4 % and finally recrystallized at 
400° C for 2 hours. 

In order to produce birefringence, the testpieces were extended by 
approximately 6 %. Distinct birefringent bands became visible in most 
crystals, either in one, two or three parallel sets. Figure la and 6 give 
two examples. Following Nys’s method, for some of the larger crystals 
both the direction of the bands and the angle of tilt of the testpiece at 
which the bands became sharpest (in this position one views the glide- 
planes edge-on) were measured under the polarization microscope. Com- 
bining these data with the lattice orientation as deducted from a Laue- 
photograph gives the crystallographic character of the corresponding 
glide-planes. As pointed out by Nyn, this method can only be applied 
if the glide-plane is not too far inclined to the normal of the crystal plate. 

In five cases, for which unambiguous results could be obtained, the 
pole of the glide-band was close (within 5°) to the pole of a {111}-plane, 
the accuracy being insufficient to state with certainty whether this 
deviation was real or not. In four of these cases, if {111} was considered 
to be the glide-plane, the value of the applied shear stress was maximal 
for the <110)-direction within this plane which enclosed an angle with 
the crystal plate close to 45°. In one crystal (fig. 1b), the normal to the 


348 


plate was close to a ¢100>-direction. The poles of the two perpendicular 
systems of birefringent glide-bands lay in this case intermediate between 
{111} and {110}, although closest to the latter. In the position of the 
crystal both these planes contain ¢110>-directions, making an angle with 
the direction of extension within 10° from 45°. 

The above results agree mainly with those found by Ny. Concerning 
the possibility of glide along {110}, it is remarked in his paper that for a 
crystal, extended parallel to <100>, the resolved shear stress would be 
maximal for glide in a {110} plane in a <110>-direction. This was found 
by Srepanow (1934) under such conditions. Our crystal 2 lies close to 
such an orientation. 


Annealing experiments 

For the annealing experiments the (Leitz—-Panphot) microscope was 
provided with a small furnace mounted on the specimen table, so that 
the testpieces could be viewed while heated. This limited their maximal 
dimensions to approximately 5 mm. In part of the experiments (figs. 2-4) 
deformation was brought about by a slight torsion and stretching with 
the aid of a pair of tweezers. In another series (figs. 5-6) the testpieces 
were extended by means of a small extension apparatus. 

Figs. 2-4 refer to testpieces with fairly large crystallites. They clearly 
show that differences exist in the rate of stress reduction between 
individual grains. In all three cases it is fastest for the crystallite indicated 
by the arrow. For these crystallites the birefringent bands have already 
practically vanished after heating up to 190-220° in about 10 minutes. 
whereas in other crystallites in the same testpieces the birefringence is 
then still distinctly visible and disappears only after heating up to 
considerably higher temperatures (270-350° C). 

Fig. 5 shows the effect of annealing testpieces of varying grainsize, 
increasing from ‘‘polycrystalline”’ in a to a single crystal with one small 
inclusion in d. The heat treatments are, although not identical, approxi- 
mately the same. 

Comparison of the four series of photographs shows that the stress 
reduction is faster in the coarse-grained testpieces c and d as compared 
with the polycrystalline pieces a and b. Whereas in c and d the bire- 
fringence has practically disappeared after annealing up to 200° in 
10 minutes, it is then still clearly visible in the polycrystalline pieces 
a and 6. 

To make the birefringence disappear in these testpieces, heating up 
to 250-300° C in 30-40 minutes was required. This result fits in with 
Nye’s observation that the stresses in the testpieces used by him were 
still visible at 350°C and only disappeared after a further period of 
heating at 415° C, 

It must be remarked, however, that even two “single” crystals may 
show a marked difference in rate of stress annealing as is clear from 
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fig. 6. This is not strange, as one would expect dependence on the direction 
of the applied stress relative to the lattice orientation. causing a difference 
in internal mode of deformation (for example as to the number of active 
glide-planes involved, the occurrence of cross-glide etc.). 


Oriented polygonization 


The process of “‘straightening-out” of glide-lamellae, considered by 
Nye, involves a displacement of dislocations to some extent akin to what 
is supposed to take place during polygonization: in this latter case 
dislocations move, by glide and climb, to the two-dimensional layers 
between the individual “polygons”. This analogy would support the 
view that also the capacity for polygonization for individual crystallites 
in a deformed polycrystalline testpiece would depend on their state of 
deformation and on their “surroundings” in the deformed testpiece. 
If so, polygonization would be “‘orientation dependent” in the sense that 
definite oriented lattice regions could polygonize ‘“‘sooner’”’ than others. 

As to the state of deformation, for aluminium it has been experimentally 
shown (Burcers, Liu and Treprma (1951)) by studying the “break-up” 
of the asterism of the elongated spots on Laue-photographs of deformed 
single crystals, that the degree of polygonization due to a definite 
annealing treatment depends on the direction of the applied stress: it is 
much faster for extension parallel to a (100)-direction than parallel to 
a <(110)-direction. 

This conclusion is of interest in connection with current ideas on the 
process of recrystallization in cold-worked metals. The view has been 
advanced ((Beck (1949; CaHNn (1950); Buramrs and Treprema 1950)) 
that the “nucleus” of a new grain is a polygonized lattice region capable 
to grow at the expense of its neighbours. It has been suggested ((vide 
for example BuRGERS (1961)) that in definite metals the above selection 
might lead to the preferred growth of lattice regions with a definite 
orientation with regard to the surrounding matrix and so to the 
formation of a recrystallization texture 1). 

On the basis of this point of view it would be of interest to extend the 
purely qualitative results of this investigation to a more quantitative 
study of the stress-annealing effect in deformed silverchloride, in the 
hope to establish a definite relation, if any, between the orientation of a 


1) Tt is well known that apart from this conception other possible causes have 
(een advanced for this phenomenon; in particular selection by oriented growth 
bBeEcx (1953) (1954) (1961)), and a martensitic type of nucleus formation ((BURGERS 
and VERBRAAK (1957), VERBRAAK (1960), BuRGERS (1961)). It seems unlikely that 
poe and the same nucleation mechanism holds for all metals. Various factors, in 
narticular the presence or absence of stacking faults and the energetic possibility 
of definite dislocation-reactions might play an important part as regards the precise 
mechanism of such processes. 
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crystallite, with regard to both the direction of the applied stress and the 
orientations of the surrounding crystallites, and its rate of annealing. 


Delft, Laboratorium voor Fysische Scheikunde 
der Technische Hogeschool 
March 1961 
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MECHANICS 


ON THE DIFFRACTION OF SOUND WAVES IN A HEAT- 
CONDUCTING VISCOUS MEDIUM 


BY 


J. B. ALBLAS 


Technological University of Hindhoven 


(Communicated by Prof. C. J. BouwKamp at the meeting of January 28, 1961) 


l. Introduction 


The solution of the problem of the sound diffraction by a screen, first 
given by SOMMERFELD [1], is based on a number of simplifications. The 
screen has been taken to be infinitely thin and the medium is non-viscous 
and non-heat-conductive. As a consequence of these simplifications, the 
velocity at the edge of the screen has a singularity. It is now well known 
that the problem formulated by Sommerfeld has not a unique solution. 
Bouwkamp [2] and Merxner [3] have pointed out that, to guarantee 
uniqueness of the solution, we have to prescribe the singularities in the 
form of “edge conditions”’, which are given by the physics of the problem. 
In Sommerfeld’s solution the natural condition, that the field energy 
remains finite, is satisfied and the form of the singularities can be 
determined. 

In a previous paper [4] we have investigated the (generalized) Sommer- 
feld problem for sound waves in a viscous medium, without consideration 
of the heat-conduction. One of the results of that paper is that, even for 
an infinitely thin screen a solution exists with finite velocities every- 
where. This solution is not unique, but it can be made unique by prescribing 
suitable edge conditions. 

In a review of the paper [4], Liguruity [5] pointed out that the 
inclusion of finite heat-conductivity is imperative, if one wants to get 
quantitative results, as for gases the influences of viscosity and heat- 
conductivity on the propagation of sound waves are of the same order of 
magnitude. Now it is known [6] that the two phenomena can easily be 
combined for the problem of the undisturbed sound propagation. But, 
as was first pointed out to us by Carrier [7], the combination of finite 
viscosity and finite heat-conductivity in the diffraction problem leads to 
a rather complicated mathematical problem. 

It is the purpose of this paper to complete the previous one [4] by 
deriving a method for the evaluation of the unique exact solution of this 
problem. We restrict ourselves to give a simple approximate solution, 
which remains finite everywhere and which can be considered to give a 
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close description of the velocity and temperature distributions. This 
solution may be of some interest for the study of higher approximations 
in acoustics, as the methods of LignTrHrin [8] and Kuo [9], based on a 
first approximation with singularities, sometimes give rise to great dif- 
ficulties, as has been pointed out by Tsren [10]. 

The mathematical treatment given in this paper is based on the dual- 
integral method, given by Nose in his book on the Wiener-Hopf tech- 
nique [11]. 


2. Mathematical formulation of the problem 


We consider the diffraction of a sound wave, which propagates in a gas 
of finite viscosity and heat-conductivity, by an infinitely thin screen, 
occupying the half-plane y=0, x<0. The screen has an infinite heat- 
conductivity and the gas obeys the ideal-gas law. We restrict our con- 
siderations to the case of two dimensions and normal incidence. 

The equations which describe the sound propagation are 


0 ! 
cael 7 = {) 
oc div (ov) ; 
dv Ov Ov Ov 
ae ak ee ae pon 
ary) o( AB Bye 5) 
1 
(1) = — grad p + wAv 4 (5 H') grad div vy, 
OT 


eco = + oc,(v-grad T) + pdivv = AAT + x, 
ime Tol, 
where the following notations are used: 


t is the time, x, y, 2 are the Cartesian coordinates, 9 is the density, 
v the velocity, p the pressure and 7’ the temperature. The coefficients of 
shear viscosity and bulk viscosity are 4 and ww’, the heat-conductivity is 
4, the specific heat at constant volume is ¢y, R is the gas constant and 
is the dissipation heat. 

The system (1) is a non-linear system of equations, which is hard to 
solve. The physical conditions, however, often allow a linearization of 
this system, which can be performed by expanding the physical variables 
into series in powers of a small parameter. 

Such a parameter is 


where w is the frequency of the incident steady-state wave (of. eq. 11), 
6 an amplitude and ¢ the velocity of sound. We assume very small 
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disturbances, which assumption is expressed by a <1. We then have the 
expansions 


CP C0 ais One X02 S= aoon 
P= Po + apr + op. + ..., | 
AYE hires Neal Mp = Ee foe 
y= AV, + a2vo 
from which we derive 
Po = RooT. 
and 
c2 =. vRT o, 
with 
Y = Cp/Cv, 


where cy is the specific heat at constant pressure. 
We restrict ourselves in this paper to the first-order approximation 
and we find 


/ 0 ; 
— + 09 div v; = 0, 
Ov] a. lu ; ; 
00 = — grad pi + wAvi + (= +4 yp’ | grad div vy, 
(2) ot 33 
goto + po div v1 = AAT, 

\ A = RooT == RT 01. 
We always can pose y, in the form 
(3) vi = grad » — roty, 
with 
(4) div p=0, 
by which the equations (2) become 

a + 00 Ay = (0, 


0) ) 
005; grad p — eo FY; rot ) = — grad 
(5) 


+ wgrad Ag + wrotrotrotp + & af W) grad Aq, 
oT; 
UCase = po Ay SA ae 


where we have made use of 


(6) Av, = grad div v; — rot rot v1. 


OP 4 ; 
(7) A Ga + ~i — (50 + Mt ) 4y) = 0, 
0 
(8) rot rot (05 + wrot rot ») == 1G). 
or, using (4) and (6) 
(9) A (0 4+ wrot rot ») = 0. 


The solutions of (7) and (9) are 


4 
eos! = — p+ (50 + “) Aq + f, 
oe = —pwrotrot + g, 
with 
Af = 0, Ag = 0. 
We can take 
j=0, g=0, 


as it is possible to incorporate f and g into » and y, respectively. 
The complete system of equations becomes 


a + eAg = 0 
eos! =— 1 4 (50 «| Ag, 
(10) eos = — wrotrot y, 
goto + py Ag = AAT), 


Pi = 00 RT; 4. Roi To. 
In a two-dimensional problem we only have one component of the vector , 


pb = p-e:, 


where e; is the unit vector along the z-axis. 
We restrict the considerations to the steady-state wave solution with 
the time dependence 


(11) eit, 
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The equations for the amplitudes of 7 and y, which will be indicated by 
the same notation, become 


1M00 Cy 00 C2 
A Ae 
1 7 1 yh Ag, 
to(4/3u + mw’) y 2 . Ry 
o0c* Cc C2 
from which we derive 
vies yoo? |e —— wwRy = 
: ii io]au+myy ? ee to(4/su+bm')y tt 
Qoc® ) Q0c? 
100 Cy wy == I 2 
AAR = pram rcminerees \ eee pee eee Oe Se 
cage | ( io area | : A(1 =a ea ?. 
00c? 00 Cc? 
We introduce the abbreviations 
A= k 
A => ; 7 io(4/su+p') y > x= 00 Cp 5) a a/c 
00 Cc? 


and the equations become 


(12) Ap + kp Ayp = — es AT, 
: 1@ = ww? 
(13) Ae ey Ata Ag. 


Equation (12) is important for the formulation of the boundary condition 

for T; at y=0, x<0. We eliminate 7, from (12) and (13) and find 
toy tw ior Ay ; 

(14) AAg + Ag boy A ae tL A) sige o= 0. 


mw 


Equation (14) can be written as 
(A +k?)(4+m?)o = 0, 


where the coefficients k2 and m? are given by 
: SA 
jem =O" A gtyA +—(1-A), Bm? = eee 
% 19 HC 
For small viscosity and heat conductivity (of the same order of magnitude) 
we have the expansions 
ipAxuky?  im(1—A)?* a 


k2 = Ako? + Ae Day A = cong 


m2 = me Sr rest beet 4 ss. 
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Every solution of (15) can be written in the form 


(16) P= Pi + P2, 
and the problem is reduced to 

(17) Amn +ko. = 0, 
(18) Age + mg = 0, 
(19) Ay +2y =0, 


with [?=imgo/u. 
The parameters k, 1, m are all complex. We define the square roots of 
k*, 2, mm? by 
Imnk>0, Iml>0, Imm > 0. 


The auxiliary conditions are 
(a) the functions gy, g2 and y have to be continuously differentiable in 


the region under consideration; 
(b) the tangential velocity w vanishes along the screen: 


op oy 
any oY y 
(c) the same holds for the normal velocity 1 
op oy | 
eee a ee a =), <0; 
v dy 7 Sa ¥y v 


(d) at the screen the temperature 7; is zero. According to (12), this 
condition is formulated as 


Ap + kpyAgp = 0, 
which after introduction of (16), (17), (18) reduces to 
agi + go = 0, y=0, <0, 

where the parameter a is defined by 

7. ko? yA — k2 

kg? yA — m2? 

(e) the incident wave is 

pi = eo thy; 


(f) near the edge we assume the temperature to be bounded. 


3. The formal solution 


The solution of the problem is found by means of the method of dual 
integral equations. We introduce the abbreviations V1, Y2, ys, by 


(20) “= Ve, 
(21) v2 = VA2 — m2, 


(22) vg = VA2—22 , 
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and write the unknown functions in the following form 


. 1 @+ic 
(23) Qi = e thy + ee f PA emt extwd), ay =O; 
ZIG —o0o+1¢ ; 
5 ee l co + ie ot 
(24) Oi em eee. Vix J BA) e-#% evry dA, =, 
Z2It —co+ic 
_ il co +i¢c on 
(25) g2=—= f ClAje*te-ny J, y > 0, 
) 27 —0o+%c x 
] c+ic 
(26) m= Diijetemyds, <0, 
27 —co+ie 
oon ] cotic ; 
(27) y= ios f{  E(A) e-@* e-vsy dA, y > 0, 
ZIt —co+ic 
l oo +ic - 
(28) p= f F(A) e-# evsy dA, “y=; 


J 2a —oot+ic 
The square roots (20), (21), (22) are defined in the A-plane, cut from 
k to k+itoco, —k to —k-ioo, etc., such that 
Re (V22—k?) > 0, Re (/#2—-P)>0, Re (V22—m2) > 0, 
along the path of integration —co+ic to w+ic. 
For A4=0, we find 

Yi=—tk, ye=—Im, yz =— dU. 

The constant c will be determined later. We have 
Imc <0. 


The continuity of the functions gi, yz and y at y=0, x>0 leads to the 
equations 


oo+ic 


(29) J eda {A(l)— BA} =0, > 0, 

(30) Bike ea da (O(a) — D(a)}=0, 4 >0, 

(31) a ez dj {H(A) — F(A)} = 0, y=, 

and the continuity of se 2. st at y= 0, e= 0 to 

(32) DT exe dg A Bay 0, 2S 0) 
Seat 

(33) i ce tediyn dee D 0, 12> 0: 
—oo +ie 

(34) ic ead y3{ H(A) + F(A)} = 90, «> 0. 


—co-tic 
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The condition w=0 for y=+0, x<0 gives 


(35) cs eter dd {iA A(A) + tAC(A) — ys H(A)} = 0, «<0, 
oa 
(36) ‘ie ee dA (4A B(A) + 1AD(A) + ys F(A)} = 0, «<0. 
—o0-+ie 
The condition v=0 for y=+0, x<0 gives 
(37) th + = Ric ett dd {y, A(A) + yo C(A) + tAH(A)} = 0, 2 <0, 
(38) —tk 4 : hy ez dd {1 B(A) + yo D(A) —iAF(A)} = 0, «<0, 
27% —co+ie 


while we derive from the condition 7;=0 for y=-+0, «<0 


LA ? ‘ 
(39) a + = J ered, {aA(d4) + CiA)} = 0, 2 <0, 
256 —co+ic 
1 co+itc ; 
(40) @ + — J edd {aB(A) + D(A} =0, 2z<0. 
\ 2% —co+iec 


We proceed with the introduction of the functions 


P,(4) = A(A) — B(a), 

P2(4) = A(A) + B(a), 

P3() = C(A) — D(a), 

ae) Pa(A) = O(a) + D(A), 
P3(4) = E(A) — F(a), 

Po(A) = E(A) + F(A), 


and we add or subtract the equations (35), (36); (37), (38) and (39), (40). 

The system (29) to (40) of 12 equations is in this way separated into two 

systems of six equations each. The first system for P;, P3 and Pg is 
co+ic 


(42) J Pi(A)e-#« da = 0, x> 0, 
—co+ic 


co +ic 
(43) J Ps(A)e-“* dd = 0, z>0 


—oo+ic 
co+ic 


(44) J ee da (iAP, (A) + iAP9(A) — ysPo(A)} = 0, 2 <0, 


—oo-+ic 


Cot t 


(45) J Pe(A)yse-**dA=0, x«>0, 


—co+ic 
if co+itc 
V2a —co+ic 


co+ie 


(47) if ex dj {aP (A) a5 P3(A)} = (0); =a Oy 


—co-+ic 


(46) e-tt dA {y; Pi(2) + y2Ps(A) + tAPo(A)} = —2ik, 2 <0, 
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The second system consists of Ps, P4, Ps. We find 
oo +%¢ 
(48) f y1 P2(A) e442 dj, = (@). ie) 
—co+ic 
co+ic 
(49) f y2 P4(A) e 42d, = O, 2p (0), 
—co-+ic 
co+ic 
(50) il P5(A) e-#42 da = 0, G0) 


(51) i) ex dj {tAP2(A) + tAP4(A) — v3.P5(A)} == (i), @ =, 


(52) i ea dA {y71 Po(A) ae vo P(A) == wAP35(A)} = 0, Mp <0), 


(53) 


fe da {aP2(A) + Pa(d)} = —2a, en) 


y| 
2 
: 


V: 
We consider the first system. From (42), (43), (44), (45) we obtain 
(54) WAP (A) + 1AP3(A) — y3Pe6(A) = 0, 

and from (42), (43) and (47) we find 

(55) P3(A) + aPy(A) = 0. 


From (54) and (55) we can express P3(A) and P(A) in terms of P;(A). 
We obtain 
yea) = —aP,(A), 
(56) iA 
j Pott) = ae — a) Px(A), 


and we have for P(A) the equations 


i co +ic 
(57) oa P(A) ear di = 0, 1G 0, 
ZIt —cotic 
il co +ic 
(58) ee f w1Pi(A) K(A) e- dd = — 21k, AN, 
27t —oo+ic 
where K(A) is defined by 
(59) K(Aa) = pe iy 
VAL y1iy3 


An important property of the function A(A) is 
lim K(A) = 0, (|A| = co). 
The second system for P2, P4 and Ps; can also be reduced. From (48), 
(49), (50), (52) we find 
(60) y1 P(A) + yoPa(A) + tAP5(A) = 9, 
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from which we derive after the introduction of 
Qi = yi P(A), Qe = ye Pa(A) 
the relation 


PsA) = 5 (Qi + Qa). 


For the functions (A) and Qe(i) we obtain the equations 


oo +-4¢ 
(61) = [ OsA)e- di — 0, ro: 
2IU CO +ic 
(62) = f : Qo(A) e-4* dA = 0, e's>O), 
} “aIt —co—ic 
Manet Se A A , y : ae 
(63) on i e eal (A) ' — We(/) = (Qld) t Qs(A) | = (), 
ZIT —co+ic ae v2 A 
erin 
(Gayo genen | Se yal dug an alg lee a 
V2 -cot+ic \y71 yo | 


We first discuss the solution of the system (57), (58). In equation (57) 


we transform the x-variable into r+ & and in (58) into «—€, with 
We multiply (57) by a function Ny(£) and (58) by N2(é). 
We integrate the equations from zero to infinity. We have 
0o-+-i¢ oO 


] : ‘ ; 2 
(65) Van f fy N,(é) e-t4 dé P\(A) e- thr dA = {). eS Ap, 
Z4It —co+te 0 


] OO--ie oO oO 
(66) —— f ff No(&)e- dé Py(A) 1 K(A) eH dA = — 21k f No(E) 
\V 27 -co+ic 0 0 


E> 0. 


In the equations (65) and (66) the functions Ni(€) and N9(&) are arbitrary 


functions such that the integrals have a meaning. 
The function 
[ Mi(8) edz = N_(A) 


0 


is a regular function in a lower half-plane, while the function 


[ Na(§) e-P8dé = N(2) 
0 
is regular in an upper half-plane. We have 
1 oo +ic 


J N(a)PiAe-*daA=0, a>0, 


V 2a —cotts 


(67) 


oO +-4¢ 


(68) = J N4(A)y. K(A) Pa(d) e- da = — 26k 4(0) 


\ 2 —o-+t¢ 


: Urea) 
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Now we shall determine the functions N (4) and N_(A) so, that we can 
invert equations (67) and (68). We split 


(69) ee ee 
and 
(70) K(a) = K_(A) K,(A) 


where the functions K4(4) and K_(A) are regular and free from zeros in 
an upper, and a lower half-plane, respectively. 

As K(A) +0, if |A] + 00 the possibility of the factorization (70) is not 
trivial. We can however show, that this factorization is possible. We take 


and the equations (67), (68) become 


(71) _ to VA—k K_(A) Pi(4) e- da = 0, G0), 
| Z2It —co+tic 

(72) = ie VA=h KA) a ez dj = —2ikN ,(0), Baa ()s 
/2It —cot+ie +T 


where 2; is an arbitrary complex constant, satisfying 
Im A; > 0, |Ai| sufficiently large. 


We multiply (72) by e-=* and differentiate with respect to x. We 
obtain 


= VA—k K_(A) P(A) e-#2 dd = —21kN'4(0)-Ay 


270 —cotic 
(73) ia en 1 


From (71) and (73) we find for P;(A) 


2 Vk I ] 
(74) False 2 Kam ata 
From (56) we find 
2 Vk l | 
(75) P(A)= a x Kx(0) Wink KO) 
Mj. ve (1—a) 1 
(76) P(A) = -il? (0) (fab Pak EK) 


First we see, that we have to take for the constant ¢ a complex number 


with 
time |= Im h, 


Ime = 0, 
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if we assume that 


Imk <(Iml, Imm). 
In the second place we remark, that the solutions (74), (75) and (76) 
satisfy the equations (42), (43), (44), (46) and (47). The integral in (45) 
diverges. It has to be interpreted as 


oo + ie a) , 
: 6(A) rs_e-t dA = 0, x> 0. 


an 
v 


Oe oasaka 


This equation is satisfied. 
To solve the system (61), (62), (63) and (64), we introduce the function 


R(A) = Qi(A) + @2(A) 
and proceed with the system for R(A) and @:(A) 


1 cot+ic 


(77) = f Qi(A)e-# dl = 0, x> 0. 
/ZIt —cotic 
co+ic 
(78) = J Rd) eda = 0, x> 0. 
ZIt —oco+ic 
l co+ie Oud ] L(a) e- da 2a 1 
== “EN OA. = : ee 
(79) V2x wis : = l—a l=% 
cf ; A 
; ~ f SNA) hed: a = 0; 
V2 -cotic 2 
ni i(1—a) d+ 2Q,(a)_, 
(80) pe f  ysR(A)e-@ da = o(f— 4) e . Q1( ) ote di. en. 
2x —oo+ic | on dx = ean V1 


where the abbrevation 


_ (yi/y2)—a@ 
(81) L(a) = a pe 


has been used. 
Equations (77), (78), (79) and (80) can be solved (formally) by means 
of an iteration method. We pose 


(82) R(A) = 0, 
and determine @ (A) from 


(83) Qld) = 2\, a eee ae 


“lal x Ll) a La 
where L,(A) and L(A) are the factors in the decomposition 
(84) L(A) = D(a) D(A). 


A second approximation is found by introducing (83) into (80), solving 
equations (78) and (80) for R(A) and proceeding with this function in (77) 


b] 
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(79). In this way, we obtain some series in powers of a small parameter 
for the functions R(A) and Q,(A). In this paper we restrict to the first 
approximation. We find from (82), (83) 


Lu 2 Vk i 1 
(a OO Sa IE L4(0) aVa +k L(A) 
86 Pill) = ar 2 Vk Va—k 1 
ee) (4) l—al x Li(0) ay72—m2 L(A) 
(87) Ps(A) = 0 


4. The factorization of the kernels 
For the decomposition factors of the kernels K(A) and L(A), we use the 
well-known formulae 


ee Le a Logan) 

log K_(A) = — —= ———— dé, 

og K (A) Vani oi rary | S 
Le a oe KE) 


V2mi -o+iqg §—A 


dé, 


; 1 ~+# log L(é) 
log L_(A) in =, J = de 
/ ZIT —co+ip 


il er Now Le) 


V2ni —co+iq =A 


log L(A) = + dé, 


where the numbers p and q lie in the strip of regularity and 4 satisfies 
Qaim A= p: 


If Im /A<0, we have 


; A log K (Ee) 
= Sey 
(88) Ka) = exp| — =; J nae | 
while for Im A>0, we find 
: A © log K(&) 


From (88) and (89) we derive 
(90) K (4) = K(—A), 
and with the help of this equation we find the value of A4(0): 


(91) K (0) = K(0) = 


From the asymptotic behaviour of K(A), 
tt al 


— 5 ell? + (1a) Bam), ([A] +0), 


(92) Ki) ~ 
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we derive 


(93) K_(A) eae 3, dee 
(94) iipee a pe ee iia}, Secet: 


The signs in (93), (94) are determined by contour integration (cf. [4]). 
From the equations for log L(A) and log L_(A) we obtain 


1 k—am 
(95) L.(0) a L_(0) = on pees 


and the asymptotic formulae of the kernels L.(4), L (A) are found from 


 me—k2 | Oy net 
(96) Lp) col + ical i (|A| > co), 
to be 
(97) E(A)~1, (|A| +00) 
(98) L4(a) ~~] |A| —> co) 


5. The velocity and temperature fields 


From the results we have obtained, we derive expressions for the 
velocity and temperature fields. We find from (23) to (28), (41), (74) to 
(76), (85) to (87), (91) and (95) for y>0 


. a k ih l l = 
u = — —— —— — ez. 
27 V/k—am -co+ie VA—k K-(A) 


(99) ¢ - {e-7'¥ —ae-”¥— (1 —a) e-” dZ 
l «@ Vkm etic |] | 
-— ——_- - iAx envy =— my us| Jae 
| 2% V1 —a Vk—am -co+ie VA 4—(A) ye : 
OO+t 7] 
D ik e iky = h f | A k l 6 tag 
“% Vk—am-co+ic 4A K_(A) 
d | ¥1 v2 A2 , | a 
(100) env —a —e-*y —(] —q) e-vy\ dj 
1 yiy3 | 
L View VIRB e 
+ > —— —_—_ —— et Je-v1y —_ a-vey\ dq} 
| 22 V1 —a Vk—am -coo+ie A L(A) e 7 tit 
toRy A 
| @ keypA ame Tt = OP + 90 = 
- al k ort 1 l 
101 We ae tky euies i = e-tAx Jetty =e —vay} 
sa) 20 Vk —am —oo+ie AVA—k K—(A) : 
= is bik yicg 7 ae : eta yilye yee ele Je! we "ay 
| Ve — am. corie AVA +k L(A) l—a ie 
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For y<0 we have 


/ 


v k cotic 1 
a LL  _—_—_ J === e-tat . 
a7 Vk — am —co+ie VASE kK _(A) 
(102) « . {en —aen¥—(1—a) env dd 
ay a Vkm co +ie l 1 
= [SS | e= y— 4 An J 1y ya ; | 
2x Vi—aVk—am —o+ie VAL k L(A) © ie ra ery | d 


v= —tke-tky . — finer te e-taz 
2% Vk—am-corie & K-(A) 
if @) 3 : . { bel v2 VARY =e 0 | 
(103) e a—e 7s | la) evsy! dj 
| Vi yiys | 
1m Vim VA 1 ide Seviy — id 
SE ——— — e 40 Jey ery 
2% V1—aVk—am-—aoric 4 L(A) . ae 
| agi + 2 = ae thy + nal sey it ; meee a ether . 
21 Vk—am —cotie AVA—k K-(A) 
(104) ¢ 0 {evry — ey} d, 
uw yk Ves a 1 1 re ne Vee wer u 
— ——— | ] — ea —+—_______! dj. 
2m Vk —am —cotic AVA+ hb L-(A) | l—a | 


From the equations (99) to (104) we derive the integral representations, 
which hold at y=0, 


. l a Vkm tie J VA2—k2\ e-tte 
(105) w= — ——— == (1- ara 1A, 
7 \/ ‘1—a Vk—am - Lie VA+k V 22 — m2 L(A) ; 


k 00 -+i¢ | ye k 
[ 


106 v= —itk+ — K (a) e- 2 da, 

: 2a Vk—am —o+ic A ook = 

(107) BO Pe Ee = a =: (ee wae et dA. 
4 Vk —am —~cotic AVA+K 


At the edge x=0, the velocities and temperature remain finite. This 
follows from the asymptotic expansions (93), (94), (97), (98). More precisely 
stated, the velocity wu tends to zero like x3, while v and 7; tend to zero 
like Vx. 

For the case of small viscosity and small heat-conductivity the con- 
tribution to the diffracted field of the exponential terms with yey and 
y3y can be neglected at some distance from the screen. As can be expected, 
the influence of viscosity and heat-conductivity is restricted to the 
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boundary-layer regions in first approximation. The kernel function A_(A) 
also contains terms with / and m and the form of the far-field will in second 
approximation depend on these parameters. A complete discussion of the 
fields will not be given here. We refer to [4] for the case of viscosity only. 

For the study of higher approximations, we need the numerical evalua- 
tion of the integrals. This can be simplified by taking |1| — oo, |m| — co 
in (99) to (104) for |y| > 6> 0, and using the integrals (105) to (107) for y=0. 


6. Verification of the formal solution and uniqueness 

The analysis we have given is a formal one and we have to prove that 
our solution satisfies the equations (2) with the boundary conditions 
(a) to (f). 

The functions wu, v and 7; are derived from qi, gy2 and y, which conform 
to (17), (18) and (19), as differentiation under the integral sign is permitted, 
because of the exponential character of the integrands. The functions w, 
vand 7; are continuous at y= 0, as follows from the integral representations. 
The derivatives of these functions are also continuous at y—0, for x>0. 
The functions dv/dy and 07;/dy are not continuous at y=0, «<0. 

From (106) and (107) it can easily be seen that the boundary conditions 
v=0, for y=0, «<0 and 7, =0 for y=0, «<0 are met with. The boundary 
condition w=0 is not exactly satisfied, as we have used an approximate 
solution of the system (77) to (80). The value of w can be expected to be 
very small everywhere and we can improve the result in the way in- 
dicated. However there is no need of such an improvement. 

The question of the uniqueness is important. We can add to our solutions 
other ones, obtained by differentiation with expect to x of the formulae (99) 
to (104). The new solutions meet all the conditions, with the exception 
of the edge conditions: the velocity and temperature fields become 
infinite as x0, y=0. The energy integral however remains finite. 
The stringent condition that in this problem the temperature must be 
bounded everywhere, leads to the exclusion of the fields with singularities, 
as the temperature and velocity fields are coupled. 
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Series B 


PHYSICS 


THEORY OF TRANSPORT AND RELAXATION PHENOMENA 
IN A GAS MIXTURE. I 


VELOCITY DISTRIBUTION AND TRANSPORT PROPERTIES 
BY 


M. J. OFFERHAUS 


(Communicated by Prof. J. pe Borer at the meeting of January 28, 1961) 


9. Introduction and Survey *) 

In a previous article [1]?) we have given a theory of relaxation 
phenomena in a single monatomic gas. In the present article, we will 
give an extension of this theory to the case of an N-component monatomic 
gas mixture. Like in the one-component theory, it will be shown that 
the kinetic stage, of a duration of a few ‘‘times-between-collisions’’, is 
characterized by relaxation phenomena leading to the hydrodynamical 
stage. In describing both stages, stress will be laid on those elements of 
mixture theory which differ considerably from the treatment of the 
single gas. 

In part I we give a short discussion of the velocity distribution and the 
transport properties in the hydrodynamical (HD) stage (ch. 10), summa- 
rizing the treatments first given by CHAPMAN (see [2]), ENSKOG [3] and 
HettunD and Urniine [4] for a binary mixture and the one given by 
HIRSCHFELDER e.a. [5] for a mixture of more than two components. This 
survey is necessary for the derivation of the basic equations for the 
velocity distribution in the kinetic stage (ch. 11), as well as for the theory 
presented in chs. 12-17 of this publication. 

In part II we discuss the entropy production in a mixture during the 
kinetic and HD stages (ch. 12). This discussion is followed by some 
remarks on the use of a variational principle for the approximative 
determination of the transport coefficients (ch. 13). 

Finally, part III gives an account of the properties of a multi- 
component Maxwellian gas mixture (characterized by a r-4-repulsive 
potential between all pairs of molecules). This account begins with the 

") Some abbreviations used thoughout parts I-III: 

BE Boltzmann equation; HDE hydrodynamical equations; HD hydrodynamical; 
MB Maxwell-Boltzmann; ChE Chapman—Enskog ; LHS left hand side; 
RHS right hand side; w& of the order-of-magnitude of, or: approximately equal to. 


3) We will occasionally refer to the chapters of [1], numbered 1-8, and to the 
equations of [1], numbered (1.1)-(8.17) 


369 


eigenvalue theory of the linearized collision operator 2;; (ch. 14). An 
outline of the transport properties of the Maxwellian mixture is given; 
the diffusion in a binary mixture is fully described (ch. 15). After this 
treatment of the HD stage, the kinetic stage is discussed (ch. 16). The 
relaxation of the velocity distribution and the quasi-relaxation of the 
transport currents are characterized by a time scale, in which the 
“diffusional relaxation time” tp is of special importance; this depends 
on the relative concentration of the mixture components (chs. 16, 17) 


10. Velocity distribution and transport currents in the hydrodynamical stage 


In the following, the subscript i (or j, h, k) will denote a property p; 
of the i-th component of a mixture; all equations for p; will hold for 
a=1,..., N and all summations di will extend over i=1,..., N, unless 
otherwise stated. 

We start from NV Boltzmann equations (BE), one for each component /: 


(10.1) defr=—e- HF Bt SI fein) =dsfet refs 


indicating that the velocity distribution f;(¢; r, t) changes due to streaming 
and to acceleration by external forces, F;, which may vary with i and 
depend on c: 


(10.2) Valea ee Rtas ak reo 


as well as due to collisions: 


( dcfs= YI (fifn) = J de J dy sin x f der oeajlg, x) {fale’) fer’) — 


( — file) fi(er)f- 


c’ and ¢,’ denote velocities-after-collision of particles 7 and 7 with initial 
velocities ¢, ¢ (|c—¢;|=g). All interactions Vj; are supposed to have 


(10.3) 


spherical symmetry: Vij= Vij(r). 

It is sometimes convenient to consider a function p; of ¢, which has a 
different character depending on the value of i—e.g. a function which 
is connected to the molecular mass m;—as a function p;(c) of the con- 
tinuous variable ¢ and the discrete variable i(=1,..., N). In this sense, 
one can state that there are N + 4 independent collision-invariant functions 
741 of « and ¢, with the property 


(10.4) gile) + x5(e1) = yale’) +(e’), 
VIZ. 
(10.5 a, b, c) if) = Ome; Y= mic; 08) = din, C?. 


(10.5a) implies that a particle of species 7, when colliding, retains the 
amount of mass of species / it carries —7.¢., mi if k=i and zero otherwise. 
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The integral properties analogous to (2:9,51) include a summation 


over 7 and 7: 


(10.6) 2 2, fide yl%\(c) I(fifn) = & J de x cfi=0 
(10.7) SE Ide x Me) (J (fign) +4 (gifi)} =. 


The average value of a function G;(¢) for the molecules of component ¢ 
and that of a function G@;(c) for the mixture as a whole are defined as 
follows: 

(10.8a, b) éG Eee , GH BO va 
defr—=n yini=n : 

(n; is the number density of component 7; 2;=7;/n is its relative con- 

centration). In particular, we consider the moments of the collisional 


invariants +!); 
lef 0a) lef 
(10.9) yi a Ni (0A) st pal ue nO, 
The N+4 values of p4@(A=1.1,..., 1.N, 2,3) are of course connected 


to the N+4 basic moments: 
(10.10) Or=mim:; Co=me/m; T=(3k)1mC2, 


and inversely; here, the thermal or peculiar velocity is defined with 
reference to the ‘‘mass-averaged”’ velocity ¢o: 


(10.11) C=c—c, so mC=0. 


Equations for the f-derivatives of the y™ are found by multiplying the 
BE (10.1) by yi and taking >; de: 


(10.12) de pAl= —P p+ 4 > Fi, 


i 
in which the “step-up” and “‘step-down” moments ~™=+ are defined in 


analogy to eqs. (2.10a, b). Remodelling these equations, we get the usual 
hydrodynamical equations (HDE) for a mixture: 


(a) / d9¢@e= —V-(01¢0)—V-(0%<C):) 
(aa) | 0% @ = —V-(0¢0) (0 = > oi) 
cos) O) Pee Pee 9P-Piot Sok 
(ec) fdeT = —e9- VT —(2/3kn) (P: Voot+V-q— > or Fe: <C>i)+ 
+(T'/n)V - 2 <C>. 


E 1 » " ine c : y 
The RHS contain 04, ¢o, 7 and F;, as well as the transport currents: 
<C>, and 


(10. 14a, b) P= > q=nimC2C, 
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In the case of a spatially uniform mixture, the basic idea of the ChE 
method— that of a distribution depending on r and ¢ through the N +4 
moments y™—of course leads to a t-independent MB distribution fic), 
and we shall not repeat the argument given in ch. 2. 

The description of the HD stage in a non-uniform gas mixture is rather 
different from the single-component theory. Some important features of 
this mixture theory shall be given now, because we need them for reference 
in the chapters on entropy production and on the properties of the 
Maxwellian mixture. The following treatment is similar to that given in 
[2], [5] and [6], but our notation shows some differences. 

Once more, the ChE solution of the BE is based on the concepts of 


1) a distribution depending on r and ¢ through the N+4 moments 
yl and their r-derivatives of various order: 
(10.15) (A510) = Jaclyn 5)) 

2) series development of /;, with the local MB distribution /;(e) 
as the first term: 


(10.16) i= S fi) =f O(L+ ¥ bi) =f(1 +44); 
r=0 iil 


3) normalization of these terms, such that the first of them leads 
to the correct moments py”! and all higher ones have zero moments: 


(10.17) yoy; yawM—_..—0 (pir) = >. f de f,™ {1 ; 
4 


4) separation of groups of terms in the BE with the use of a para- 
meter 0. 

Application of these principles leads to a series of equations closely 
resembling (2.31), the difference being that 

a) J(ff,) is replaced by Di J(fsM fn); 

b) all other f (or ¢) receive the subscript 7. 
Now we introduce, instead of 2, operators Y and Z, such that 


(10.18a) Yidi Sf de f dz sin x f der fj (er) aaj(g, x) {hile) — dile’)}; 
(10.18b) Zio; S f de f dz sin x J dex fj) (e1) waj(9, x) {hi(€1) — F(e1')}. 


Clearly, their sum: 


(10.19a) 2b SF Vidit Zui 
and, in particular, the sum-over-): 
(10.19b) 26 SY Qud 

v) 


are found in the following approximation of the collision term for 
small ¢;: 
(10.20) DI (fifa) ~ — fi 209, 

j 


372 


as well as in the (exact) equations for successive terms of the >ifi” 
series, which we can write in the form: 


(10.21a) 24 i) f5, ) =0; 


(10.21b) | FO =2:4 | (r>1). 


Here, the LHS F;“" depends on all f;@ with 7=1,...,N; o=0,...,r—I, 


as follows: 
r—1 rl ; 
(10.22) Fy, (1/fi) {ds fi) — > 2-1-2) f,(e) 4 > > J (fi fir! @)) >. 


e=0 j e=1 


> 


The “partial t-derivatives” are of course defined by: 


def st Ouj 


(10.23) d¢) ug(elp) = D4 (8) yl, 


in which 0;) py! is the extra contribution to the RHS of the HDE (10.13) 
caused by calculating the currents P,q,<C>; with the use of s a 
instead of Sho r=0 


r=0 
Apart from 2;; and 2;, we introduce the operators 


(10.24a) Wish S Zi bi + du > Vind 
; 
and 
(10.24b) Wid SY Wi h=2i4: 
j 


apparently, in the basic integral equations (10.21b), Q; can also be written 
as W;. The use of W lies in the property: 


(10.25a) Wij(ad)=aj;Wis hs Wi(ad)= > ajWis ¢d, 
i 
as compared to the “‘linearity behaviour” of Y4;, Zi; and Qj, which reads: 
(10.25b)  Vij(aidi) =ai¥ i his Zij(aj hy) )=ajZij $j; Qi(ad)=ai¥ bit a; Zij $j. 
Next we define scalar products: 


(10.26a,b) (Fi, @i)i Ff de fOFiGy; (F,G)& Y (Fi.Gy, 
i 
and matrix elements: 


[F,G]S (F,QG) => (Fi,2:4); => > (Fi, Qy@): 


4 a j 


(10.27) ¢ = (f, W@)= d( Fo WiG) => > (Fe Wu 
ae] 


a 


= 22 YisGi)it (Fi, Zig Gi}. 


So 
The scalar product is, of course, symmetrical in F and G: 
(10.28) (fi, Gi)i= (Gi, Fi); (FG) =(G, F); 


the symmetry relation for [F,G@] is derived from those holding for the 
Y and Z matrix elements: 


TBS, ViiG; = G;, ive A 72 
(10.29a, b, c) (Li, Vag Ga)i = ( ag L'a)e 


(1,2 Gj)i=(Gy,2 Fai; (F,G]=[G, £]. 
The diagonal matrix element of 2 is non-negative: 
(10.30) [eee 0, 


=holding only if F is one of the N+4 collisional invariants 7, or 
if F=0. 

For the system of N coupled integral equations (10.21b), at a fixed 
value of the index r in the Hilbert series ¥,/;, there are N + 4 integrability 
conditions: 


(10.31) (FO, 7%) — 0, 


To the solutions ¢;(i=1,...,. NV) one may add a linear combination 
of the N-+4 solutions to the corresponding homogeneous equation 
2,6 =0, v.e., a combination of the N+4 collisional invariants 7;!. 
We have already used up this freedom by applying the normalizing 
conditions (10.17): 


(10.32) (2, 704) =0. 


Again, in proving that the conditions (10.31), at a certain 7, are fulfilled, 
one makes use of the normalization (10.32) of all ¢ up to r—1. 
We quote the explicit form of the first equation (10.21b), r=1; ef. eq. 
(10.22): 
° (=~) 


(10.338) FQ) = —2C,* C,*: Veo + (§—Ci*2) C-V In T — C- (di/az) = 244%, 


in which 

(10.34) dj 24 Va; —aif(mi/m) —1} 7 In (nkT) — (0i/nkT) {F;— (mF /m)}, 
so that 

(10.35) ¥ d;=0. 


Thus, the expression (10.33) for F;@) contains N gradients dj, ..., dy, 
which are related by (10.35). As a consequence, the solution to (10.33), 


which reads: 


(10.36) $i =AO-7 In T+ AM: Foot > Adj, 
) 


contains N vector coefficients A;%(j=1,..., N) which are dependent (one 
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could always add a constant vector to all of them). In order to fully 
determine these coefficients A, we arbitrarily require that 


(10.37) > x Aid =0 
j 


for each 7. 
For later use, we mention that the coefficients 4A must have the following 


vector or tensor shape: 


Aji = Ay(O,*2) Cs Ag) = Ay (Oy*2) C,* C;*: 


a 
(10.38) Aj 9) = A; 9(C,*2) C. Geek stcaihh) 


Now from (10.36) we eliminate one gradient, say d,, using (10.35). 
By substituting (10.36) into (10.33) and equalizing coefficients of the 
fe) 


independent gradients: Veo, V In 7’, d; (i Ah), we obtain separate integral 
equations for A: 


(10.39a) Gi = (§—O,*2) C=2, A; 
(10.39b) G,) = —2¢,*°C,*=0,AM): 
(10.40) (C/2:)(Sin — 043) = (AD — AM), 


of which we simplify the last one by taking Sjx; over both sides and 
applying (10.37): 

(10.39¢) Gi) = C{1 — (dy/a;)$=Q, AD (j # h). 

Of course, from AY with 7 +h, the coefficient A” is found by using (10.37). 


Kgs. (10.39a, ¢) are supplemented by the normalization (10.32), which 
leads to: 


(10.4la,b)  (AM,mC)=0; (AW, mC)=0. (J== 1, o.05 ZV) 
If we introduce “transport forces” as follows: 


ie} 
(10.42) XM=—-7InT; X“)=—eg; XD =—dy, 
we can write: 


(10.43a, b) PYM = —LG4OXO; gM=—F Aw Xe, 

In this notation, the integral equations (10.39) and the normalizing 
conditions (10.41) are summarized as: 

(10.44) GO=2:A); 

(10.45) (A, mC)=0. 

In ‘all this, 2=4, 4; 1, ..1, JV. 


For the most useful definition of “transport currents” in a mixture, 
we introduce the unweighed average velocity of all particles: 


(10.46) u=Cc, 
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and the peculiar velocity of a particle referred to u: 


(10.47) ae ots 


Apparently there is the following difference between u and co, and 
between C and y: 


(10.48) €—y=u—aq=c_q=C. 


In the following, we do not consider the heat current q (eq. 10.14b) but 
the enthalpy current T which, providing that the distribution fi does 
not depart too much from /;, coincides with the heat current defined 
with the use of y instead of C: 


(10.49) Ss q—3nku NN tiny? y. 


The three “transport currents” are now given by definitions which ensure 
that each product J“ X® has the dimension of an energy density: 


(10.50) JO=T; JM=P; JO=nkTLy),. 

It is easily seen that, in the approximation f; » f;(1+ 4), 
(10.51a) J®=—kT(G®, 4%); JM = —kT(GM,6); JO= —kT(GO, 4), 
or, in short, 

(10.51b) J® = —kT(G, 4) Gris lee) 


Substituting (10.43b), we get the phenomenological equations connecting 
the currents J to the forces X by means of coefficients L: 


(10.52) JO=bT(GO, YAM XO) = LOOX, 


in which, using (10.44), 


Le) =hTo,(G, A) =k x,,[A@, Aw] 


(A= 041 =. — ON oy =8), 


(10.53) 


because, on taking a vector X out of the scalar product (e.g., X” out of 
(G?), A®-.X)), one has to put in a factor 4, and likewise one includes 


1 for a tensor X. . 
On the other hand, the transport coefficients are usually introduced 


by writing down the phenomenological equations in the form: 
(10.54a) Ps —2u Veo; 
(10.54b) T= —’VT—nkT ys Dd; 


(10.54c) (i= —DrjyV In T— > Dydi. 
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Comparing (10.53) and (10.54), we can express the L“*? in terms of the 


usual transport coefficients: 


intl A | 328 N 
ite 2M 0 (0) ate 0) 
A 0 VT WL Dir zee MT Dong 
(10.55) Lee) = ] 0 nkT Dy nkT Dy, Ty nkT Dyn 
N \ 0 nkT Dorn nkT D4 nie nkT Dyn / 


Inversely, using (10.53), one can put down expressions for the transport 
coefficients linear or quadratic in the A: 


viscosity coefficient: 

(10.56a) = GkT (GH),AM) = ZT [AM, AWM] 
heat conductivity coefficient: 

(10.56b) N= tk (G®,A®) = tk [A®, AM] 
diffusion thermo-effect. coefficient: 

(10.56c) Dy, = (1/3n)(G, AM) = (1/3n)[A®, AW] 
thermal diffusion coefficient: 

(10.56d) Dyj = (1/3n)(G%, AM) = (1/3n)[AM, AM] 
multi-component diffusion coefficient: 

(10.56e) Dj = (1/3n)(GM, AM) = (1/3n)[AD, AM], 


These results imply a number of relations between the transport 
coefficients, which follow from (10:29, 37): 


(10.57) Dy =Dris Diy=Dy; 
(10.58) »tUDri=0; SY xDy=> xjDy=0. 
i i j 


Kgs. (10.57) bear out the symmetry of the L matrix (Onsager relations). 
In (10.56b), 4’ is not the coefficients of heat conductivity as it is usually 
defined. To bring the heat current equation (10.54b) in a form containing 


the usual A, we introduce multi-component thermal diffusion ratios kr; 
according to: 


(10.59a) Dy= z Dykri: 
(10.59b) O= > kri, 


(here (10.59b) has been added because the NV eqs. (10.59a) are dependent) 
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and we replace d; in the RHS of (10.54b) by its value as following from 
(10.54¢); this leads to the hybrid expression: 


(10.60) (Pes —AVT +nkT » kori Gy i 
in which ; 
(10.61) =i =TWs d kei Dri 


is the heat conductivity in the usual sense. An elegant expression for / 
in terms of the coefficients A is found by introducing 


(10.62) a) SAM — YF kp, A 


Using this definition and substituting (10.56b, d) into (10.61), we find: 
(10.63) A=1k[A®, 0] 
and by substituting (10.56d, e) into (10.59a), 
(10.64) [A®, a4 | =0. 
Combining (10:62, 63, 64), we obtain 
(10.65) A=tk[a®, a] =2'’—4k[> hry A 2 kp, A®] = 1-2". 
i 


A, ’ and 4” are all >0 on account of (10.30); A’>A implies that, for a 
given V7’, the heat current in the absence of concentration gradients is 
greater than the one holding in the absence of diffusion currents. In fact, 
the former situation will develop into the latter if a fixed V7 is maintained: 
gradients Vx; will establish tending to reduce T by means of the diffusion 
thermo-effect (Dufour effect). 

Introduction of kp; in the diffusion equations (10.54c) gives: 


(10.66) CY, f= — > Dyilkri V In T +-dj;). 
For use in the expressions giving the irreversible entropy production 
(ch. 12), we now give a formulation of the phenomenological | relations in 


terms of independent forces X, currents J and coefficients L, instead of 
the present X, J and L which are subject to 


(10.67) YmJO=0; > XO=0; det L=0. 
For this purpose we define N—1 forces x@ and N—1 currents J by 
(10.68a,b) X@=—d;; JO =nkT yy, = nkT (yi —- yy) (= 1, NI) 


and N—1 thermal diffusion coefficients Dj as well as (N—1)? concen- 
tration diffusion coefficients Di by: 
(10.69a) Dri=Dri—Don / 


iNet 
(10.69b) Diy = Di—Din—Dyj+ Dyn \ 
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Instead of (10.59a, b), N—1 independent thermal diffusion ratios hyp; 
(which coincide with the old ky; for i=1,..., N—1) follow from 


(10.70) Drjy=> Dyukri; 
i 


»y indicates a summation over 1=—1,..., N—1. . 
Now both forms of the heat transport equation (10:54b, 60) remain 

formally unchanged: 

\r = —A' VT —nkT ¥ Day d; 

(10.71a, b) 2 7 

= APT +nkT S key <3; 

\ j 


the same holds for the equivalent forms of the diffusion equation 
(10:54¢, 66): 
\ <y>3= —Dry Vin T—> Dy di 


Des 2a) }  =-Sdplkrvin? +d). | 


j=1,...,N—1. 


In particular, for a binary mixture, 7 and j take the value 1 only; eqs. 
(10:70, 71, 72) now read: 

(10.73) Dm=Dykr; 

(T= —’7T—nkT Dr dy 

10.74a = a“ 
pee ag, ( =-AVT+nkT kr: <Y)1; 
(<y1=—DnvinT—Dy d, 


(10.75a, b) ( = —DathenY in FS d;). 


Of course, <y>: is the usual diffusion velocity: Dy, and kr; are the usual 
he A = 

coefficient and ratio of binary thermal diffusion; D,, differs from the 

binary diffusion coefficient D by a factor 2 22: 


(10.76a, b,c, d) <y1=<C : «Co ; Da= De® kn ~kyp ; Pe D Jay x2. 


Finally, we should mention that a formal solution of the integral 
equations (10.21b) with the use of eigenvalues and eigenfunctions of the 
collision operators shall be discussed in connection with the Maxwellian 
gas mixture (chs. 14, 15), 


Ll. Velocity distribution in the kinetic stage 

We shall now tackle the set of BE (10.1) for an N-component mixture 
from an starting point more general than that of the ChE theory, like 
we did for a single gas in chs. 3 and 4. First we shall discuss the case of a 


spatially uniform gas, without external forces; there are NV coupled BE 
reading 


(11.1) defi= > I (fifi) 
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and the N+4 moments y™ are not only uniform in space, but also 
constant in time, on account of the integral property (10.6). 

In ch. 10 we have already stated that the ChE supposition of a 
distribution depending on r and ¢ through y™ leads to the equilibrium 
solution f;(c); this MB distribution contains the given values of p4 
(or g@, €¢9 and 7’) as parameters. In the more general method of solving 
(11.1), there is room for an initial condition: 


(11.2) fi(0) = fi(e: t=0). 


We write both the solution and its initial value as f; +a correction term: 


(11.3) file; t)=fi(e; yp) + file; t)=fs(1 + bile; £)) 
(11.4) fi(0) =f, (1 +4;(0)) 


so that both ¢;, for which we are looking, and its initial value ¢;(0) have 
the normalization conditions: 


(11.5) (d, vy!) — 0); 
(11.6) (A(0), 74) =0. 


As we intend to give here a combined discussion of the uniform and 
the non-uniform case, we immediately write down the general form and 
the initial value of the general, non-uniform solution to (10.1): 


~ 


(17) file; r,t)=fc-2(clyl4) + fic; r,)=fpO{1 + bit bile; ee 
(11.8) — f,(0) = f,(e; r; t= 0) =f,-2 +f, hi (0). 


Once more, both the given function 4;(0) and the function which we are 
going to determine, ¢;, have the normalization (11:5, 6). 

For both the uniform and the non-uniform case, we now follow 
chs. 3 and 4: the functions f; (11:3, 7) are put in the BE (11.1,.10.1); the 
BE for the ChE parts of f; are subtracted; a series development of di is 
carried out and a parameter # is duly inserted for separating groups of 
terms. This leads to a number of integro-differential equations: 


(11.9a) oi Y = —2;6M; 


(11.9b) de hil = Fi) — 2:46 


Here F; has different meanings, depending on whether the uniform 
or the non-uniform mixture is considered; we need not put down the 
complete expressions, for they can be derived from (3.14) and (4.18) in 
the same way as F;(") (10.22) was derived from F'® (2.34b), t.e., by 


a) replacing J(f%A©) by dF (fi fn); 
b) adding the subscript 7 to all other f and ¢. 
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We assign the whole initial value ¢i(0) to the first term pil: 


~ ~ 


(11.10) bi (t=0)=4:(0) ; $i@(t=0)=...=0. 
Again, we normalize the NV functions biM(i=1, ..., N) by N+4 conditions 
(11.11) (Br), 0A) — 0; 


between these equations and the N +4 solubility conditions for (11.9b): 
(11.12) (Pir, yf) — 0 


there is the usual relation: (11.11) for 1, ...,7—1 is used in the proof of 
(11.12), during which, of course, use is made of the properties (10:6, 7). 

As we did in the preceding chapter, we postpone a solution of the basic 
equations which we have formulated (11.9), until the eigenvalue theory 
of the Maxwellian mixture has been given in ch. 14. 


(To be eontinued ) 


PHYSICS 


THEORY OF TRANSPORT AND RELAXATION PHENOMENA 
IN A GAS MIXTURE. II 


ENTROPY PRODUCTION AND THE VARIATIONAL PRINCIPLE 
BY 


M. J. OFFERHAUS 


(Communicated by Prof. J. pe Borer at the meeting of January 28, 1961) 


12. Entropy production 
We start from a definition of the entropy density in a gas mixture in 


terms of the distribution, analogous to eq. (5.5): 


(12.1) s(rj=k > J defi{1 +3 In (m;/h)—In fi}. 


Like in the single-component case, 2s can be divided in two ways: 
(12.2) a) 48 =d7(S +52) +58) 4+.,..), 


in which s® is based on f;); s@ contains, in its integrand, the square of 
the deviation 


(12.3) Gt=(faffe)— 1, 

and so on; 

(12.4) b) d¢8S=ds8+ 0C8, 

v.e. 8 changes due to streaming and to collisions. A comparison between 
the terms of (12.2) and those of (12.4), for the kinetic and HD stages 
leads to results completely analogous to those of ch. 5 (see eqs. 5: 29, 34, 37). 


In particular, the irreversible production of entropy is given by an 
expression similar to eq. (5.29): 


(i225) des ~ kD (bi, Qh) = kd, $1. 


We will now make one remark about the “streaming” term dss, which 
equally applies to the N=1 case. This term contains contributions of 
streaming proper and of external forces: 


(12.6) ws=—k Z (de Ce oe ; of) {3 In (mi/h)— In fi}. 


or c 
If we take f; ~ /;, and if F; is independent of c, the F; contribution 
vanishes, leaving 


(27) ss Y—V'S, 
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in which the entropy current s is logically defined as 
(12.8) s=k > J de fi{1+3 In (mifh)— In fe}. 
This gives, as the s-balance, 
(12.9) uWs+V-s + kd, o]. 


Comparing this to the HDE for the ¢-derivatives of 9;, ¢9 and 7’ (10.12), 
which read (providing we leave out the F;-term on similar grounds): 


(12.10) yy) +7 -pal+ —0, 


we see that in both cases the streaming terms are comparable (y+ is 
a “‘y-current” defined in a way analogous to (12.8)), but in the s-case 
there is a “source of entropy” to the amount k[¢, 4], whereas the HDE 
are conservation laws, which express the fact that there are no “‘y-sources’’. 

The result (12.9) holds both in the kinetic and in the HD stage (in 
the kinetic stage, 4; is the complete deviation (12.3), i.e. 6; +d; in the 
notation of ch. 11). In particular, the irreversible entropy production 
during the HD stage can once more be translated in terms of the forces X 
and currents J. We shall indicate two equivalent methods to obtain 
this result: 

a) starting from the RHS of (12.9), taking 6 ~ ¢@) and successively 
applying (10: 21b, 43a, 51b), we find: 


EdD, $M ]=h(GD, FO) = KG, — ¥ Ge Xe) = 
=—k ¥ X(6@,Ge)=T. YIMX, 


(12.11) 


This result provides a justification for the choice of forces X and currents J 
made in (10: 42, 50). (For the relation of the entropy production to X 
and J in the thermodynamics of irreversible processes, see also 
Dk Groot [7]. This reference gives dces= SIX, which is obtained by a 
slightly different definition of X or J). 

b) One can also take the LHS of (12.9) and consider this from the 
thermodynamical viewpoint, following PRIGOGINE [8]. The computation 
is much longer and we give only the essential points of the reasoning. 

According to (12.2), s ~ s (up to terms with fi? and higher); this 
implies that s can be expressed in terms of fi, ve. in terms of nj, ¢9 and 
7, but ¢9 drops out. In fact, 


(12.12) sO —=T (e+ p—> me 11), 
in which : 
(12.13a, b) e=snkT and ws=kT In {ni(h2/2ami kT) 


are the internal (kinetic) energy in unit volume and the thermodynamic 
potential of component i, respectively. The differential relation corre- 
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sponding to (12.12), if we use n; and 7 as independent variables, reads 


k i 
(12.14) ds(0) — 7 aT + 2 (3 Ce — rr) dn; 


or, with n; and e as variables. 


. 1 i 
(12.15) ds) = mde — dni, 
a 


which is the usual thermodynamical expression (Gibbs’ relation per unit 
volume). Now (12.14) or (12.15) can be used to determine the derivatives 
of s(r,t) with respect to both r and ¢; Cite. 


(12.16) d,s = oe > (gk — dE. 
i 


Next, (12.3) is put in the definition (12.8) of s and the first two terms 
are retained, of which the second contains one power ¢; in the integrand. 
This gives 


(12.17) sw s%+sO=T (e+ p)co—q—> pani Cedi}. 
Now we substitute the HDE for 0,7; and 47’ (10.13) into (12.16). This 


leads, together with (12.17), to the following ¢-derivative of s, corrected 
for the outflow of entropy: 


o“stV-s & 048 + V-(s 4 s0)) = 


~ FP: Veo—2q-VInT+ikn€-V nT — 


Gf 
1 
hk SO Vath S ecFe- KOs 


ee if 1 n) V (x 


(12.18) 


\ 


Kgs. (12.11) and (12.18) show that both sides of (12.9) yield approximately 
T1S JX. Apparently, the ChE character of ¢;, which was essential for 
the derivation of (12.11), has also been used in deriving (12.18): only 
for a ¢; of the ChE type, ds), V-s@) and higher terms may be neglected. 

In order to translate d¢cs into a quadratic form in the forces X or in 
the currents J, we make use of the alternative forms of the equations of 
heat conduction and of diffusion (10: 60, 66). We abbreviate these 
equations by introducing a heat current I’ instead of T and a diffusion 
force d;’ instead of d; as follows: 


(12.19a) Tl’ ST — ak Y keri Cys; 
(12.19b) di Sd; tk ViInT. 


Now egs. (10: 60, 66), together with (10.54a), form a set of diagonalized 


phenomenological equations: F 
[o) 


(12.20a) P= —2u Veo; 
(12.20b) ag a 
(12.20c) i= — > Didi’: 


26 Series B 


384 


The choice (12.19) has the merit that 7~! $JX does not change when 
we replace T by I’ and d; by d;’; instead of (12.18), we get 


(221) ALIx= eF Vo—-I"-V In T —nk 2 sy i: di’. 


Now we can substitute either (10.54) into (12.18) or (12.20) into (12.21): 


7 haa yy LOK Kw 
a i ie (Veo)? +4'(V In T)?+ 2nk ¥ Dri(V In T'-di)+-nk ¥ > Diy di-d; 
i tw 
(12.22b)| = 7 (eo)? +A(V nT)2+nk ¥ ¥ Didi’ -d;’. 
i 


The RHS equals k[¢,¢“] and it is therefore +0 even if Peo and 
VT'—0; apparently the D-matrix must be positive semi-definite. 

In fact, det D=0 and we must recur to a system of independent forces 
and currents with a positive definite L-matrix, like we did ad the end of 
ch. 10, in order to invert the eqs. (12.20c). So we introduce ‘Y ; and Dy 
by (10: 68b, 69b); we retain the first N—1 diffusion forces d;’ as in 
(12.19b) and also the heat current I’ as given by (12.19a), which is 
unchanged, because 


(12.23) D—ukT ¥ ker <y>:=T —nkT ¥ keri Yt 


Kqs. (12.20a, b) do not change, but (12.20e) is replaced by 


(12.242) cys —> Dyas’; 


a 


this Dj; is positive definite and we can use the inversion: 


> (B>Yy <yy. 
s] 


(12.24b) tea 


Now we substitute (12: 20a, b: 24b) into (12.21) to obtain: 


Finally, we give the forms of eqs. (12: 22,25) for a binary mixture, 
using (10.76a, ¢, d): 


dos x (Peo)? + AV nT _ nkD ) 
T eo —— (di ke 7 ine 


v1 ve 


1 oe 1 
] — aT P25 cet ie {T—nkTk7(<C>, —€C)»o)}2 + seal (<C>) — Co), 
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13. Hydrodynamical solutions by the variational principle 


We would like to give an analysis of the use of the variational principle 
for obtaining approximative solutions to the integral equations (10.44) 
for the coefficients A in the HD distribution function frOl—>,AwxX), 
By introducing the obvious notation: 


(13.1) Gi) =GO-G%; Ao = A,)— A, 
we can write the eqs. (10.44) in the form: 

(13.2) Gi%=02,A%; GwH=0,AM; G,0=0,;A0, 
or, together, 


(13.3) GiO=2,A® (w=A,p,1,..., N—1). 


Two methods are found in literature for obtaining approximative 
solutions to the eqs. (13.3): series development of A (CHAPMAN) and the 
variational method of Enskoc. Later in this chapter we will compare 
these methods and discuss the literature on this subject; it is convenient 
to postpone this discussion until we have given a survey of the variational 
treatment. 

By the variation method, coefficients A are determined which correspond 
to maximal values of the transport coefficients wu, 2’ and Dj; under certain 
auxiliary conditions. The physical implications of this maximum require- 
ment are seen by successively taking three situations, in each of which 
only one type of gradient is present: 

a) Veo only; 

b) VT only (all d;=—0, cf. (10.34); ¢.9:, Vas—0; F=0;V Inn=—V in7, 
so that the hydrostatic pressure is uniform); 

ce) Va; only (components other than 7 and N have uniform concen- 
trations). 

For each case, we write down the irreversible entropy production (12.22a): 


eo} 


( (a) dos ~ (2u/T) (Veo)?; 
fo) = tee aA, lol"); 


“~ 


( (ec) des & Dy(Vxj)?, 


as well as the “transport current”: 


(a) = —2uVeo; 

: (b) T=—-A'/VT; 
Pe) (ce) <€y;—<On = Py 
= — Dy Vay = —(Dy—Djy—Dyj + Daw) Va; (10.72). 


Clearly, the requirement that pu, 2’ and De shall be maximized, implies 
that, for given values of the “transport forces” Veo, V7’ and Vx;, one 


obtains: 
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1) the greatest possible dissipative entropy production, and 

2) the greatest possible yield of the transport process. 

Now yp, 4’ and Dj; are connected to the coefficients A by means of the 
eqs. (10.56a, b, e): 

(a) 1l0u/kT=(GM, AM)=[A™, AM]; 
\ (b) BA [k=(GO,AM)-—=[AG, AGT 
l(c) 8n Dy=(G, AD) =[AD, AM), 


(13.6) 


or, comprehensively, 
(13.7) c= (G), AM) = (Ae AT 


This transport coefficient 7, then, we are going to maximize. If we 
denote the approximative value of A by A and leave out the x index, 
we have the requirements that A shall be an approximate solution to 
(13.3)? 


(13.8) G ~ OA, 


and that the approximative value of + obtained by taking, in (13.7), 
A instead of A, shall be an extreme: 


(13.9) 0(G,A)=0 or 6é[A, A]=0. 


Now, instead of (13.8), we introduce the following assumptions, the 
second of which serves as an auxiliary condition for the extremal value 
problem : 


a) the solution A shall have the form of a finite development in a 
series of € known functions of 7 and C: 


a €=-1 
(13.10) A;(C)= > x (2) (2) (C), 
e=0 


in which the coefficients ~, which we are going to determine by applying 
the variational principle, do not depend on 7; 

b) A shall be orthogonal to the difference between the sides of 
eq. (13.8): 
(13.11) (A, @—QA)=[A, A—A]=0. 


The last three equations: (13:9, 10,11) are placed at the base of the 
approximation method. 


As a first consequence of eq. (13.11), we see that one approximative 
value 7 of the transport coefficient t follows from both expressions (13.7), 
with A instead of A, one linear and one quadratic in A: 


(13512) t=(G, A)=[A, 4]. 
Secondly, it follows from (13.11) that 


(13.13) [A, A]=[A, A]=[A, 4]. 
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a ee ener ees ee 
This implies that the approximation 7 is smaller than the exact rT: 


(13.14) tt (|A, Al {A.A (A — 4, A—Ale0 


because, for any F, [F, F]>0 (10.30), but the equality sign in (13.14) 
would only apply if A—A were a linear combination of collisional 
invariants, but these are absent from both A and A on account of the 
usual normalizing conditions, or if A=A, and this is, of course, 
excluded. 

_it we first look for a solution (13.10) consisting of one term only: 
A; =a o,(C), the orthogonality condition (13.11) alone is sufficient 
to fix the one coefficient «©: 


(13.15) ox (0) — (og), G)/(o, Qo) = GO) /Q(00) , 
in which we have introduced the notation: 
(13.16) Ge) = (o®, G); Qe’) = (c@, Na’). 


When the trial function A is to consist of two or more terms, its 
coefficients «@ are determined by the maximum requirement (13.9) 
together with the condition (13,11). According to the Lagrange theory, 
this is equivalent to determining the following free extreme: 


(13.17) d{(@, 4)+4(A, @-QA)}=0. 


Now the coefficients «@(o9=0,...,£—1) and the Lagrange multiplier 
follow from (13: 11,17). We write these equations using the notation 
of (13.16): 


(13.18) Y «GO — FY ae) x) Qlee) = 0; 
Q @ oe 
BS) df(1+A oc (2) Ge) — A 0 (2) ov (o’) () (02) — 0, 
l i 
Q o Qo 
tC, 
0 


f1=(1+A)GO—22 ¥ a O@e=0 (e=0,...,é—]). 


(1 


(13.20) = 


If we multiply each equation (13.20) by the corresponding «@) and sum 
over o, we get (13.18), provided we take A=1. This determines the value 
of the multiplier, and now we rewrite (13.20) as ¢ equations for the 


coefficients «): 


(13.21a) Ge) = ¥ oe) (ee’) 
oF 

or 

(13.21b) (o@, G—QA)=0. 


Thus, for €>1, the extremal condition (13.9) narrows down the orthogo- 
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nality condition (13.11), to the effect that G—Q24A is now perpendicular 
not only to A, but to each separate function o®@ in the series (13.10) 
for A. 

We will next discuss the choice of the functions ai2(C) to be used 
in the development (13.10). Of course, if the exact A is required to have 
the vector or tensor form given in (10.38), we choose all o to have this 
same character. The further dependence of ¢ on C will now be discussed 
for a simple gas. It is by no means necessary to use an orthogonal or 
complete set of functions o for the development of 4A; one does not 
generally go further than 1 or 2 terms and if these are suitably chosen, 
t will approach +t fairly well. For the coefficient connected with heat 
conductivity, one usually takes a series: 


ee! : 
(13.224) AM = ¥ a4 mCSL+Y(C*2) = mC{a (5 —C*2) 4 ar 


o=0 


and for the coefficient connected to the viscosity : 


ee | O 
vac ; : : ° int Face mec 
(13.22b) At) = 2 x2) C¥C* Si2) (¢ *2) = tr 4% (9) T DPS 
o=0 


These are developments in eigenfunctions of the operator Qy for a 
Maxwellian potential, arranged in order of increasing eigenvalues wy; 
in both cases, the first term is proportional to G, the corresponding LHS 
of the integral equation (13.5), so that, of course, one term only would 
give an exact solution if we have a Maawellian intermolecular potential. 
In order to see which choice of o must be made if we treat the non- 
Maxwellian case with a development (13.10) consisting of one term, we 
write down the transport coefficient 7, which follows from (13: 12, 10, 15) 
if we use o only: 


(0) @y2 0) Me 
(13.23) = Ae ee = (o(, @) . 
(a0), Qa()) [o(), G1) ] 


We now compare the 7 values as found by successively using, as a 
one-term trial funetion o@), the function G occurring in the LHS of (13.8) 
and a different function FP: 

(GG) (FQ? (G,@)2 LF, F]-(F,@)2[G,G] 


(13.24) At?=tg—tp =a _ LEG _ } 
) [GG] [FF] [F, F](G,@] 
G constitutes the best trial function, when this expression is positive for 


all 4G. Sufficient for At to be ~0 is, that Schwarz’s inequality : 
(18.25) (F, F) (G,@) > (F,@)2 

is supplemented by 

(13.26) (f, QF) (G, 2G) 


es a ee 


(Ff, F) (G,@) 


This last condition is clearly fulfilled, for both the 4- and -equations, 
in the Maxwellian case, if only because t@ now equals the exact t and Tp 
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does not, so that (13.14) applies. For other potentials we have argued in 
ch. 7 that (@, QG)/(G, @) can be used as an approximation to the lowest 
positive eigenvalue of © corresponding to any eigenfunction of the 
required vector or tensor character, and therefore, for any F of this 
same character but different from (, (13.26) will continue to hold. 

In the case of the N-component mixture, the functions o used in (13.10) 
also depend on 7, and, instead of (13.22), one might try to use: 
(13.27a,b) Aj) = Ya) miCS0(O;*2); Ai) = ¥ xl) C;*C,* 80) (C,*2), 

Q 0 
These developments are satisfying as the C-dependence goes, but they 
are not subtle enough to render the i-dependence of A. In fact, the use 
of (13.27b), for N =2, leads to an expression for the viscosity of a binary 
mixture: 


(13.28) ee ee ee ae 


which amounts to combining the cross-sections for the three types of 
collisions in this mixture according to the relative frequency of each 
collision; this could only serve as a very first approximation. 

For the description of binary diffusion, a choice of A as in (13.27a) 
cannot even roughly serve, because here it is essential that A) and A®) 
should, for each C value, have opposite signs. That is why we improve 
(13.27) by taking a development in €N functions: 


(13.29) a; = 04: )) =648;4(C) (P==0), ese Lal. ee I), 


which gives: 


(13.30a) AC) = > Y 0%) aj I(C) 
pa) 
(13.30b) = > aj”) s47)(C) (ats) = a,), 


Using functions « which contain Kronecker 0’s is clearly equivalent to 


allowing the coefficients « in (13.10) to depend on t. 
The scalar product of o;7(C), with another function F of « and C, 


is of course defined according to (10.26): 


(oes U) (oD), F)=> by (96, Fi)i= (81, Fa): 
therefore (13.21b) can be further narrowed down to read: 


(13.32) 


for each separate 7. . in 
This final result could also have been reached by applying the conditions 


(13:9, 11) in order to determine the EN coefficients a; in (13.30b), 
instead of the & coefficients «®) in (13.10). Such a treatment is found in 
the original formulation of the variational method by Enskoe [3], 
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as well in the accounts of this method given by HELLUND and UBnHLING [4] 
and by HiRSCHFELDER e.a. [5]. However, the discussions in references [3] 
and [4] are limited to single gases and binary mixtures, whereas in the 
treatment of [5], which holds for an arbitrary number of components, 
the variational problem (13.9) for 7 is unnecessarily subjected to N subsidiary 
conditions 


(13.33) (Ay, G;—2; A),=0 


instead of the one condition (13.11); hence one multiplier 2; has to be 
introduced for each 7, and, of course, they all prove to be = 1. The present 
treatment shows that one condition (13.11) is sufficient; from the result 
(13.32) it follows that the N conditions (13.33) are automatically satisfied, 
and it is not necessary to introduce them a priori. 

On the other hand, the approximation method for A, valid for N=1 
and 2, as given by CHapmMan and Cow1ine [2] does not explicitly 
introduce the variational principle, but it leads to identical results by 
way of cutting off determinants of infinite size. Perhaps these authors 
lay too much stress on the possibility of reaching accurate results by 
using series developments of the type (13.30b), for the rate of convergence 
of such a series will markedly depend on the potential field V(r); it is 
very cumbersome to compute coefficients a; with r>1, and, when the 
number of them, &, is increased, they all have to be calculated over 
again. Moreover these higher terms are not even of great interest, because 
the 7 calculated at €=1 or 2 are optimal and, in general, they compare 
favourably to experimental data. 

We shall finally give the explicit form of the result (13.32). To this 
end, (13.30b) is substituted into (13.32). As we also did in (10.24b), we 
replace the operator 2; by W;, which has a more simple linearity property 
(see 10.25). (We would like to point out that the rather special “‘linearity 
character” of the 2; operator has not played a role whatever in the 
variational treatment). This leads to: 


(13.342) (5, Gii=> ¥ as) (86, Wis’) 
ck 

or 

(13.34b) GO =F >¥ aj) Wir, 
ey 


if we use an obvious notation analogous to (13.16). This provides us with 
EN equations which determine the coefficients a;"), Of course, G; has a 
different meaning, depending on whether we consider the f-, A- or 79- 
equations (see 10.39), and so have s;), In fact, the coefficients a; should 
bear an additional index it, A or j, which we have omitted for convenience 
(moreover, the a vary with the degree of approximation 6). 

We mention for completeness that, in the A- and j-cases, the system 
(13.34) contains N dependent equations, viz. those with r—0; one of 
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them is superfluous and must be replaced by the normalizing condition 
(10.45): 


(13.35) (mC, A)=0 or DS aj aj;mj;=0. 
j 
The transport coefficients are expressed in terms of the a by sub- 
stituting (13.30b) into (13.12): 
(13.36) T=(4, AN= YY (G41 8) => Y aM. 
Eqs. (13: 34, 35, 36) summarize the results of the variational treatment. 
The choice of an orthogonal set of functions s“, which includes G, is 


mainly useful because all but one of the coefficients G;”") = (84), Ga)s, 
occurring in both (13.34) and (13.36), are zero. 


(To be continued) 
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IN A GAS MIXTURE. III 


THE MAXWELLIAN MIXTURE 
BY 


M. J. OFFERHAUS 


(Communicated by Prof. J. pe Borr at the meeting of January 25, 1961) 


14. Higenvalue theory 


In ch. 7 it was shown that the theory of transport and relaxation in a 
simple gas is considerably simplified in the case of the Maxwellian 
intermolecular potential. The fact that the differential cross-section 
x(g, ~) does not depend on the relative velocity g=|¢—«| of the colliding 
particles, makes it possible to find the eigenvalues and eigenfunctions of 
the 2 operator. 

In a mixture, a similar simplification takes place if one supposes that 
the interaction between all pairs of molecules (i, j) is of the Maxwell form: 


(14.1) V as(r) = xu,/r4. 


As in the single-component theory, the exact results for this Maxwellian 
mixture will be closely analogous to those holding for the general potential 
field in first approximation. Therefore we shall, in chs. 14-17, consider 
the Maxwellian mixture only; the symbols will not carry the subscript mM. 

In this chapter we shall discuss the eigenvalue problem of the operator: 


\ Qi P(e) = J de J dy sin x f der fj (e1) vag, x) - 


a) - {hi(€) + Pj(€1) — file’) — hj(e1’)} 
= wij di (€). 


The only property of the Maxwell interaction (14.1) of which we shall 
make use, is the fact that a;; does not depend on g. This property is 
characteristic of a wider class of interaction potentials, and the present 
discussion applies to all these. 

As to the Maxwell case, aj; has the form: 


(14.3) ois) = (85/Mas)*a*(z), 
in which wv; is the reduced mass of the i-) pair: 


(14.4) Mig =m mj/(mi +m) 
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and «*(y) is a dimensionless function of y common to all pairs 7-7, which 
was already determined by J. C. Maxwetu himself: for its explicit form, 
we refer to [9]. 

We shall now solve the eigenvalue problem (14.2) by the method of 
generating functions, which has been applied by WaLpMANN [6] to the 
case of a single Maxwell gas. It is clear that the solutions will contain 
the parameters nj, ¢o and 7 occurring in the MB distribution: 


(14.5) [5 (€1) = 3(mj/2nakT)"* exp {—(mj/2kT) (e1 — €0)?}. 
Of these, n; will only enter as a factor in the eigenvalue @,;;. As for ¢9, we 


shall first suppose it to be =0; it is easy to see that, once the eigenvalue 
problem for this special case: 


(14.6) Qij $(e) FQ, $(€) = ey bile) 


has been solved, we can use these eigenvalues @;;, along with the 4; at 
the value ¢—co of the argument, as solutions for the case ¢)40: 


(14.7) LOFF $(¢ —€o9) = mi; di(€ — Co). 
The temperature 7 will occur in the eigenfunctions ¢;; it disappears when 
we introduce dimensionless reduced velocities as follows: 
(toa Py 6 (cc “Jc (ig ae Py “cep er fey (m7) 2k7)?. 
In this notation, the 2;; operator for cy=0 reads: 
( 24 g(c) = a~" ny f de J dz sin x os;(x) [ der exp (— c1?) - 
- {Pil€) + by(€1) — Pile’) — Pi(er')}. 


The motion of particles 7 and 7 with velocities c,¢, before collision 
can be separated into: 
a) motion of centre of gravity: mass m;+mj;, velocity 


(14.9)* 1 


(14.10a) go= e+ Prrer; 
b) relative motion: mass mij (14.4), velocity 
(14.10b) §=—C+ Ci. 
Here we have introduced: 
(14.11) B={mil(mit+m)}? 5 Bi={mj/(mit+my)}* (62+ Bi2=1). 
After collision, 7 and j have velocities ¢’, «1; go has not changed, but 
(14.10c) g’=—c’+¢,’ with g’=g. 


In the course of the following derivation, we shall temporarily introduce 


1) * after the equation number indicates that all occurring velocities are 


reduced ones. 
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a fixed velocity vector v; velocities from which v has been subtracted are 
indicated by corresponding capital letters: 
(14.12) C=c—v,etc.; Go=go—v 
Earlier we used this notation to denote peculiar velocities, but there is 
no danger of confusion because v will disappear from eq. (14.42) onwards 
and the peculiar velocities will not return until eq. (14.46b). 

C and C’ are reduced as in (14.8a), Cy and C;’ as in (14.8b), whereas: 
(14.13a) 80" /go = Go*/Go =v* /v = {(mi + mj;)/2kT }S; 
(14.13b) g*/g=g'*/g' =(uij/2kT)}. 
The following relations hold between the reduced velocities, instead of 
(14: 10, 12): 
(14.14)* go=fe+fia; g=—fic+fa; g’=—fie’+Pa'; g'=g; 
(14,15)* CO=cO—fv; G%=a%—-fiy; Go=go—v. 

The rotation of the g vector !), which takes place in the course of a 
collision, is described by 
(14.16)* g’—gcos 7+hsin y cose+k sin x sin e, 


in which g, h and k are three mutually orthogonal vectors of equal length. 
The direction for which the azimuth e=0 can be chosen arbitrarily in 
the plane | g; a choice on the intersection of this plane with the plane 
of g and Gp is convenient, because now 
(14,17)* Go-k=0; Go-h=+|Go A gl/=+|C A Cj}. 
We use these relations to derive expressions for C’2 and C;’2 in terms of 
C and C,. First, C’2 is expressed in Gp and g’ (see 14: 14, 15), then g’ is 
expressed in g, h and k (14.16): 

Ae = (Go — Big’)? = B2G'o2 — 2881 Go: (g cos y +h sin x cos €+ 
(14.18)* Peale | 
( +k sin y sin e) + i292. 


Now we return to C and C, by means of (14.17) and 
\ Gio? = B2C2 + 2681 C- Cy + Pr?C 712; 
g= fi?C2 — 268, C. Cc; + B2Cy?; aie 


(14.19)* 
( Go. g=—fpr0r-+( + (BP — Br?) C-Cy + parr. 


This leads to: 
(14.20a)* C%— PC2+ P, 24+ 2Q0C-C, + 2R\C A Ci| cos e; 


> 


(14.20b)*  Oy2= P1024 PO2—2QC.C, + 2RIC A Gi] cose, 


1) From now on up to eq. (14.45), all symbols in the text are also supposed to 
be reduced. 
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in which we have used the abbreviations: 
(14.21) \ P=1— 46262 cos? y ; P, = 46262 cos? y; 
‘ ( Q = 2661(62 — B12) cos? y; h=286) sin p cos y. 


Here y is used instead of y: 


(14.22) y+ 2y—7. 
Evidently, 
(14.23) P+Py\=1; PP\=Q?+ R?. 


We quote the following integrals for future reference: 
(14.24a) f dx exp (— Ax? + 2B-x)=7°sA-*s exp (B2/A); 


(14.24b) ( f de J dx exp (— Az? + 2B-x + 2|D A x| cos e)= 


( = 27° A~*h exp {(B2+ D2)/A}-Io(2|B A D|/A). 


Here Jo(z) is the zero-order Bessel function of the imaginary argument 72. 


Now we are ready for the determination of the eigenvalues and 
eigenfunctions of Q;; as given by (14.9). To this end, 2;;© is applied to 
the generating function of the Sonine polynomials with argument C?, 


which reads: 


(14,25)* Fs) == (ls) exp (oC), 
in which 
(14.26) o = s/(1—s). 


Omitting the numerical factor (1—s)~'~** for the time being, we obtain: 


(14.27)* 21; exp (—oC?) = ah n; J dy sin x aij(x) (+1 —1 —yv); 


of the four integrals s, the last three can be handled by using (14.24a, b): 


(14.27a)* 4 =f def dey exp (—c1?—oC?)= 2x exp (—o0?); 


= 20 { dey exp {—(1+.0) ¢1? + 20f1 v-¢,—of12v?\ = 


= Qq"2(1 — 8)" exp (— shi? v); 


\ a= (del dcezp(—a—ol7)= 
(14.27b)* 


7 =( del dC exp (— gC.) 
=f de f dC, exp {—(CGit+fiv)? 


_o(PC2+ P1012 + 2QC- C1 + 2R|C A G;| cos e)}= 


(14.27¢)* = On lf(1 —s)/(1 Seay a exp {s(1 —sP)-1 : 


-(PC2 + P; Bi2v? + 2Qfi1C-v)}- 


Ip{s(1 —sP)7-2Rfil€ Avi}; 
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J =f de f dC, exp (—c2—o0}'2)= 
\ = 2ar's{(1 —a)/(1—sP3)}"* exp {s(1—sP;)- 


| . (PC? + PB:2v2 — 2Qf, C-v)}- 
\ -Io fs(1 —sP;)1-2Rpi|C /\ vi}. 


Now we return from C to ¢ by means of (14.15); if x denotes the angle 
between c and v, 


(14.28)* ¢-v=cu cos x; |C Av|=|e Av|=cv sin x. 


This gives: 
24; exp {—o(c?—2fev cos x + f2v?)}= 
= 2an,; J dy sin x xij(x)- 
- [exp {—o(c?— 2fcv cos x + B2v?)} + (1 —s)exp (— sf 2v2) — 
—{(1=s)|(1—sP )y"*. 
-exp {—sP (1—sP )-1c24 28(1—sP )-1. 
(14.29)* -(1— 261? cos? yp) Bev cos x—V v?}. 
-Lo{ 4s (1—sP )-! Bf)? sin yp cos y cv sin x} — 
—{(1=s)/(1—sP1)}- 
‘exp {—sP(1—sP})1c?+ 4s(1 —sP,). 


-BB1? cos? p-cvcosx — Vyv 


2 


Ais 
5 
‘Io{ 48 (1—sP1) Bf? sin y cos y cv sin x}], 


in which V and V; are coefficients which we do not write down explicitly, 
because they will not be needed. We write the x-dependent factors in the 
terms of the expression between square brackets | ] as follows: 


\ K(x) “exp (0 cos x); 
(14.30) ‘K'(xn) & exp (0’ sin 6 cos x)-Io(0’ cos 6 sin x); 
Ky'(x) & exp (@1’ cos p cos x)-1o(@1' sin y sin x), 
in which 
(14.31) tg 0=(1— 2f1? cos? y)/(2f12 sin yw cos yp); 
g y y ¥ 


\ G = 2opev; 
(14.32)* oO = 2a —sP)-1 BevPt; 
Oy’ = 2s(1 —sP;)-1 pi cuP;}. 


The three functions K(x) shall now be developed in Legendre polynoms 
Pi(cos x). As for the first one, 


(14.33) K(x) =exp (0 cos x)= > Pi(cos x)(2U+ 1) i! jx —1Q), 
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in which 7; is a Bessel function of half-odd order: 
(14.34) juz) = (o0/2z)T144(2). 


As for K’ and ky’, Jo came in during the integration (14.24b). We go 
back one step and put 


(14.35) [o(z)=(22)44 f dd exp (z cos A); 

0 
this gives 
(14.36) Ky) =a) i di exp (@1’' cos &), 

0 
in which & is the angle between two vectors with polar coordinates 
(0,0) and (x, 4). Now we develop exp (@;’ cos é) in P; (cos &) ag in 
(14.33), then develop P; (cos &) in sperical harmonies of (0, 0) and (x, 4), 
and finally integrate over j: 


(14.37b) ky aS Pi(cos x) Pi(cos w) (21+ 1)i'71(—i@1'); 


(14.37a) K(x) = ¥ Pi(cos x) Pi(sin 6) (21+ 1) i*j( — i’). 
=0 
Substitution of (14: 30, 33, 37) into (14.29) gives: 
(14.38a)* ¥ Li(21+1) Pi(cos x) =0, 
1=0 


in which JZ; is the following operator: 
in = Q2,© [exp {—a(c2-- 620") k:(O)} | — 
—2nun; f dy sin x x«;(z) - [exp {—o(c2 + B2v a ie )}+ 
+(1—s)"* exp (—sf12v?) di9 — 

(14.38b)* {(1—s)/(1—sP )}* exp {—sP (l—sP )-1c2—V v2}. 
-P(sin 0) ky(O" )— 

{(1—s)/(1—sPi)}* exp {—sPi(1—sPi)1¢e?— V0}. 
-Pi(cos w) ki(O1’) J. 


| 
| 


Here it is understood that 2;; acts on [ ](2/+1) Pi(cos x); k; is an 
abbreviation of: 


(14.39) ki(z) & 4 jy( —iz). 

Now we develop Pi(cos x) in spherical harmonics Y{;) of the directions 
of ¢ and v, which we call (9g) and (9’y’); next we multiply by Yum(0'¢’), 
and then integrate over all directions ()’9’). It appears that (14.38a) 
is narrowed down to read: 


(14.40)* L1Y im(Pp) = 


for each separate J,m. Each term of L; contains a factor k;(const. v); 
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ki(z) shares with 7,(z) the property that a series development in powers 
of z starts with 2z!: 
(14.41) key(z) = 24 1! {(214+1)!}-1244 O(z!41), 
Using this property, we divide (14.40) by 7! 2!v!'/(21+1)! and then take 
the limit v > 0: 

243 {exp (—ac?)-(20fc)' Yim(dg)} = 
est ion oij(~) Yim(Y@) - [exp (—ac?)-(2afe)! + (1 —8)*d19 — 
—{(1—s)/(l1—sP )}* exp {—sP (1—sP )-1c?} 

-Pi(sin ) {2s(] “SP )7BeP? }— 
— {(1—s)/(1—sP))}"* exp {—sP\(1—sP)—c?}. 
-Pi(cos p){28(1—sP1)— Bi cP1}]. 


(14.42)* 


Divide by (2s)'(1—s)"* and introduce the generating function SY according 
to (14.25): 


Wei {(Bc)! Si44(c23 8) Yim(Op)} = 
(14.43)* — « = 2a; J dz sin x xij(z) (Bc)! Yim(dq)- 
{Pi44(C?; 8) +000 — S144 (c2; sP) P#! P; (sin 0) — 
— F144(c?; 8P1) Pi*(p1/f)! Pi(cos yp). 


Differentiate r times with respect to s and let s—+ 0: the -functions 
render Sonine polynomials according to 


(14.44) SiQa(0) = lim (=) S 142; 8), 
and we get the eigenvalue equation for ¢)—0: 
\ $j {(Be)' S12 4(C?) Yim(Op)} = 2an; J dz sin x oi(z)- 
(14.45)* {1+ 610 6r0 — Pr! P,(sin 6) — Py", /B)" Pi(cos y)}- 
“{(Be)' S19 4(c?) Yim(O—)}. 
Now we give the results of the following four steps : 
(a) return to natural values of ¢ instead of reduced ones; 


(b) derive the result for the general Q;; (for e+ 0) from 24; (for 
co=9), in the way indicated in (14.7); 


(c) introduce the symbol C in the original meaning of ¢—co; 
(d) multiply both sides by {2k7(mj+my,)}*. 


This gives the general eigenvalue equations: 
(14.46) Qi; yy (imr) = wi) yiimr) ; 


(14.46a) wi = an; f dy sin x 4;(x)- 


. {1 + 61000 — Pr+#! Pi(sin 6)— Pi" +41(B1/B)! Pi(cos yp); 
(14.46b) gn. od, gi) = (miC)! ST} 4(m¢C2/2kT) Yim(I¢). 
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We have not changed the notation (Jp), which was used in (14: 40-45) 
to denote the polar angles of ¢; in (14.46b), of course, polar angles of 
C are meant. 

As far as the C-dependence of the eigenfunctions 4; goes, they 
are identical to the 7" for the single gas case (eqs. 7: 3, 5). Of course, 
the difference lies in the fact that, in the eigenfunctions for N>1, m; 
occurs and the operator 2;; works on both the C- and the 7-dependence 
of y. 

In the derivation of (14.46), all RHS contain four terms corresponding 
to the four terms: 4;, dj, —¢:’, —¢j’ in the definition (lae2) ol Qin 
If we had accordingly split up Q;; into four parts from the beginning: 


216 = Jde f dy sin x f der fj(€1) oaj(z) ile) 

(14.47) pap = Pi(¢1) 
Qnd = di(cr’) 
pa z piler’), 


we would have found four separate eigenvalue equations: 
(14.48a, b, c, d) 921, 11,101,1v 4(€) = or, 11,1, 1v xi(¢), 


in which oy,...,ry, in the given order, are the four terms of the eigenvalue 
oi” (14.46a) and x; is the eigenfunction given in (14.46b). The first 
two equations (14.48a, b) are trivial (at least in the case of a Maxwellian 
mixture, in which « does not depend on g), for Q; is a mere multiplication 
and 2 applied to any of the functions (14.46b), apart from ;(0% — 1, 
yields zero on account of the orthogonality properties of the Sonine 
polynomials and those of the spherical harmonies. 

This also means that we have solved the eigenvalue problems of the 


operators 

(14.49) Yy =24+0m; Zy = On+Qy. 

The solutions read: 

(14.50) Vy a hie) = yi) yyhim) ; Lij yg bmn) = 24) nr), 

Here, y:(’"” are the functions given in (14.46b); the eigenvalues y;; and 
zi; together yield mj as of (14.46a): 

(yg) = 2a; J dy sin x ois(x)- 

eee) ( -{1—(1—4f? Bi? cos? p)t+#"P;( [1 — 261? cos? p][1 — 46261? cos? y]-*)} ; 
(14.50b) 217) = —2nn,; f dy sin x xaj(x) (26/1 cos p)?"+! Pi(cos yw); 

(14.51) yay) + 245) = wy"), 

Moreover, the 7; are also eigenfunctions of the W,; operator defined in 
(10.24): 

(14.52a, b) Wag bmn) = wig) 44.0 5 wig) = 245) + Og D Yew. 


k 


27 Series B 
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In the multi-component case considered, simple recursion formulae 
between the eigenvalues, like those holding for the single Maxwell gas, 
are not obvious. Here we shall mention only those eigenvalues y, z and 
eigenfunctions y, which correspond to an eigenvalue mjij=0: 


(a) 1=0,r=0: 7, =1; yj 0) = — 24,00) — 0, 


(b) [= if r—O(: iy A900) — mMCYim(Iq); Yiz) — — 24; (10) == 


(14.53) ( = 2an; Bi? J dy sin x «i;(y)(1—cos x). 
(c) I=0,r=1: x, =§— ae yi OD) = — 240) = 


= 4an; B? Bi? f dy sin x x4(z) (1—cos x). 


We have found only five eigenfunctions 7; with w;;=0, whereas we stated 
in (10.5) that there are NV +4 collional invariants 7;. This is because the 
method used above generates functions 7 of C which only depend 
parametrically on 7. On the other hand, for the derivation of the NV 
conservation laws for the mass densities 9; (10.13a), we have, instead of 
the one function 7; =1 of (14.53a), used the N collisional invariants 
yl -K] = bi~mx (10.5a), which are all eigenfunctions of both Yi; and Zi; 
at the eigenvalue yijj=2;;=0. This last property is shared by each 
function y;(m1, ..., my) of all N masses, and conservation laws could be 
formulated for the moments of all these functions, but the conservation 
laws (10.13a) for the moments of the 7;!!-4] are of special importance. 
We may remark that, for the /-r-combinations (00), (10) and (01), 
corresponding to wij =0, the eigenvalues of Wj; take an especially 
simple form: 
\ t=5: Wye bm Yir= — >’ 2x; 


cei kei 


(14.54) 
( b5e9: Wij = Zig = — Yij. 


The yi” as defined in (14.46b) are mutually orthogonal but not 
normalized. Normalization is of course possible, but it should be noticed 
that, if we use normalized eigenfunctions: 


(14.55) Xi lmMr) = ogy (lr) y,(Umr) | 

in which 

(14.56) {org (ir) 2 = (i7)() = (yalemr*, yim); , 

some of the corresponding eigenvalues are different from those given in 


(14: 46a, 50a, 50b, 52b). This is seen by looking at the linearity properties 
(10.25): 


( Yap xi=yer xi > Yigloi yi) =yi ai x3 
(Zig x4 = 24g 14 —> Zagloxg 43) = (cog | ora) + 245 4 H0- 


Evidently, when we replace the eigenfunction y;(””, with eigenvalues 
(y, 2, w, @)y™, by the normalized eigenfunction (14.55), we get eigen- 
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values, some of which are now different: 


4 Fg) = yg : Zig) = (06307 fxg) 24,0) ; 

14.58 = = See = eve = 

oe, Oi" = Jag Zig 5 De) =z) + by, > Fie. 
k 

In ch. 16, where the normalization of the x plays a role, we do not use 

these 7 and @, etc.; the old values will be used and (i7)() is introduced 

when necessary. 


15. Transport properties 

We shall now indicate how the eigenvalue theory given in ch. 14 ig 
put to use for determining the transport properties of the Maxwellian 
mixture by the procedure outlined in ch. 10. 

We recall that the relative deviation of /;,°2 from a local MB distribution 
fi is given by 4; as of (10.36); the coefficients A occurring in this ¢ 
follow from the integral equations (10.39a, b, c) supplemented by the 
normalizing conditions (10.41a, b). The special advantage of considering 
a Maxwellian mixture lies in the fact that the functions G, which are 
the LHS of the eqs. (10.39), are found among the y;, the eigenfunctions 
of the operators Qi, Yij, Zi; and Wi; of the Maxwellian gas mixture: 


(15.1) Gi) > 70m); GH) > yql2m) 5 GI) —> y,(1m0), 


In each case x, the — sign indicates that the independent components of 
G, are linear combinations of the 2/+ 1 degenerate y;(’") (J and r fixed). 
Therefore the exact solutions to each eq. (10.39) will be given by its 
LHS, G, times an i-dependent coefficient: 


(0) 
(15.2a, b, c) A; =a, (2—O;,*2) mC; A; =a; C,* C,*; AyD =a, mC, 


in which the coefficients a are found by substituting (15.2) back into 
(10.39) and (10.41). The resulting equations are brought in the nicest 
form by using W; instead of 2;, so that the linearity property (10.25a) 
of the W operator can be used: 


(15.3a) (G—O1*?) C= Wifa (5 —C,*2) C} = (@—C*2) ma CS wis a), 
or I/m= Y wy aj”: 
(15.3b) 26," C,* =Wia C°C;*) =O," C;* Yay aj, 
or ~2= YS wy 2 aj, 
(15.3c) ( (1—(Sie/2)} C= Wi(a me) =m ES wy a,®), 


t] 
or (1/mi){1 = (Oix/a) } = > wy a; *) : 
j 


(15.34) 0=(a'®) mC, mC), 
or O= Sa; mja;), 
j 
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The normalization of AM (10.41a) does not give rise to a separate equation, 
because A® as given by (15.2a) is already orthogonal to mC. The 
normalization of A® (10.41b) is expressed by (15.3d); this is a useful 
supplement to the eqs. (15.3c), because the matrix w;j" is of rank V—1 
(cf. eq. 14.54). 

The transport coefficients are expressed by eqs. (10.56) in terms of 
G and A; substitute G and A from (10.39; 15.2): 


(a) w= kT > (—26;* C;*, a; C;* C;*); = —43nkT > x; 4;" 
ry ri 
(b) A he Lk > (ls 2 -C; i*2] C, aj; [3 —C;*2] m;C),; = 3nk2?T ¥ ajay; 
¢ j : j 
(c) Da=(1/3n) > ([1— (dey/x2)] C, ay [5 —Cy*2] m; C); = 0; 
y 


B)- 


(d) Diyx=(1/3n) > ([1—(dey/22)] C, ay mj C) = kT(> aja; —a,®)). 


7 
The thermal diffusion coefficients are all zero; according to (10.61), 
A and 4’ coincide. 

Eqs. (15: 3,4) give the framework for an exact calculation of the 
transport coefficients of a Maxwellian mixture. They are analogous to 
the approximate expressions (13: 34, 35,36) for the general potential 
field. The eqs. of ch. 13 contained matrix elements of the Wi; operator 
between eigenfunctions of the Maxwellian 2;; and Wj;; instead of these, 
the eqs. of this chapter, of course, contain the eigenvalues w;;. In both 
chapters, the coefficients a and the transport coefficients (rt, or ju, 2, D) 
have similar meanings. 

To get explicit values of the transport coefficients, the w;; values given 
by (14:50,52b) should be inserted in (15.3). We shall illustrate this on 
the simple case of binary diffusion. First we introduce 


(15.5) a= my/(m + m2) 3 My = Mme/(mM1 + me) : M12 =™,Ms/(m1 + me); 
(15.6) r= o/{n(m + me)} =a mM, + Xe My. 


From (14:54, 53b) we see that the elements of the w;; matrix for 
N=2,l—1,r=0 can all be expressed in one quantity containing the 
1-2 differential cross-section only: 


(15.7) w 27 anus J dy sin x x12(x%) (1 — cos x), 
as follows: 
W111 = Y1o= WL2/ My, W12= — Y12 = —WXe/™M 
(15.8) \ / y 2/my 
( Wa = —Y21 = — way/mMg W2= Y= wx/Me. 


For k=1, the two equations (15.8¢), which are identical and (15.3d) can 
be written as follows: 


( 1/aiw= 7,0). yt 


(15.9a) ” 


0= Xy M, a1) + 2 My ag), 
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with the solution: 
(15.10a) aq = —x u,/a,rw; a2) = om, /rw. 


Similarly, for k=2, 


(I/azw= a, — Ag?) 
(15.9b) 5 2 

( O= 2 Mm, a) + a9 Mag) ; 
(15.10b) a) = u,/rw; ag) = — a m,/xorw. 


Substituting (15.10a, b) into (15.2c) and thence into (10.36), we conclude 
that, for a binary Maxwellian mixture with a concentration oradient 
Vx=Vx1= —Vx2 only, the distribution gi reads: 


(£5711) Sie Mah git ee Dy fe 


rw Hial Tw v2 
According to (10.76d), the binary diffusion coefficient reads: 
(15.12) D=2 a2 Dy =212 (Dir — Dip — Doi + Dog). 
Substitute the Diz as given by (15.4d): 
(15.13a) D=kT x1 %9(— a, 4+ a, ® 4 ag —az@)), 


or, inserting the a-coefficients from (15.10a, b) and using (15.6), 
(15.13b) | D=kT |w=kT]/{2am12 f dy sin x o12(7)(1 — cos x)}. 


We compare this result to the usual ChE formulae for the first approxi- 
mation to D: 


(15.14) [D ]i = 3kT /(16 npt2£219(1+1)) 
with 


(15.15) QA D= zt i e-9*2 g*4dq* { dy sin yaie(g*, ¥)(1—cos x) 
0 i) 


(g* as defined in (14.13b)). Introducing the property of the Maxwell 
interaction: a2 independent of g*, we get: 


(15.16) Qy20)) = zt - 2 - (w/ 2712) 


and (15:14,15,16) together yield (15.13b): the ChK first approximation 
leads, for a binary Maxwellian mixture, to the exact value of the diffusion 
coefficient, as it should. 

If we substitute the A coefficients of (15.2) back into (10.36), the 
resulting ¢;%) is an exact solution of the integral equation #;® =02;6 
(10.21b, r=1). Thus, for a Maxwellian mixture, the formal developments 
of F;® and 4; in eigenfunctions of 2 (cf. the developments (2: 46, 47) 
for the N=1 case) break off very soon. 
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16. Relaxation of the velocity distribution 


In the preceding chapter it was shown that, for a Maxwellian mixture, 
the development of the first non-MB term in the ChE velocity distribution, 
di", in eigenfunctions of 2;; consists of a few terms only. We now turn 
to the equations derived in ch. 11 for the non-ChE or kinetic terms of 
poe bi (r= Leip eats 

a) integral equations (11.9); 

b) initial conditions (11.10); 
‘) 


fe) 


normalizing conditions (11.11); 

d) normalization of the initial condition (11.6). 
The formal solution is again based on development of the ¢; in 
eigenfunctions, y;/""), of 2;;. As has been mentioned at the end of ch. 14, 
we shall use 7 which are mutually orthogonal but not normalized: 


(16.1) (ygltmr* re ide) =Oimr, Um’ r'(t1)M 


(cf. eq. 14.56). In the following, we shall suppose that all 7; are real (to 

this end, spherical harmonics with real g~-dependent parts are used); 

now the complex-conjugate sign * can be omitted from all scalar products. 
We start with the first equation, (11.9a); develop 4;™ in ;: 


(16.2) bi) = a zy llnr) gy tmr) 
lmr 
and 
(16.3) bi Ut =) d;(0) = > iim py (lmr) (0), 
Imr 


in which the coefficients 


(16.4) pir (0) = (;(0), yim) /(ii)—) 


are known. Substituting (16.2) into the equations (11.9a), we find, for 
each Imr, N coupled first-order differential equations for g;(m"): 


> Pg she! es 
(16.5) ot py (tmr) ss $3 wy) (mr), 
The solutions read, of course: 


k 


putting this back in (16.5) and equalizing coefficients of each exp (—wy,ft), 
we find: 


(16.7a) Pine UMP) wy llr) => win Dix lmr) 
or 
(16.7b) > (wij) — dep We) Djglmn) = 0: 


for fixed Ir, each wz occurring in the exponents of (16.6) is an eigenvalue 
of the w;j')-matrix, and the eqs. (16.7) determine the components 
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pik(t=1,...,.N) of the corresponding eigenvector, together with the 
N initial conditions: 


(16.8) gi llmr(0) = S pine, 
k 


where the LHS are given by (16.4). 

The key role of inverse relaxation times is played by the eigenvalues, 
we”), of the wz”) matrix. There are »(< N) of these; we number them by 
k(=0,...,.y—1). If wz has a multiplicity mz(S.m,~=N), the N equations 
(16.7b) for fixed Imrk at i=1,...,N have a coefficient matrix of 
rank N—m,. So the number of independent equations supplied by eqs. 
(16: 7, 8) amounts to 

=i yl 
(16.9) > (N-m:)+N=Nv— Ym+N=N>, 
k=0 k=0 
which equals the number of coefficients qj to be determined. 

We now recall that the matrix element [F, F] = (Ff, WF) is positive 
unless F is a collisional invariant, when it is 0 (eq. 10.30); if we put 
Fy=7,%™), we find: 


(16.10) >» > (20) ) wy) > 0. 
ie 


It follows that, for general 7, the w;;-matrix is positive definite and 
all its eigenvalues w,>0. An exception is formed by the collisional 
invariants. For 7,©%—1, all w4©9—0; for y,0™ and 7;©, the wy;- 
matrix is given by eq. (14.54): 


211 — > 21j =12 “1N 
j 
221 222, — > 22) 
(16.11) (wij) = 4 
ZN1 ZN2 2B) 8 aN > ZNj 
: j 


Clearly det w=0, as is seen by adding all columns. The matrix for NV 
components, to be indicated by the symbol w(N), is positive semi-definite 
and of rank N—1. This is proved by induction as follows: for N=1, w(1) 
consists of one element 0. Suppose w(2) to be positive semi-definite of 
rank 1; ...; w(N—1) to be positive semi-definite of rank N—2. For each 
of these matrices, this implies that all principal submatrices are positive 
definite and the corresponding determinants (principal minors) are > 0. 

Now consider the N—1 principal submatrices of w(N) found by taking 
only the elements with i,j<n(n=1,...,N—1) and call them p(n). Such 
a p(n) differs from w(n) in that its diagonal elements are greater by an 


amount 


(16.12) piu(n) —wi(n) = — > = aun), 
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which is >0 because all a4 are >0 and 1—cosy¥>0 (see 14.53b, c). 
If one adds one positive term gi(n) to the ii-element of w(n), det w(n) is 
increased by an amount qii(n)mi(n), where mu(n) is the minor of the 
vi-element in w(n). Both gq and m are >0, so the replacement of one 
diagonal element in the positive semi-definite matrix w(n) by the corres- 
ponding element of p(n) already yields a positive definite matrix, and 
continuing in this way we find a complete p(n) which is certainly positive 
definite. 

Evidently, det w(N)=0 and all its principal minors of increasing size: 
det p(n) with n=1,..., N—1 are >0. This proves that w(N) is a positive 
semi-definite matrix of rank N—1, Q.E.D. 

We conclude that the inverse relaxation times of (16.6) have the 
following signs: 

( Ir=10, 01: wo =0; wi, ..., wi >0; 


(16.13) (i+ r = : wo), oe Ww, 6) >(), 


The only terms in (16.6) which do not die out in time are the constant 
ones at Jr=10, 01; k=0. These, however, are ruled out by the normali- 
zation of ¢;, as can be shown in the following way. 

Apply (11.11) to the solution ox) given by (16: 2, 6): 


~ 


(16,14) (6G), ytimr)) = 5 (Gi) LMS igtmn exp (—w,t)=0, (lmr=1m0, 001) 
i K 

which for each separate k value yields: 

(16.15) > (ii) Or) Hig (tmr) =) (lmr = 1m0, 001). 

For each of the four /mr values concerned, this gives » additional condi- 

tions on the Pik, but we have already argued (see 16. 9) that these Ny 

coefhicients are fully determined by the eqs. (16: 7,8). Now it can be 

shown, to begin with, that the eqs. (16.15) for k= 1,...,»—1 are dependent 

on the set (16.7). To this end, we multiply each eq. “(16.7a) by (t)( and 


take >;: 

\ WKEM SY (41) UY iglm = Y jx ltmry > wis") (ti) On = 

| ge ¥S pjx(tmr) > Wij — 0, 
| i 


in which 


(16.16) 


(16.17) Wi — (ya (emr), Wij y(imr)), — wy” (vi) Cr), 


These Wi; can be expressed in terms of Zi in the same Way as we 
wrote wij" in terms of 24;(”) (16.11), viz. 


(16.18) Wig) = Zig — diy > Zin, (r= 10, OL) 
k 
It follows that 


(16.19) > Wit) = iS Zn y) Zi") 0, 
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which vanishes on account of the symmetry of the Z matrix element 
(see 10.29b). Consequently, (16.16) gives: 
(16.20) WE SF (GUM Diy lmr) =O : (lmr = 10, 001) 
i 

the pix with lmr=1m0, 001 and w,~0 (v.e., k>1) as determined by 
(16: 7, 8) do satisfy the normalizing condition (16.15). 

In order to prove (16.15) also for k=0, we=0, we make use of the 
normalization (11.6) of the initial values ¢,(0), which yields the equi- 
valent of (16.14) at f=0: 


(16.21) > (a) Pi Gmr)(O)= > > (2) pink”) =). (lmr=1m0, 001). 
a k 4% 


Comparing this result with (16.15), which, as we have seen, is satisfied for 
k=1,...,y—1, we find that (16.15) does hold for k=0: 
(16.15a) > (41) pio) ead), (lmr = 1m0, 001). 

v 
This equation is not dependent upon the N equations (16.7) with 
k=0, wr=0, which read: 
(16.7c) 0= > we pol. (mr = 1m0, 001). 

j 

The matrix of these equations is of rank N—1, so (16.15a) and (16.7c) 
form NV independent homogeneous eqs. which have the zero solution only: 


(16.22) piol™ =0. (lmr = 1m0, 001). 


This implies that e-powers with exponent zero, corresponding to infinite 
relaxation times (lr=10, 01; k=0; w;() =0) do not occur in (16.6). 

We sum up the preceding argument: the requirements (a, b, d) mentioned 
at the beginning of this chapter lead to eqs. (16: 7, 8, 21), from which 
the Pik can be determined. These gj, automatically satisfy the conditions 
(16.15) derived from requirement (c). In particular, coefficients corres- 
ponding to non-negative exponents vanish, which enables us to write, 
combining eqs. (16: 2, 6): 

(16.23) biD = ¥ yom pa Pik exp (— wyl£) (w> 0). 
ime 

It is possible to extend this treatment to higher approximations: 
¢;2), ete., like we did for a single gas in ch. 3. The occurring relaxation 
times will of course be sums of the positive wz, the eigenvalues of the 
wij matrices. We will not do this explicitly. 

To conclude this chapter, we shall now illustrate the r=1 theory by 
the following simple example: in a binary mixture, an initial deviation 
of the ChE distribution is given which, at ¢=0, consists of a term 
proportional to 7; =m, C only, and which satisfies the normalization 
(16.21): 


(16.24) r(0)= (cpa); bx9(0) = — (pla) € 
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(~p is a constant vector). The w;j-matrix for the N=2 case is quoted in 
(15.8); its eigenvalues tollow from 


(ware/m1) — We —wee/my 
(16.25) . ay 
—way/ms (way/me2) — Wr 
which gives: 
(16.25a) wo = 0; Wy =rw/ re = W(x1mM, + L2Mz2)/M 1 Me. 


(Definitions of r and w as given in (15: 6, 7)). 
It is now easy to find the solution 4, of the form (16.23): 


(16.26) bi) = (ep/ay)-Ce-mt ; pg) = — (paz) Cet. 


~ 


One can show that this 4“) does in fact satisfy the equations (11.9a). 


17. Quasi-relaxation and the kinetic time scale 


We shall now discuss the ‘‘quasi-relaxation”’ of the diffusion current, 
as we did for the transport of heat and momentum in ch. 6. The discussion 
will concern a binary mixture in both the kinetic and HD stages. In 
either case, ¢ will be used as a symbol for the relative deviation from a 
local MB distribution, like in the entropy treatment (12.3). 

From eqs. (10: 50,51), we see that the average velocities (y>i,2 of 
components | and 2 relative to the mean particle velocity u are given by 


(1751) nkT <y>j=(H®, 4), 
in which we have replaced G;‘) by 
(17.2) H,o) = —kTG;9 =kT {(64;/2%) —] Ci 


like in (6.2). Subtracting the two eqs. (17.1), we get the following expression 
for nkT’ times the diffusion current: 


(17.3) Jo=nkT (<C>: — (C2) =(H® — H@), J) = (H, ¢), 
in which 
(17.4) H = HO) —H®), i.e. Hy=kTClay; He= —kTC/xo. 


Now we calculate the collisional part of Jp, introducing the approxi- 
mation ¢? <1: 
\ IcJo= an de Hdcfi= > SY f de Hid (fi fi1) 

i : ¥ 


/ ww — > Dd J defi WH; Q5¢=—[H, 4]. 

\ 3 

Comparison of (17:3, 5) gives a result equivalent to (6.12b): 
(17.6) dc(H, ¢)= —[H, 4]. 


Like in ch. 6, we can develop H and J in eigenfunctions y; of Q;;, and 
equations analogous to (6.15) will show that each Imr-mode of Jp decays 
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with relaxation times which are the inverse of the eigenvalues of the 
wij -matrix. We shall demonstrate this for the HD distribution @ in 
a binary mixture with a concentration gradient Vx only; this ¢ is given 
by eq. (15.11). We use it to calculate Jp and dcJp from (17: 3, 5): 


oe » kT 
Were ie ae 
( ) nkT 1%. w 


in agreement with D=kT/w, (15.13), and 


ae WcJb ya r 
TTS eee lL —; 
( ) nNkT H1 Xe eee 
SO 
(17.8) dcJo/Jo= —rw/m2= — wy, 


as defined in (16.25a). Like the corresponding equations for uw and A 
(7.32), this implies that collisions try to break down the transport (in 
this case, the transport of particles), a process with a characteristic time 
I/w;. In the HD stage, of course, this tendency is exactly compensated 
by the streaming terms, (dg+dq) Jp. 

We now compare the following expressions, which we have derived 
for a binary Maxwellian mixture: 

a) distribution for the HD stage in the case of a gradient Vx: (15.11); 

b) distribution for the kinetic stage in the case of an initial 4; proport- 
ional to m,C: (16.26); 

c) quasi-relaxation of the particle current: (17.8). 
In all three equations, we encounter the following characteristic time: 


17.9 — 1 pie  mi+me 1 
(17.9) ED = wi Tw 0 22 J dy sin yo12(x) (1—cos x) * 


Introducing this tp, we write the expressions (15.11; 16.26; 17.8) as 
follows: 


(17.10a) M=—tpC-Vaji1; da=tpC-Va]xe; 
(17.10b) di=p-Ce-n/a, 5 do= —p-Ce-"*n/x9; 
(17.10c) IcJo/Jo= —I1/tp. 


Apparently, tp plays the role of 

a) coefficient in the diffusional part of the HD distribution; 

b) relaxation time in the diffusional part of the kinetic distribution; 

ec) quasi-relaxation time for the diffusion current; 
it is clearly comparable to rt, and 71, defined in (7.13). 

From ch. 14 onward, we have only considered the special case of a 
Maxwellian mixture. As in the single component treatment, many results 
hoiding for a more general intermolecular potential in first approximation 
can be obtained from the exact result for the Maxwellian case by taking 
matrix elements of the general Q;;, in a representation of Maxwellian 7, 
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instead of eigenvalues of the Maxwellian 2;;. Arguing in this way, it is 
easy to see which quantity takes over the role of the characteristic 
relaxation time: expressing tp=1/w, in w by (16.25a) and w in 2)24,)) 
by (15.16), we find: 


3 7 3 


For a binary Maxwellian mixture, this is exact; for other potentials it 
holds in first approximation. (17: 9, 11) show that tp is inversely proport- 
ional to the mass density o, as well as to the temperature-averaged 1-2 
cross-section. 

Incidentally, the diffusion coefficient D can aiso be expressed in terms 


of tp, in the way (7.16) connects A to 1, and yw to T,: 


(17.12) D=kT |w=kTrtp/ ur, 


but this is not a very significant expression, because the x-dependent 
factors in r and in tp cancel. 
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List of symbols 


(Missing symbols are explained near to the place where they are used. 
Superseripts and subscripts are mentioned separately at the bottom. In 
the equations quoted, the symbol is defined, first mentioned, or best 
elucidated.) 


(10.36) A) coefficient in 4) 
(15.2) a; coefficient in 4; 
(13.30) ay”) coefficient in A; 
Cc Cc—Co 
(14.12) Cc c—v (ch. 14 only) 
(10.54¢) D; multi-component diffusion coefficient 
(10.54¢) Dr; multi-component thermal diffusion 
coefficient 
(10.34) d; generalized concentration gradient 
(10.39) Ge) function of C, LHS of integral 
equation 
(14.12) Ge go—v 
(14.10b) g velocity of relative motion: ¢,—¢ 
(14.10a) Zo centre of mass velocity 
(17.2) H(i) function of C ae 
(10.59) hry multi-component thermal diffusion 
ratio 
(15.5) M,, M, relative masses of I-2-pair 
N number of components in mixture 
(14.21) te tea functions of y or x 
15.6 i 
S Sonine parameter 
(10.46) u average particle velocity 
(14.12) v constant vector 
(10.24a) Wi integral operator 
(14.52) wy) eigenvalue of Wy4j;-operator 
(16.7) wy) eigenvalue of wij(”)-matrix 
ie Ww 
tae Li relative number density of comp. 7 
(10.18a) Vy integral operator 
(14.50a) yg eigenvalue of Ya 
(10.18b) Lij integral operator 
(14.50b) eu) eigenvalue of Zi 
(10.3) Xi differential cross-section = 
(14.11) B, Pr roots of relative masses of ¢-j-pair 
(10.47) Y particle velocity relative to u 
(10.3) E azimuth of g’ 
(14.31) 6 function of w or x 
(14.40) — polar angle of cor Cry 
(10.60) A heat conductivity coefficient in ab- 
sence of diffusion currents 
(10.54b) Us heat conductivity coefficient in ab- 
sence of concentration gradients 
(10.54a) lt viscosity coefficient 
(127 12,-13) [Mi thermodynamic potential 


(14.4) bay reduced mass of 7-7-pair 


(14.22) 
(10.9) 
(10.19a) 


(10.26a) 
(10.26b) 


(14.56) 
(10.27) 


(10.8a) 


(10: 68, 69) 
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number of different eigenvalues of 
wy") -matrix 

s/(1—s) 

diffusional relaxation time 

transport coefficient 

azimuth of ¢ or C 


coefficient in 4) 
coefficient in ¢@ 


coefficient in g,(mr 

angle between g and g’ 

eigenfunction of ;; for Maxwellian 
mixture 

eigenfunction of Qj; at eigenvalue 
zero (collisional invariant) 

& (%— x) 

moment of f 

integral operator 

average cross-section 

eigenvalue of Q;;-operator 

e.g., fi, a: indices of components 
+, 9,k 

e.g., we’: index of eigenvalues of 
wii) matrix 

e.g., yi, pi’: indices of eigen- 
values. of 2;;-operator for Max- 
wellian mixture 

e.g., X®,J@: index of transport 
phenomenon (=w, 4, 1, ..., N) 

e.g., C;*: reduced (dimensionless) 
value 

e.g., (Fi, yi): scalar product within 
component 7 

e.g., (Ff, 7): scalar product over all 
components 

norm of y;(/mr) 

e.g., [d, A]: matrix element of Q- 
(or W-) operator 

average of function Gi(e) over com- 
ponent 7 

average of function Gi(¢) over all 
components 

©.g., T: approximation derived by 
variational method (ch. 13) 

e.g., X, J, L: independent X and J 
and corresponding L-coefficients 

e.g., I’, di’: alternative definitions 
of forces and currents 


PHYSICS 


ON THE MECHANISMS UNDERLYING DE LANGE’S RESULTS 
BY 


F. VERINGA 


(Communicated by Prof. G. W. RarHenau at the meeting of January 28, 1961) 


In a previous paper [1] more evidence was offered for the statement, 
made by Dr Lance [2, 3, 4], that linear low-pass filter action occurs in 
the visual system prior to non-linear processes. Moreover the successful 
interpretation by LEvINsoN [5] of the experiments by BRowNn and ForsytH 
[6], based on Fourier-series, stresses the fact that the mechanism respon- 
sible for the fusion of on-off flicker behaves linearly. 

It is now intended to show that this low-pass filter action, found by 
Der LANGE [2, 3, 4], as well as by other authors [7, 8], can be explained 
on the basis of a mechanism located inside the retinal receptor cells them- 
selves. In doing so it is taken for granted that a photon absorbed by the 
visual pigment inside the rod or cone, produces, by photochemical reaction, 
a molecule (or ion, or radical, for that matter), which is specifically able 
to fire the receptor if and only if it reaches the wall of the cell. In that 
case the process bringing about the shifting of the molecule from the 
place where it is produced to the wall can be taken to be diffusion. On its 
way, however, the molecule may be rendered inactive by some reaction, 
which might be regenerative or recombinative. If the fraction of molecules 
lost by this reaction is important, i.e. if the quantum efficiency @ of the 
process described is low, or again, if the mean life of such a molecule is 
short as compared to the time it needs to diffuse to the wall, then 
clearly only molecules produced near the wall will have to be considered, 
and the shortest path to the wall will be favoured. 

Apart from the fact, proved by other authors [9, 10, 11], that the over- 
all quantum efficiency of the visual process is a few per cent only, it will 
be shown from Dr LancGe’s [4] data that the conditions described do 
apply. A reasonable model for calculating the mean behavior of a large 
number of receptor cells seems to be obtained by considering diffusion 
along a straight line from the place of production to a certain sensitive 
point, symbolizing the cell membrane, where the molecules under con- 
sideration are inactivated on arrival. For the sake of mathematical sim- 
plicity all photons are supposed to be absorbed at a distance xo from the 
wall, zo being smaller than the radius of a real receptor cell, perhaps con- 
siderably so. The value of the diffusion coefficient D for the kind of mole- 
cules expected (such as protein molecules) in the viscous cell mass will be 
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somewhere between 10-7 and 10-9 em2sec.-!: the characteristic time 


>) 
xo- é : 
— on of the process might be of the order of magnitude of 10-1 see. 


a 


It is now desired to calculate the flux of molecules at the wall. which is 


0 
D Ge °) when the origin is taken in the absorption plane, the x-axis 
Ov L=X0> = 
perpendicular to this plane and ¢ represents the time- and space-dependent 
value of the concentration of travelling molecules. 
The calculation has to be based on the diffusion equation 
»2 ) 


Cc gts 


oc + D = 
é oar ot 


o describing the rate at which molecules are rendered inactive on their 
way. The boundary conditions are 


: n) se ; 
Ce-2,=0 at the wall and D é | = 1+ Me'*' in the absorption plane 
z=0 


when the retina is illuminated with sinusoidally modulated light at a 
frequency of w/2a and a ripple ratio M. 
The solution is 


P) 1 M : 
D & | = 1 giot. 
ov Jez cosh V20t cosh aay 


where 


ax =\or + V m?t2 + 927? 4 iV — ot +) wt? + 9272. 


.. M cosh | 2ot 
So the ripple ratio on the wall is Sen 05 


. and the reciprocal of the 
cosh ax 


cosh ax 


complex filter factor Re” 
cosh |/ 957 


From this 


tany = tanh Vor +) wt? + 9?t?- tan ) ot +) w?t? + 9272, 


The amplitude of the ripple is diminished by a factor 


/4(cosh2 Vor +) w?t2 + 6212+ cos2 Ji ot + V w2r? + ot? ) 
ge een te capensis unt ae ee 


Re 
cosh | 2or 


This last result can at once be compared with Dr LANGr’s experimental 
points by plotting R as a function of wt on log-log paper, taking ot for 
parameter. The experimental points [4] for 0.375 troland are then nor- 
malized to give R=1 at w=0, inserted in the graph and shifted in the 


wt-direction until a reasonable fit is obtained. They are found to agree 
best with the curve for ot 10, 
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From this the quantum efficiency Q, evidently being is found 


8 cosh Vip 
to be some 2} per cent, thus satisfying the condition for validity of the 
calculation. The precision of the or- Seon is admittedly small, 
allowing a range of Q-values from some 7 per cent downward. The relation 
of this result to numerical values for the overall quantum efficiency 
will be an object of continued study. 

Moreover, comparison of frequency scales shows 1 —0.37 sec., yielding 
e=27sec.1, thus fulfilling the condition 1/o <t. 

When these values of @ and 7 are inserted in the formula for tany a 
prediction on phase-relation is obtained. 

The reason for selecting the experimental points at the lowest luminance 
available in Du LanGn’s [4] paper to determine parameters 0 and 7 is that 
at higher luminances an unknown non-linear mechanism diminishing the 
low-frequency response comes into play, thus causing the sensitivity 
maxima around 10cps in DE LanGn’s [4] other curves. The latter statement 
is justified by inspection of a plot of Dn LAanGn’s [4] measurements with 
log L and ML? as axes and the frequency as are (L stands for 
luminance). Statistically optimal behaviour is represented by straight lines 
parallel to the log L-axis; deviations from this are largely confined to the 
points for high luminances and low frequencies at the same time. 

The said non-linear mechanism must have a complicated character, in 
view of the fact that low-frequency behaviour is strongly influenced by 
the luminance of the surrounds [12, 8]. It is presumably responsible for 
lack of fit in Rosk’s [9] plots, as well as for Broca—Suuzmr’s [13]. 
Britcke’s [14] and related phenomena. Part of it might be attributed to 
or at least described as, a lengthening of the recovery time after firing a 
rod or cone as a function of the firing rate. Such a lengthening might be 
brought about by a bottle-neck in the energy supply necessary for 
rebuilding an electrical double layer. 

Returning to the linear diffusion mechanism discussed here, a final 
remark may be made about the boundary condition cz-, =0 stating that 
stimulus-molecules are inactivated on reaching the wall. 

From the solutions of the diffusion equation given above, the response 
R; to a 6-function at the input can be ee by means of Laplace- 


ret f 
transforms. It is found to be R;=e-%- — - dy’ m (0%) where 7? stands for 


the Jacobi #-function, or again, using #1’ (0, 1/k)=k?/2 dy" (0, k), 


i 27 \ 3/2 Die 
Safi, eo wae | W.== N\. 
HS Ar @) v1 ( a 


F 20 
Using #41’ (0, &) = 2(r—3r8+ 575 ...) with i= Coe Sard, rt and 


considering that for t> 0.25 sec. e-%<2.10-%, it is seen that for the deter- 


98 Series B 
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mined values of 9 and rt only the first term of the series need be considered 
to obtain a 1 per cent-approximation. 
The result is 


Qt 
Ry & . $-3/2. g—et. p—T/2t. 


A 


This is a well-known solution for the diffusion equation, which is obtained 
by neglecting the boundary condition ez, =0. 

This, of course, means that so few of the stimulus-molecules produced 
reach the wall that it does not matter much whether or not they are 
inactivated too. 
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PALEONTOLOGY 


MIOGYPSINA IN HUNGARY 


BY 


C. W. DROOGER 


(Communicated by Prof. G. H. R. von KoENIGSWALD at the meeting of Febr. 25, 1961) 


Abstract. The discovery of Miogypsina septentrionalis in northern Hungary, 
together with other foraminiferal species more or less characteristic of the Chattian, 
justifies the assumption of some direct way of marine connections between Hungary 
and northern Germany during the Late Oligocene. 


INTRODUCTION 


In order to clarify the position of the Chattian-Aquitanian boundary 
in their country, Hungarian paleontologists started a detailed investigation 
of a continuous sedimentary series near Novaj, some 10 km from Eger 
in northern Hungary (fig. 1). This boundary has been a constant problem 
in this part of Europe, as may readily be seen from the wealth of papers 
on this subject, for instance those of Fucus (1894), Szérs (1956) and 
CSEPREGHY—MEZNERICS (1960). 

The lower part of the investigated series consists of sandy marls, the 
microfauna of which is a Clavulinoides szaboi association. Such associations 
are generally considered to be of Middle Oligocene (Rupelian) age. 
Overlying these marls, glauconitic, clayey sands occur, which in their 
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Fig. 1. Sketch map of Hungary showing the location of Novaj and Eger (after 
R. M. Nyrr6). 
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higher part contain a level with indurated limestone lenses with Lepido- 
eyclinids. Still higher in these sands, towards the top of the unit, a sample 
from the vineyard Nyarjas was found to contain a number of Miogypsinids. 
On top of the glauconitic sands the series continues with a grey marine 
clay (for instance that of the brickworks Wind), which in turn is overlain 
by sandy deposits, partly of continental origin. With rhyolitic tuffs, 
covering these sands, the series of Novaj comes to an end (fig. 2). 
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Vig. 2. Schematic section of the sedimentary series near Novaj (after T. BAwop1). 


The arrow indicates the position of the Miogypsina-bearing sample. 


Rhyolitie tuffs. 


Grey, limonitie sandstone with plant remains. 


Coarse sand with gravel, micaceous. 


Grey, more o1 less limonitie clay with molluses (Wind brickworks). 


Glauconitic, clayey sands with limestone lenses. 


Be Rae S, ts a . yo « * a 
RSS Sandy marls with Clavulinoides szab6éi association (Rupelian). 
MN ac SS 


A detailed publication on the stratigraphy and fauna of these deposits 
near Kger is to be expected in the near future in the Annals of the 
National Museum of Natural History of Hungary (T. BAwpr, a.o., 1961). 
The stratigraphy and the molluscan fauna will be treated by TamAs 
3ALDI, the smaller Foraminifera by Rika M. Nyrré, and the Lepido- 
cyclinids by Trsor Kecskemért. The author feels gratefully indebted 
to these Hungarian colleagues, and especially to Dr. Nyrré, for sending 
the Miogypsinids to Utrecht, and for kindly supplying him with strati- 
graphic and faunistie details for this publication. 


THE MIOGYPSINIDS 


Altogether some fifty specimens of Miogypsina were available. Most 
of them are small, consisting of only the embryonic and nepionic stages 
with but a few equatorial chambers. 
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Bigger specimens are roughly fan-shaped with a broadly rounded 
apical portion. These individuals are commonly unequally biconvex, in 
some instances one side being completely flat. The observed diameter of 
the test ranges from 0.7 mm to 2.4 mm, the corresponding thickness from 
0.3 mm to 0.9 mm. The greatest thickness occurs at the place of the 
embryonic-nepionic stage. Pustules on the surface are distinct. In larger 
specimens they are fairly regularly distributed with the bigger ones over 
the thickest portion of the test. Their size range is from 35 to 200 wD, 

Transverse sections of bigger specimens, five in total, show a nepionic 
spiral which is only slightly trochoid. The median layer was found to be 
strongly variable in height: 125-250 wu (internal diameter). This is mainly 
on account of the inflated character of the median chambers themselves. 
The side walls contain distinct lateral chambers, relatively few in number, 
and in section strongly different in size and shape. Their arrangement 
is variable, commonly very irregular and not in distinct layers or tiers. 
Their number between median layer and the surface of the test varies 
from place to place, from zero to four. The horizontal walls, especially 
those of the median chambers, are fairly thick and with distinct pores 
(fig. 3, 4). 

For the first approximation of the specific characters twelve successful 
horizontal sections of macrospheric specimens could be used. They all 
show embryonic chambers of about equal size and a simple nepionic 
spiral of about 14 coil (fig. 5-7). The equatorial chambers are generally 
relatively few in number. A maximum of six rows of these were found. 
Their shape is arcuate to ogival, and the greatest observed dimensions 
are 175250 uw. In some specimens the median layer has been damaged 
by burrowing organisms. 


Fig. 3-7. Miogypsina septentrionalis Droocmr. x 36. 3, 4, transverse sections ; 
5-7 median sections. Upper Glauconitic Sands, vineyaid Nyarjas near Nova), 
Hungary. 
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The features of the embryonic-nepionic apparatus of these twelve 
individuals are as follows: 

Diameter protoconch: Rr=95-160 uw; Mri=119 uw; om,=4.8 up. 

Diameter deuteroconch: Ry=90-140 uw; Mn=118 wu; Om, =3.9 m- 
Mr/My =().99. 

The angle between the apical-frontal line and the line connecting the 
centres of both embryonic chambers in the section: R,= —340°— — 180°: 
M,= —266°; om,=12.1°. | 

The number of nepionic chambers: Rx = 22-14; Mx=17.8; OM, = 0.6. 

The number of operculinid nepionic chambers (after which distinct 
equatorial chambers begin to develop) is in many specimens hard to 
decide. A second, distal stolon was evidently often present already in 
chambers adjoining the apical border of the test, each one giving rise 
to but a single, usually very small equatorial chamber. The probable 
number of operculinid chambers ranges from 9 to 3.M=5.4:; om—0.4. 

In the sections the nepionic chambers usually increase fairly regularly 
in size, and the decrease afterwards is not abrupt either. An average 
value for the largest nepionic chamber had to be chosen in several of the 
specimens. This increase in size of the nepionic chambers occurs up to 
the 9th to 14th chamber (mean value 12.3; o4=0.6). A number of three 
to seven chambers follow this largest chamber, before the spiral comes 
to an end. 

Furthermore, some seven sections of specimens from the same locality 
had been made at Budapest. As far as their features are sufficiently clear, 
six of them are within the range of variation indicated above. The 
seventh, a slightly oblique section, is clearly beyond these limits with: 
diameter I and II both 200 yw; y= —70°; X=10: only one operculinid 
chamber, and increase in size up to the ninth chamber. It is unlikely 
that this specimen could have been reworked from older deposits, since 
its features represent rather those of a more advanced type. As yet it 
had better be regarded as an extreme variant of the population at hand. 

Including all available sections, also two which were made later and 
again fall within the primarily found range, and also that of the individual 
with X= 10 (total N=18 to 20), we obtain mean values for the enumerated 
features which differ by about one standard error from those given above. 
All combined values give shifts towards the individual values of the 
aberrant specimen. 

A single microspheric individual was succesfully sectioned. Both its 
protoconch and deuteroconch are about 15 f# in diameter. X =27; the 
number of operculinid chambers amounts to 14: increase in size occurs 
up to the 23d or 24th nepionic chamber: there are only two rows of 
mainly ogival equatorial chambers. 

Finally we had at our disposal two specimens from the claypit of the 
brickworks Wind, near Eger, and derived from the basal part of the 
grey marine clay, which overlies the glauconitic sands. The sections are 
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hot very good, but the number of nepionic chambers, the angle y, and 
other visible features are comparable to those of the majority of the 
individuals from Nydrjas. For this reason these few specimens could 
have been reworked from the glauconitic sands, but for us there is no 
need for such an explanation. According to Dr. Nyiré reworking is 
unlikely because of the filling of the specimens by pyrite instead of by 
glauconite. 

As to the specific determination of the Nyarjas assemblage, the out- 
standing point is the peculiar combination of a very high average number 
of nepionic chambers with the presence of distinct, though often poorly 
developed lateral chambers. Very similar assemblages were recently 
described from Chattian deposits in northwestern Germany as Miogypsina 
(Miogypsina) septentrionalis DRoocER (1960, p. 41, pl. 1, fig. 1, 2, pl. 2, 
fig. 1-8). Comparison of the Hungarian and German assemblages shows 
distinct differences. Some of these differences are statistically very 
significant (diameter protoconch), others somewhat less, of the order of 
2 o to 3 a, such as those for Mx, M,, and the mean number of operculinid 
chambers. In all these features the Hungarian assemblage is less advanced. 
The average diameters of I and II are smaller, the mean value of the 
angle y is greater, on the average the nepionic spiral and its operculinid 
part are longer, and the maximum size of the nepionic chambers occurs 
slightly later. Regarding the principle of nepionic acceleration the 
Hungarian assemblage could be slightly older than those from Astrup 
and the Doberg in Germany. Just as it was shown for the latter, the 
Nyarjas population must be considered as closely related to the populations 
of Miogypsina (Miogypsinoides) complanata SCHLUMBERGER, which 
occurred in the more southern parts of the Mediterranean region. 

Although statistically significant differences could be shown to exist 
between the Hungarian and the German assemblages, known so far, 
their exceptional combination of the very high Mx values and the 
beginning development of a system of lateral chambers, makes it certainly 
preferable to group them together in one species: M. septentrionalis. 


SOME OTHER FAUNAL ELEMENTS 


The microfauna accompanying the Miogypsinids at Novaj, could be 
evaluated from part of the original sample, kindly sent by Dr. Nyrro. 
Without entering into much detail (see R. M. Nyrr6 and T. KecskeMsrt, 
in press), some of the species may be mentioned, which enable us to give 
an insight into the association, and which are considered to have some 
stratigraphic or paleogeographic value. Such value being considered 
doubtful, we did not try to determine specifically the accompanying 
species of Amphistegina, Operculina and Heterostegina. They merely stress 
the conclusion on environment that can be drawn from the Miogypsina 
species, which is open marine and at shallow depth, probably not 


exceeding 50 meters. 
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The few complete Lepidocyclina specimens of the sample are of moderate 
size; diameter 0.7—3.5 mm. Their central part is inflated and there is 
a thin flange which has usually been damaged. The centre has one or more 
pustules of considerable size (up to 150 «); the remainder of the surface 
has very tiny pustules, it is smooth or it shows a meshwork of the lateral 
chambers’ walls. Horizontal sections show that the individuals belong to 
the European Nephrolepidina lineage. Several specimens have been 
damaged by burrowing organisms. In five of the median sections the 
number of accessory auxiliary chambers on the deuteroconch could be 
counted (fig. 8, 9). In two individuals three such chambers were found, 
in the three others there are four. These numbers are considered to be 
primary features for specific determination. They would fit in best with 
the diagnosis of Lepidocyclina morgani Lemoine and R. Douvitik (1904, 
Soc. géol. France, Paléont., Mém. 32, p. 17, pl. 1, fig. 12, pl. 2, fig. 4, 
pl. 3, fig. 2; DRooGeR and Socry, 1959, Micropaleont., vol. 5, p. 420). 


Fig. 8,9. Embryonie-nepionie stage in median seetion of two specimens of Lepido- 
cyclina morgani LEMOINE and R. Dovuvitie. & 48. Upper Glauconitic Sands, 
vineyard Nydrjas near Novaj, Hungary. 


The shape of the equatorial chambers is ogival to short spatulate- 
hexagonal. Their arrangement in concentric circles is not distinct at all 
ontogenetic stages. Several individuals show irregular growth, demon- 
strated by the intercalation of rings of too small chambers. Especially 
among the later equatorial chambers the diagonal stolons are not exactly 
in the same plane, though their distal ends often coalesce. Annular basal 
stolons were observed, but not as a regular feature. VAN DER VLERK 
(1959) and GRimMsDALE (1959) in their recent studies of the Lepidocyclinidae 
stress the importance of the stolon pattern. Our species would belong to 
the Y lineage. Calculating the evolution factors, as advanced by VAN DER 
VLERK (1959), we find a range of 35 to 44 for factor A. of 44 to 65 for 
factor B, and of 84 to 100 (average 91) for the total. These figures are 
again best comparable with those of L. morgani. 


Some flattened fragments in the sample evidently belonged to very 
large Lepidocyclina specimens. Their surface shows a thickened meshwork 
of the lateral chambers’ walls. In sections the shape of the equatorial 
chambers (fig. 10) is long-spatulate. The chambers are arranged in 
concentric rows. The thin diagonal stolons are clearly crossed and there 
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are equally distinct annular basal stolons. The thinness of the test and 
especially the shape of the equatorial chambers indicate that we are 
dealing with a species of the Hulepidina series. Probably we may determine 
these fragments as belonging to Lepidocyclina dilatata (MicHELorTTt). 
GRIMSDALE (1952) has shown that most of the stated differences between 
described Hulepidina species are invalid for specific determination. 
Whether the European L. dilatata can really be differentiated from 
GRIMSDALE’s L. ephippioides (JONES and CHAPMAN) must be considered 
doubtful, if one accepts the wide variation he gives to the latter species. 


atte 


Fig. 10. Horizontal section of part of the median layer of Lepidocyclina dilatata 
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(MicHELorti), showing crossed diagonal stolons and basal annular stolons between 
the equatorial chambers. * 48. Upper Glauconitic Sands, vineyard Nyarjas near 


Novaj, Hungary. 


Among the smaller benthonic Foraminifera may be mentioned the 
presence of Almaena osnabrugensis (ROEMER). This species was originally 
described from the Chattian deposits of western Germany. It has been 
reported from other European areas as well, and it is said to range into 
the Aquitanian, and in Central Europe even into the Burdigalian (CrcHa, 
1960). This upward extension of the range cannot be doubted, but so 
far the species has not been found together with Miogypsina species of 
Aquitanian or Burdigalian age. 


Yet another interesting benthonic species of the Nyarjas sample is 
Bulimina kasselensis Batses (1958, p. 127, pl. 5, fig. 4-6), also originally 
described from the German Chattian deposits. In all their details the 
Hungarian individuals are identical with Bavsus’s specimens. Never- 
theless, not too much stratigraphic value can be given to this species, 
its range not yet being fully known. Similar forms occur in the Miocene 
of Germany (B. dingdenensis Batsus, 1958, p. 128, pl. 5, fig. 1-3) and in 
the Miocene of northern Italy. The latter may be closely related to the 
species T'ritubulogenerina dertonensis GiaNort1 (1953, Riv. Ital. Pal. 
Strat., Mem. 6, p. 271, pl. 12, fig. 7) from the type Tortonian. Among the 
individuals from the Miocene of these countries, stored in the collections 
of the University of Utrecht, especially some of the smaller ones are not 
clearly different from the Chattian specimens. 


In the Nyarjas sample representatives of the Globigeri nidae are frequent, 
but mainly small (up to 200 y). All belong to the genus Globigerina and 
probably no more than a single biological species is represented. The 
majority of the specimens are referable to Globigerina globularis (ROEMER) 
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(DRooGER and BatyEs, 1959, p. 174, pl. 1, fig. 2). Adjacent types seem 
to be present, but their determination had to be based on very few 
specimens: Globigerina cf. suteri (Bouut), G. ef. venezuelana HEDBERG, 
G. obesa (BOLL), G. wnicava (Bouur, LorspiicH and Tappan) and possibly 
G. ciperoensis Bouut (Bout, 1957). This association resembles those of 
the German Chattian deposits, though it might be claimed to resemble 
those of other deposits of different age as well. The only distinct con- 
clusion is that the character of the planktonic association, consisting of 
Globigerina types only, would better fit in with an Oligocene than with a 
Miocene age of the Nyarjas sample. 


CONCLUSIONS 


Considering the age of the fauna from Nydarjas, the joint occurrence of 
Miogypsina septentrionalis and Lepidocyclina morgani certainly supplies 
most of the evidence, since both species belong to rather well known 
evolutionary series. They have been found together in the deposits of 
the German localities Astrup and Doberg. The latter sediments are 
irrevocably of Chattian age. 

The evolution of the European series of Nephrolepidina is not yet 
sufficiently known. As a consequence the occurrence of L. morgant does 
not give the restriction of a Chattian age. 

However, the evolution pattern of the Miogypsinidae has been somewhat 
better established. Although the assemblage of M. septentrionalis of 
Nyarjas is considered less advanced than the German ones, the difference 
in age is probably slight. No arguments can be found in the faunal 
composition to suppose a Rupelian age. Also in the North Sea basin the 
difference between the faunas of the Rupelian and the Chattian deposits 
is strongly influenced by the change in facies, which makes the geologic 
time factor and the Rupelian-Chattian boundary practically elusive. 
On the other hand the Miogypsinids preclude a younger, Aquitanian age 
of the Nyarjas sediments. The equivalent of septentrionalis in the 
Aquitaine basin is M. complanata, which occurs in distinctly pre-Aqui- 
tanian sediments of this region. In the typical Aquitanian it had been 
replaced by the more highly developed species of the same stock: 
M. gunteri and M. tani. It cannot be accepted as probable that the 
Hungarian stock of M. septentrionalis remained unchanged for a long 
time, reaching into the Aquitanian. The occurrence of a population in 
Slovakia, which is morphologically intermediate between the two Aqui- 
tanian species M. gunteri and M. tani (Papp, 1960; CrcHa, 1960), makes 
it very likely that also in the Hungarian—Slovakian area the Miogypsinidae 
had a development, comparable to that in other regions. 

The conclusion of a Chattian age 1s supported with more or less weight 
by the encountered Globigerina globularis association, and by the presence 
of Almaena osnabrugensis and Bulimina kasselensis. As has been shown, 
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none of these additional elements gives conclusive evidence, but they 
strengthen the comparison with the Chattian deposits of Germany. 

The presence of an Hulepidina species at Nyarjas opens up possibilities 
for comparison in southern direction. In the pre-Aquitanian deposits of 
southwestern France (Saint-Géours, Saint-Etienne d’Orthe, CuCl 
ephipproides occurs together with L. morgani, M. complanata and A. 
osnabrugensis, only the last species of which has been claimed recently 
to be still present in younger deposits of this area. In itself Hulepidina 
(L. dilatata or L. ephippioides) is not a proof of Chattian age. Represen- 
tatives have been found in Early Miocene deposits of the Mediterranean area. 

About a year ago it could be concluded from the discovery of Miogypsina 
and Lepidocyclina in the Chattian deposits of Germany, that at least 
one wave of these larger Foraminifera had succeeded in penetrating from 
the Mediterranean area into the North Sea basin. One of the groups, the 
Miogypsinids, was considered to have developed a somewhat aberrant 
northern offshoot, M. septentrionalis, but still well comparable to the 
Mediterranean stock of M. complanata. The assemblages of the latter 
species from Christus (SW. France) and from the Bric del Duca (N. Italy) 
were mentioned as closely comparable, though already slightly more 
highly developed in their nepionic features. One may not expect, not 
even in an area such as the Mediterranean, in a narrow sense, that the 
change from Miogypsinoides into Miogypsina s.str. took place at but 
one level of nepionic development. Local deviations certainly occurred. 
Nevertheless, it is striking that the German M. septentrionalis, which in 
this respect is highly aberrant, is found again, and still more deviating 
in Hungary. Of course, one might assume that the German and the 
Hungarian M. septentrionalis populations had developed completely 
independently of one another from a similar stock, but then one might 
expect differences in minor features. The fact that they nicely fit into 
a single evolving lineage, added to the consideration that the separate 
rise of so very similar aberrant populations not very far from one another 
looks very unlikely, makes it much more plausible that both population 
groups were closely related. 

This conclusion immediately requires some idea about the marine 
north-south connections during the Late Oligocene. The more so, if not 
only L. morgani, but also A. osnabrugensis and B. kasselensis of the 
German Chattian are considered to be immigrants from the south. The 
more or less classical solution, which assumes a connection from the 
Kassel area and the Mayence basin, along the Rhine depression and the 
Swiss Molasse basin towards the Lower Rhone area, has lost much of 
its attraction. One reason is that in the Swiss Molasse foredeep non- 
marine conditions prevailed during the Chattian (age based on mammals), 
resulting in the “Untere Siisswassermolasse”. This blocks marine con- 
nections to the south in between the zones of the present Alps and Jura 


Mountains. 
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It is well known that the Oligocene marine faunas of Upper Bavaria 
very closely resemble those of Hungary (Haan and Hoéuzu, 1952), and a 
direct marine connection between these two areas may be considered as 
well established. Admittedly, the distance from Bavaria to Kassel and 
Mayence is not great, but it is difficult to bridge this gap, since marine 
sediments of Oligocene age are not continuous. They are at the least very 
indistinct. One might think of some short-time marine connection from 
Upper Bavaria in the direction of Basel and from there northwards 
along the Rhine depression and Alsace, but the indications are not very 
clear (see for instance OERTLI, 1956). Part of the fully marine sediments 
of this strip may still be present, buried under younger deposits, but 
another part of them has probably been removed by later erosion. It is 
certainly worth while to compare again the faunas of both ends to see 
whether such a connection could more easily explain the geographic 
distribution of still other species. 

Whether the marine connections from Hungary to the more southern 
Mediterranean area during the Late Oligocene went directly south cannot 
be decided either. Miogypsinids are known to occur, together with 
L. morgant, in northern Yugoslavia (Papp, 1954, 1955). They belong to 
M. complanata s.1., but their nepionic stage is more highly developed 
than that of M. septentrionalis. Further evidence has to be awaited. 
Recalling the discoveries of Miogypsinids in Slovakia, Austria (Papp, 1960) 
Hungary and Germany during the last few years, it becomes evident 
that this group was fairly well distributed in and along the northern 
parts of the Alpine area. Further discoveries may therefore certainly be 
expected. 
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1. Introduction 


“Voleano-tectonic” depressions are a much discussed subject of volcanic 
geomorphology. The name of these, often rather irregularly shaped, 
depressions indicates that they are supposedly the result of both volcanic 
and tectonic phenomena. They are consequently considered as an inter- 
mediate category between the more symmetrical calderas formed by 
pure-volcanic collapse on the one hand, and the normal tectonic features 
—graben, down warped basins and the like—on the other. The leading 
thought is that blowing off or outflow of the magma results in the collapse 
of the roof of the magma chamber, which process results in the formation 
of graben and downwarped basins. The volume of the issued voleanic 
products generally is much less than the volume of the depression 
concerned and the latter then apparently are partly due to other, tectonic, 
causes (V. BEMMELEN, 1949, p. 212). The relative importance of both 
these factors should be carefully studied for every individual case. The 
tendency exists, in the author’s opinion, to emphasize especially the role 
of volcanism in the collapse mechanism. 

Many “‘voleano-tectonic” depressions have been the scene catastrophic 
eruptions of the Katmaian type. Fiery tuff-flows issued from fissures on 
those occasions and the huge amounts of deposited ignimbrites now form 
extensive plateaux, first described from New Zealand by MarsHALu ( 1935). 


1) Geomorphologist, International Training Centre for Aerial Survey, Delft 
The Netherlands. Formerly: Geographical Institute, Indonesian Topographical 
Survey, Djakarta. 
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The island of Sumatra is the classical area of these “voleano-tectonic” 
depressions. A great number of larger and smaller features of this type 
exist, the largest one being the Lake Toba depression. They were described 
as such by vAN BEMMELEN (1929, 1939, 1949) who introduced the 
“emptying out-and-caving in’ hypothesis. He also emphasized the 
existence of associated “volcanu-tectonic’” uparching of the earth’ crust, 
caused by the rising magma of a near-surface batholith. The so-called 
“Batak-tumor”’, since is known as the largest and clearest example of 
this phenomenon. The Lake Toba ‘“‘voleano-tectonic” depression was 
formed in a later stage, according to VAN BrmMELEN. The vast ignimbrite 
sheets of Sumatra have been described by WrsTERVELD (1947, 1953). 

The author, on behalf of the Geographical Institute, Djakarta, carried 
out four larger geomorphological field investigations in the Barisan 
Mountains of Sumatra in 1954 and 1955, the results of which will be 
published soon. Special attention was paid automatically to the “voleano- 
tectonic’ depressions, because of their conspicuous occurrence throughout 
the island. The Kerintji, Singkarak, and Toba depressions were studied 
in greater detail and the author’s opinions therefore, are based primarily 
on his observations in these areas. The general conclusions regarding their 
origin, included in this paper, deviate considerably from the established 
ideas on the nature of these depressions. They may prove of some value 
for the study of similar features elsewhere. Unfortunately, the coverage 
of aerial photographs is negligible for the areas under consideration, and 
the study is based almost exclusively on field observations. 


2. Lake Kerintji 


The Kerintji plain forms part of the longitudinal graben of the Barisan 
Mountains in the southernmost part of Central Sumatra. The area was 
visited by the author in 1954 (VERSTAPPEN, 1955). It is a wedge-shaped 
complex graben lying at an altitude of approximately 800 metres above 
sea level. The plain is bordered by steep fault-scarps, having a height of 
about 400 metres and converging towards the north west. Consequently 
the narrowest and highest part of the graben is located there. The plain 
becomes broader and lower lying towards the south-east and it reaches 
its maximum width of about nine kilometres in the extreme south-eastern 
end where the “volcano-tectonic” lake of the same name, having a 
maximum depth of 110 metres, is situated. The narrow north western 
apex of the graben becomes more conspicuous due to the occurrence of 
two splinters in the northern part of the graben, which converge even. 
more than the two main faults. The graben finally terminates at the 
North-Kerintji granite and quartz diorite massif, where the further 
continuation of the faults can be ascertained only with difficulty in the 
field. The impression gained is that the North Kerintji massif influenced 
the younger tectonic events of the region, which is reflected in the poor 


430 


in ii 
aint ite 


a Ae tron 


Fig. 1. Geomorphic map of the wedge shaped Kerintji graben. 


Key: 1. Areas adjacent to the graben. 2. Splinteis in the graben. — 3. Tuff 
cone of lake Kerintji. 4. Rhyolite plug. 5. Fault-searps. 6. Hot springs. - 


S = Sungaipenuh. 


development of the northern part of the Kerintji graben. Similar situations 
were observed in several other parts of the Barisan Mountains. 

The Kerintji “volcano-tectonic” depression is only a comparatively 
small one. The volume of the tuffs issued can be roughly estimated at 
4 km, and the volume of the depression, as deduced from a lake survey 
carried out by the author in 1954, amounts to approximately 1,75 km3. 
It should be emphasized that the tuffs were deposited in the graben. 
The remainders of the tuff cone which originally surrounded the lake, 
are still present in the plain and form hills of considerable size, for instance 
the Kotapetai hill. Thus it is evident that the eruption took place in the 
pre-existing Kerintji graben. The drainage of the graben was not altered 
by the volcanic event: the main river enters the depression in its narrow 
north western end, and breaks through the eastern fault-scarp in the 
extreme south east. It leaves the graben nowadays at exactly the same 
locality as before the eruption and it has incised its course there in 
Kerintji tuffs deposited in its pre-eruption valley. The depression appar- 
ently is almost exclusively of a tectonic nature. No “‘voleano-tectonic”’ 
uparching can be traced either. 

The river was dammed temporarily by the ejected tuffs and thus the 
graben was flooded. The water rose to about 90 metres above the present 
lake level, but it was lowered rapidly due to the quick incision of the 
river in the tuffs, and the plain ran dry again. A few terrace levels now 
indicate the subsequent stages of the lowering in lake level. The volcanic 
eruption and the flooding of the plain are registered in local folk tales 
and thus the event seems to have occurred in historic times. Another 
volcanic feature is the occurrence, along the eastern fault-scarp, of a 
rhyolite plug, a tholoid, which also was formed after the graben came 
into being. Earthquakes and young scarplets point to recent tectonic 
movements along the existing fault lines. 
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3. Lake Singkarak 

This characteristically steep sided and flat bottomed “volcano-tectonic”’ 
lake is about 20 kilometres in length and its maximum width is 7 km. 
It is located in the Central Sumatran section of the Barisan median 
graben, approximately 200 kilometres to the north-west of Lake Kerintji. 
It lies at an altitude of 362 metres above sea level and reaches a maximum 
depth of 268 metres. The relevant sheets of the topographic map of 
Central Sumatra, reyeal that this lake evidently occurs also in the very 
lowest part of this section of the graben. 

The Solok plain, located further to the south-east, is only slightly 
higher lying than the level of Lake Singkarak (362 m). Its altitude ranges 
from 387 metres above sea level in the north-west to about 450 metres 
in the south-east. The greater height is more conspicuous, however, if 
the depth of Lake Singkarak (268 m) is taken into consideration: the 
altitude of the Solok plain above the lake bottom is between 293 and 
356 metres. The floor of the graben rises rather abruptly to the south-east 
of the Solok plain and reaches a height of approximately 1,500 metres 
above sea-level in the surroundings of the lakes Baruh and Atas, at 
some 35-40 kilometres distance from the southern shore of Lake Singkarak. 
The rise of the graben floor to the north-west of the lake is even more 
rapid. A height of 1,180 metres above sea level is reached at a distance 
of 19 kilometres from the north western end of the lake, where between 
the Merapi and Singgalang volcanoes, a broad saddle’ is formed in the 
young volcanic products. 

This saddle is important as a graben-watershed. The terrain to the 
north of it gradually slopes down towards the ignimbrite plateau of 
Bukittinggi, about 900 metres above sea level. The southward flowing 
rivers once entered Lake Singkarak. They were captured, however, in 
comparatively recent times, by the Anai river and since drain towards 
the Indian Ocean in the west. The pirate Anai river is deeply incised in 
a narrow gorge. The capture occurred near the town of Padang Pandjang, 
and the present, almost imperceptible divide between the Anai branches 
and the beheaded lower river courses still flowing towards the lake, lies 
at about 770 metres above sea level. 

The rivers of the remainder of the graben section just described, 
including the Solok plain and as far south as Lake Baruh, flow towards 
Lake Singkarak, which is drained eastward by the Umbilin river. This 
situation only came into being rather recently, however. VAN VALKENBURG 
(1921) mentions the existence of a dry river valley to the east of the 
town of Solok. Evidently the Solok plain once drained eastwards via this 
nowadays deserted valley and the Lasi-Pumuatan river, which river 
debouches in the Umbilin river, many kilometres further downstream. 
The height of the dry valley is about 410 metres above sea level and 
equals the height of the higher terraces at the Umbilin outlet of Lake 
Singkarak. Obviously, both passages have been in use simultaneously. 
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The present north-west drainage of the Solok plain towards Lake Singkarak 
can only be explained by assuming a recent subsidence of the graben 
section between Solok and the lake. Several geomorphic features evidence 
the extremely young faulting movements occurring in this particular 
part of the graben structure. 


Umbilin river 
(terraces) 


Lake Singkarak 


Fig. 2. Longitudinal profile of the Singkarak section of the median graben. 


The geomorphic development sketched above is essentially due to 
tectonic causes. Some voleanic activity was affiliated to the faulting 
movements, however. Young, whitish, volcanic tuffs are present in the 
graben to the north of Solok. They are not related to any young volcanic 
cone, but they originated from a fissure eruption at the eastern side of 
the graben, which accompanied the recent collapse of this section of the 
median graben, to which reference was made above. The tectonic move- 
ments continued after the fissure eruption and a number of well-preserved 
young scarplets can be traced in the tuffs. The Sumani river locally 
breaks through the tuffs in the central and eastern parts of the graben, 
even though an easy course through the flat and lower lying westerly 
zone seems indicated. This points to a continued tectonic subsidence, 
witnessed by the western side of this critical graben section. 

The fault bordering the eastern side of Lake Singkarak also acted as a 
volcanic vent. The tuffs were deposited to the east of the fault-scarp 
and they form a surface with a gentle easterly slope (fig. 3) between the 
Umbilin river and Mt. Batu Baragung. Another, smaller flow occurs to 
the south of this mountain. A third tuff deposit is a tuff cone in the graben 
and enclosing the northern side of the lake (fig. 3). 

The deposition of the tuffs in the graben to the north and south of 
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Fig. 3. Geomorphic map of the Singkarak section of the median graben. 
Key: 1, Areas adjacent to the graben. — 2. Tuff deposits (partly in slightly higher 
graben sections). — 3. Plains. — 4. Voleanic cones. — 5. Fault-searps. — 6. Dry 
river valley. — M = Merapi — S = Singgalang — T = Talang — B = Batubara- 

gung — P = Padang Pandjang — Si = Singkarak — So = Solok 
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the lake justifies a conclusion that this graben also existed prior to the 
voleanic events. The eastward directed tuff flows occurring to the east 
of the lake similarly indicate the pre-existing relief. It is beyond doubt 
that the volcanic (fissure) eruptions, geomorphologically speaking, were 
of comparatively little importance. They only occasionally accompany 
the intense faulting movements, provided that a batholith is present at 
shallow depth. They are definitely not the cause of the depression. 

The Kerintji and Singkarak depressions evidently do not fit van 
BEMMELEN’S definition of ‘“volcano-tectonic depressions”. Either this 
definition should be modified into “tectonic depressions with occasional 
volcanic activity along faults”, or the term should not be applied at all. 
VAN BEMMELEN (1949, p. 684) lists 7 “volcano-tectonic” depressions in 
southern Sumatra, an area he studied in more detail. The formation by 
volcano-tectonic collapse of two of these depressions, viz. the Suoh and 
Gedongsurian depressions, is rejected by WESTERVELD (1953). He explains 
the Gedongsurian depression by tectonic causes with subsequent erosion, 
and considers the Suoh depression as a graben where acid lavas occur 
along fault-zones. WESTERVELD’s views on these features are rather 
similar to those developed by this author concerning the Kerintji and 
Singkarak depressions in Central Sumatra. The question arises whether 
volcano-tectonic uparching and collapse, in the sense of vAN BEMMELEN, 
exists at all in Sumatra. In seeking an answer to this question, it seems 
relevant to study the “‘classical example’’ of this type of phenomenon 
viz. Lake Toba, located in Northern Sumatra. 


4. Lake Toba 
A. Critical considerations about the supposed “Batak-tumor” 


VAN BEMMELEN claims that the Toba cauldron is formed on top of a 
huge “tumor” formed by volcano-tectonic uparching and measuring 
275150 km. He states that the outline of the “tumor” is marked by 
the contour lines of 100 and 1,000 metres. The basement complex in the 
top part of the “tumor” reaches a height of more than 2,000 metres 
above sea-level as can be derived from the topographic maps (vAN 
BEMMELEN, 1949, la, p. 687). 

However, a few objections are to be raised against this conception. 
The contour lines, to start with, are without significance to outline the 
assumed “‘tumor’’ because they indicate the slope of the ignimbrite 
plateau and not the pre-volcanic surface. It is established by borings that 
near the town of Tebingtinggi, some 70 kilometres from the lake, the 
ignimbrite sheet is still 80 metres thick (v. Dp. Maren, 1948). Closer to 
the lake a thickness of several hundreds of metres is reached. If one 
takes the thickness of the ignimbrite sheet into consideration, the contour 
lines do not indicate an exceptional width of the Toba section of the 
Barisan Mountains. One could say that the parts adjacent to the Toba 
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area are exceptionally narrow. The width of the Toba section may be 
slightly above the average, but this is readily explained by the somewhat 
complicated structure of this part of the Barisan Mountains viz. the 
occurrence of the Toba graben to the east of the main longitudinal graben. 

Furthermore, the height of the basement complex in the Toba area is 
not an extreme value. Similar heights are reached further to the south 
in the Barisan Mountains near Panjabungan and Hutanopan (Mt. Sopo 
Untjin, 2,199 m). Further to the North heights of 3,000 metres and more 
occur (Mt. Bandahara, 3,030 m; Mt. Leuser, 3,381 m). 

The present topography certainly does not appear to justify the 
assumption of a “Batak-tumor’. Warped lake terraces, formed in the 
Toba tuffs and occurring along the eastern shore of Lake Toba, were 
studied by the author (p. 12). It is evident from these geomorphic 
observations that after the Toba eruption(s) parts of the area were further 
uplifted, locally as much as 250-350 metres. Immediately after the 
eruption the topography thus was even less in line with the ‘“‘tumor’’- 
hypothesis. 

The post-voleanic relief actually is not at all crucial in this relation, 
but only the topography before the ignimbrite eruptions counts. It does 
not seem imperative that a “tumor” should still prevail after the 
“emptying-out-and-collapse”” phase. VAN BEMMELEN emphasizes on the 
one hand the existence at present of the Batak ‘‘tumor’’, whereas he 
states on the other that the “tumor” of the much larger Strait Sunda 
dome, cannot be traced nowadays. He claims that in this latter case 
“the field of gravitational stresses was too strong for the bearing-power 
of the crust, so that it collapsed” (vAN BEMMELEN, 1949, la, p. 633), 
without any proof, however. The evidence of the Lake Ranau ‘‘tumor’’, 
Southern Sumatra, is not very convincing. The coastal features of the 
nearby part of the Indian Ocean, mentioned in this connection (VAN 
BEMMELEN, 1949, la, p. 680) have no certain relationship to the Lake 
Xanau eruption, neither in space nor in time. 

The original, pre-voleanic topography of the Toba-area is, of course, 
very hard to establish due to the ignimbrite cover and because of the 
important tectonic events that occurred. The existence of a radial drainage 
pattern at that time, seems likely, if a dome or “tumor” has been present. 
It was observed by the author that the drainage of the Karo Highlands, 
situated to the northwest of Lake Toba prior to the eruption, was directed 
towards the lake. A radial drainage evidently did not exist there. It is, 
in view of the above considerations, rather unlikely that the assumed 
Batak “tumor”, which cannot be traced at present, has ever existed at all. 


B. The Toba graben 


This “voleano-tectonic” lake, the largest one of Sumatra, measures 
approximately 87 kilometres in a northwest-southeast direction and its 
maximum width js about 26 kilometres. The graben itself is still somewhat 
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Fig. 4. Tectonic sketch of the wedge shaped Toba graben. 
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Key: 1. Hot springs. — 2. Fault-scarps within the main graben. — 3. Areas ad- 
jacent to the graben. — 4. Blocks/splinters within the main graben. 
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Fig. 5. Geomorphie map of the SE part of the Toba graben. 


Key: 1. Pre-Toba rocks emerging from the tuff cover. — 2. Areas emerging from 
lst Toba terrace. — 3. Ist Toba terrace. — 4. 2nd Toba terrace. — 5. Multi-cycle 
alluvial fan. — 6. Fault-scarp. — 7. Overflow of Toba tuffs. — 8. Escarpment of 


Asahan amphitheatre. 
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larger, viz. 10031 km and it is bordered by steep fault-scarps, many 
hundreds of metres high. The island of Samosir is an important feature 
in the lake which finds its counterpart in the Sibulangit or Uluan peninsula 
located between Prapat and Porsea. The greatest depth of the lake is 
529 metres and as the present lake level is at 906 metres; the deepest 
part of the graben lies at 377 metres above sea level. 

The origin of the Toba trough is a much discussed problem, which 
centres on the question whether it was a graben existing prior to the 
Toba eruption, or whether it originated after the volcanic events by 
voleano-tectonic collapse. Of the older investigators Wrxa EastTon 
(1894, 1896) is in favour of the graben hypothesis. WESTERVELD (1947, 
1953) is among the modern scientists holding similar views. He maintains, 
with good right, in the present author’s opinion, that a ‘‘shallow acid 
magma-chamber underlying a pre-existing tectonic depression, of which 
the bordering fault-fissures furnished the eruption channels” accounts 
for the phenomena observed. 

Vorz (1899, 1909) dates the graben as a post-eruption phenomenon. 
Van BEMMELEN (1929, 1939, 1949) advocates the voleano-tectonic origin 
of the trough, which assumption also implies that the trough basically 
was formed after the eruption. Van BemMELEN, however. is not very 
consistent about this matter. He states (v. BEMMELEN, 1949. la. p. 212) 
that “the Lake Toba area was probably already a graben, drained by 
the Asahan river, before the eruption took place”. He then continues 
“it is difficult, if not impossible, to estimate which part of the volume 
of the present Toba-depression belongs to the pre-Toba graben, and 
which part is due to the over-deepening by the eruption’. He concludes: 
“Because the Toba cauldron was formed partly by the doming up of the 
Batak tumor, which caused tensional stresses in the top area and the 
subsidence of a graben, and partly by the blowing out of the contents 
of the underlying shallow-seated, granite batholith, it may be called a 
“volcano-tectonie depression”. VAN BrmMELEN thus admits that the 
graben, at least in part, existed prior to the eruption. He expresses a 
quite contradictory view on p. 691, however: “No investigator who has 
set eyes on the imposing precipitous escarpments surrounding the Toba 
cauldron, will ever contend that this trough existed already before the 
outburst”. Also, no pre-existing graben is indicated in his schematic cross 
sections (p. 695) of the various stages of development of the Toba area, 

It should be emphasized that Lake Toba actually is located in a wedge 
Shaped graben having its apex in the north-west. T he southern end of 
the lake is formed by an E-W directed escarpment, which is most likely 
another fault-scarp. The graben is drained at its broadest southern end 
by the east-ward flowing Asahan river to the Malacca Straits. WestERVELD 
(1947) pointed out that the Asahan already drained the graben before 
the Toba eruption occurred, a view shared by the present author. Even 
the oldest products of the Toba voleanic vents penetrate far downstream 
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in the pre-existing valley carved mainly in pre-tertiary rocks (FRYING, H.., 
1925). These deposits consequently now form a tongue shaped valley fill 
where the river re-incised its course in a narrow gorge. The older, con- 
solidated, tuffs at present can be traced as far as 18 km to the east of 
the eastern Toba fault-scarp. 

The fact that the wedge-shaped Toba graben always drained at its 
broadest south-eastern end indicates that this area was also, most 
probably, the deepest part of the graben. It should be mentioned in this 
context that at least the more recent Toba volcanic activity centered 
mainly in the SE (v. D. Marg1, 1947, 1948a). There is evidently a certain 
degree of similarity between the Kerintji and the Toba graben structures: 
Both are wedge-shaped features with their apex pointing towards the 
NW and drained eastward at their broadest and lowest lying south 
eastern parts, where volcanic events occurred. This similarity was even 
more pronounced formerly, when the Karo highlands, located to the NE 
of the lake, still drained towards lake Toba in much the same way as 
the Siulak river at present still enters the Kerintji graben at its NE 
extremity. The Karo river system was beheaded by a branch of the 
Renun river system, after being pounded up by Toba tuffs and as a 
result of the increased vertical erosion of the Renun river system due to 
further subsidence of the Renun graben. The problem will be discussed 
at some length in a forthcoming publication of the author on the 
geomorphology of Sumatra. 

Van BEMMELEN (1949, la, p. 689, 691) also is of the opinion that the 
Asahan river already existed before the Toba eruption. Because of his 
views that the formation of the Toba trough mainly occurred after the 
volcanic events, it is inevitable for him to consider the Asahan valley as 
a relic of a river system that became beheaded by the downbreak of the 
Toba cauldron. The headwaters of this ‘‘proto-Asahan’”’, according to him, 
were located in the western part of the present Lake Toba in between 
the Renun and Batang Toru catchment areas. The great thickness of the 
Toba tuffs deposited in the Asahan valley would not appear to support 
this hypothesis. 

These tuffs reach a height of 1,160 metres above sea level at the 
entrance of the Asahan gorge to the east of Porsea. This tuff terrace 
slopes down downstream and is at 750 metres above sea level near the 
well-known Tangga falls at approximately 15 km from the eastern Toba 
fault scarp. The oldest, consolidated tuffs outcrop, according to the map 
of Frynine (1925) near the base of the Tangga falls. This outcrop is 
at a height of about 250 metres above sea-level and the total thickness 
of the tuffs at that locality consequently amounts to at least 500 metres. 
The bedding of the tuffs is approximately parallel to its surface and one 
may therefore safely assume that the thickness of the tuffs at closer 
range to the lake is at least the same and possibly more. The height 
above sea level of the Asahan valley bottom underneath the tuffs thus 


438 


amounts to 1,160—500= 660 metres or less, near Porsea. It will appear 
from the author’s observations on the terraces around Lake Toba, that 
possibly the Porsea area has been uplifted some 110 metres after the 
eruption occurred. The height of the Asahan valley bottom immediately 
after the eruption then was only 550 m above sea level at the most. The 
greatest depth in the southern part of the lake occurs to the north of 
Balige and amounts to 428 metres (vAN BrEMMELEN, 1939), or in other 
words 906—428=478 metres above sea level. Evidently this part of the 
lake bottom is only slightly lower-lying than the valley bottom of the 
pre-eruption Asahan river. One may object that the present lake bottom 
is not representative due to deposition of Toba tuffs in the graben and 
as a result of post-eruption tectonic movements. It is more justified 
therefore to state that the original Asahan valley bottom is located 
approximately 1,500 metres below the highest parts of the edge of the 
eastern Toba fault-scarp. Collapse of the Toba graben as a result of — and 
thus after—the eruption, is obviously of a negligible order, if it existed 
at all. It is beyond doubt that the Asahan river, before the eruption, 
drained the already existing Toba graben. 

Voz (1899, 1909) based his views of a post-eruption age of the Toba 
graben on the consideration that otherwise the graben would have been 
filled up completely with ejected tuffs. A detailed knowledge of the 
Katmaian type of eruption is necessary to account for the limited deposition 
of tuffs in the graben. It is suggested by the author that the — comparatively 
small amounts of—tuffs occurring on the bottom of each of the three 
“voleano-tectonic’’ depressions described in this paper, were deposited 
during the less violent phases of eruption, whereas the tuffs were deposited 
in the higher lying surroundings during the paroxysmal parts of the 
eruption. 

A study of the pre-existing relief alone is not sufficient to estimate 
the influence of voleano-tectonies in the Toba area. The crustal movements 
which occurred posterior to the eruption also should be considered. Such 
young tectonic movements are rather important, as appears from the 
author’s observations, and thus even less space is left for the assumed 
voleano-tectonic events at the end of the eruption period(s). 

The Toba graben is a complex feature. A few faults of secondary 
importance and splinters occur along the western lake shore. The most 
striking features, however, are the south-westerly tilted Samosir block 
and the Sibulangit (Uluan) block, which are separated from each other 
by the Latung Strait. The fault-scarp bordering the Sibulangit block to 
the west can be traced in a southerly and south-westerly direction across 
the Porsea Bay and to the west of Balige. The Porsea Bay itself, in the 
author’s opinion, is to be considered as a drowned part of the original 
Asahan valley. 

RurrneEr (1935) reported the occurrence of lacustrine (diatom) deposits 
on the island of Samosir as high as 1,360 metres above sea-level and—as 
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he leaves young tectonic movements out of consideration—he claims 
that the lake formerly had a level about 500 metres higher than nowadays 
and that it was drained southward by the Batang Toru river. The 1,400 m 
contour does not close around Lake Toba, a fact which led to RurtNER’s 
assumption of a regional tilt. WEsTERVELD (1947), who shares Rurrner’s 
opinion about the 1,400 metres lake level, tries to overcome part of the 
difficulty by supposing a considerable retrogressive erosion in the tuffs. 
However, this is geomorphologically not justified, as only narrow gorges 
are formed in the tuffs which cannot account for the broad gaps in the 
1,400 m contour. A former southerly drainage of the lake seems to be out 
of the question, after a thorough field study of the author. No dry river 
valley/gorge can be traced there. 

VAN BEMMELEN (1939) certainly is correct in assuming that the highest 
lake level is indicated by the upper tuff terrace near Porsea. This terrace 
is at present located at 1,160 metres above sea-level. The author agrees 
with vAN BEMMELEN that only young tectonic movements can explain 
the high location of the former lake deposits. V. p. Maren (1948a.) 
maintains that the higher lying deposits are formed in small ponds and 
not in Lake Toba itself. The author in 1955 found an extensive lake 
deposit in eastern Samosir at not less than 1,600 metres above sea-level, 
which certainly was not formed in a small pond. The deposit contained, 
apart from sponge needles, a rather large number of the following diatoms: 


Melosira distans (Khr.) Kiitz var. robusta Manguin 
Melosira granulata (Ehr.) Ralfs var. valida Hustedt 
Cyclotella stelligera Cl. & Grun. 


This combination points to a littoral-pelagic freshwater deposit, f.i. 
along the shores of a large lake !). The author is thus inclined to reject 
v. D. Marev’s hypothesis, which opinion is further strengthened by a 
study of the lake terraces. 

The Toba lake terraces are mentioned by several authors, but they are 
never studied in detail. The terraces are not preserved everywhere around 
the lake because the surf action in the lake was not vigourous enough to 
form terraces where hard rocks occur, and furthermore because slumping 
—a common phenomenon in the diatom deposits — destroyed the terraces 
in other localities. Recent emergence due to young tectonic movements 
is another possible cause for the scarcity of terrace remnants (Samosir f.i.). 
The best preserved terraces occur in the extensive tuff deposits in the 
south eastern part of the graben. These were investigated by the author, 
The study of the Toba terraces is complicated by the frequent occurrence 
of sub-horizontal rocks, where structural (pseudo) terraces are formed, 


1) Determination by A. v. p. Werrr, bio-geologist, Geological Survey, Haarlem, 
The Netherlands on the request of Prof. Dr. J. WEsTERVELD, Amsterdam. The 
author is greatly in debt for their kind cooperation. 
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which are without any significance for our purpose. Many older observ- 
ations therefore had to be rejected, as being incorrect. 

Two main terrace levels were distinguished. The highest one is 
represented by the 1,160 metres level at the entrance of the Asahan 
valley near Porsea. This tuff terrace corresponds with the upper level of 
the tuff terraces found along this river (p. 10). This terrace (I) was 
observed almost uninterrupted along the eastern fault-scarp to the south 
of Prapat and some scattered remains were found along the southern 
side of the graben to the east of Balige. It is also developed on the 
Sibulangit peninsula, where only a few small areas emerge from it, 
representing former islands. A second important terrace (II) is best 
developed along the road Prapat-Porsea, between the remnants of terrace I 
on the peninsula and along the eastern fault-scarp. Its height, near 
Porsea, is 1,050 metres. The strong tilt of these terraces can easily be 
seen in the field, even with the naked eye. Important post-eruption 
tectonic movements, resulting in terrace warping, occurred between 
Prapat and Porsea, as indicated in fig. 6. The height of the 1,160 m 
terrace varies from 1,070 to 1,420 metres. It is possible that this terrace 
—and thus the maximum lake level—was originally about 1,050 metres 
above sea-level and that the Porsea section was afterward uplifted some 
110 metres (p. 438), thus forcing the Asahan river to increased vertical 
erosion. The uparching demonstrated by the terraces evidently occurred 
after the eruption and thus can not be explained by a ‘‘tumor-hypothesis” 
(except one connected to a future eruption). Small remnants of a number 
of lower terraces also occur, which have their counterpart in the multi- 
cycle alluvial fans formed in the lower parts of the tuff deposits. 

[t is important to note that the Toba lake terraces can be correlated 
to structural terraces formed in the Toba tuffs/ignimbrites along the 


son A: Terrace I along easterly fou/tscarp. 


B:Terrace I on Sibvlangit peninsula. 
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Fig. 6. Warped lake ter:aces along the eastern shore of lake Toba between Prapat 
and Porsea. 
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Asahan valley. The intermittent lowering of the lake level is readily 
explained by the different rate of incision of the Asahan river in the 
various subsequent volcanic strata. Every hard bed in the valley 
corresponds with a lake terrace. Two more lake terraces will be formed 
when the Siguragura (Wilhelmina) and Tangga falls of the Asahan river, 
which are both due to resistant beds. will have reached the lake as a 
result of retrogressive erosion. The stratification of these older Toba 
products points to the fact that a number of eruptions have occurred 
instead of the one or two (v. p. MAREL) assumed so far. 

The important post-eruption tectonic events on and near the Sibulangit 
peninsula make it likely that the nearby island of Samosir witnessed a 
similar development. No terraces can be traced there, however. This 
can only in part be due to destruction of such terraces by the slumping 
phenomena in the diatom deposits. The author is inclined to assume 
that the island was almost completely submerged after the eruption 
and only emerged afterward. The fault/flexure bordering Samosir in the 
east, the general westward tilt of even the youngest deposits and the 
diatom deposit at 1,600 metres altitude, point in the same direction. 

The northern part of Lake Toba was probably rather stable after the 
eruption, as no terrace warping was observed there. The young crustal 
movements evidently centre in the Sibulangit peninsula and on Samosir. 
Extremely young faults were observed by the author in the north-western 
part of Samosir. This area may, at present, well be the most unstable 
part of the Toba graben. 

The above considerations lead to the conclusion that the development 
of the Toba graben is basically tectonic and is spread over a considerable 
period. It is accompanied by occasional violent volcanic outbursts. It 
does not appear to be justified to attribute the depression or cauldron 
to volcano-tectonic processes in the sense of VAN BEMMELEN’s “‘emptying- 
out-and-breaking down”’ hypothesis. 


5. Summary 


The three ‘“‘volcano-tectonic” lakes studied, viz. the Lakes Kerintji, 
Singkarak and Toba, are evidently formed in a pre-existing graben 
structure. The longitudinal graben of the Barisan Mountains is more or 
less like a guirlande, with “‘volcano-tectonic depressions” formed mostly 
in the lowest and broadest parts. The latter is probably due to the fact 
that the faults there also are best developed and reach the greatest depth. 
Thus they have a better chance to reach a near-surface batholith, if such 
a batholith is present underneath the area concerned. No “‘volcano- 
tectonic uparching” prior to the eruption, was observed, but in the case 
of Lake Toba important young tectonic movements could be proved. 

The view here expressed is just the opposite of vAN BEMMELEN’s 
‘““emptying-out-and-caving-in” hypothesis on the origin of such depressions 
according to which the cauldrons are formed as a result of and thus after 
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a volcanic eruption, occurring preferably in a culmination or “tumor” 
of the mountains. The term ‘‘voleano-tectonic depression”? should be 
abandoned according to the author, or its definition should be altered 
into “tectonic depression with occasional volcanic activity along faults”. 
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